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Abstract: The Majorana nature of neutrinos can be experimentally verified only via 
lepton-number violating processes involving charged leptons. We study 36 lepton-number 
violating (LV) processes from the decays of tau leptons and pseudoscalar mesons. These 
decays are absent in the Standard Model but, in presence of Majorana neutrinos in the 
mass range ~ 100 MeV to 5 GeV, the rates for these processes would be enhanced due to 
their resonant contribution. We calculate the transition rates and branching fractions and 
compare them to the current bounds from direct experimental searches for AL = 2 tau and 
rare meson decays. The experimental non-observation of such LV processes places stringent 
bounds on the Majorana neutrino mass and mixing and we summarize the existing limits. 
We also extend the search to hadron collider experiments. We find that, at the Tevatron 
with 8 fb~^ integrated luminosity, there could be 2a {5a) sensitivity for resonant production 
of a Majorana neutrino in the jJi^jJir^ modes in the mass range of ~ 10 — 180 GeV (10 — 
120 GeV). This reach can be extended to ~ 10-375 GeV (10-250 GeV) at the LHC of 14 
TeV with 100 fb"^ The production cross section at the LHC of 10 TeV is also presented 
for comparison. We study the //^e^ modes as well and find that the signal could be large 
enough even taking into account the current bound from neutrinoless double-beta decay. 
The signal from the gauge boson fusion channel W^W'^ ^i^t LHC is found to be 

very weak given the rather small mixing parameters. We comment on the search strategy 
when a r lepton is involved in the final state. 
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1. Introduction 

In the Standard Model (SM) of strong and electroweak interactions, neutrinos are strictly 
massless due to the absence of right-handed chiral states {Nn) and the requirement of 
SU{2)l gauge invariance and renormalizability. Recent neutrino oscillation experiments 
have conclusively shown that neutrinos are massive j^] . This discovery presents a pressing 
need to consider physics beyond the SM. It is straightforward to obtain a Dirac mass 
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term mDipiNji + h.c.) for a neutrino by including the right-handed state, just hke the 
treatment for ah other fermions via Yukawa coupUngs to the Higgs doublet in the SM. 
However, a profound question arises: Since Nr is a SM gauge singlet, why should a gauge- 
invariant Majorana mass term ^MNrNr not exist in the theory? In fact, there is strong 
theoretical motivation for the Majorana mass term to exist since it could naturally explain 
the smallness of the observed neutrino masses via the so-called "see-saw" mechanism [01 



(1.1) 

From a model-building point of view, there are many scenarios that could incorporate the 
Majorana mass. Examples include Left-Right symmetric gauge theories |^]; S'O(IO) Su- 
persymmetric (SUSY) grand unification Q and other grand unified theories models 
with exotic Higgs representations [^, ^; R-parity violating interactions (AL = 1) in Super- 
symmetry (SUSY) § and theories with extra dimensions |^]. There are other proposals 
to generate Majorana masses for neutrinos at a higher scale M without relying on the 



right-handed state [10, 11 1. According to the scheme in generating the mass scale M 



in Eq. (|1.1| ), it has been customary to call them Type I Q, Type H [|T^ or Type HI [pT]]. 



Within the context of the SM, there is only one gauge-invariant operator [12| that is 
relevant to the neutrino mass, 

IlH, (1.2) 

where and H are the SM lepton and Higgs doublets, respectively. The constant k is a 
model-dependent effective coupling and A is the new physics cut-off scale. It is a dimension- 
5 non-renormalizable operator, and leads to Majorana neutrino masses of the order kv'^ / A, 
after the Higgs field acquires a vacuum expectation value v, in accordance with the see-saw 
scheme. Higher dimensional operators that give rise to Majorana neutrino masses have also 
been constructed in a model-independent manner [^]. The challenging task is to look for 
experimental evidence to probe the new physics scale A and to distinguish the underlying 
theoretical models mentioned above. 

In the neutrino sector, besides the rich phenomena of neutrino flavor oscillations and 
the possible existence of new sources of CP- violation, lepton number violation by two units 
(AL = 2), as implied by a Majorana mass term, plays a crucial role. Not only may it 
result in important consequences in particle physics, nuclear physics and cosmology but 
it would also guide us in understanding the fundamental symmetries of physics beyond 
the SM. Although the prevailing theoretical prejudice prefers Majorana neutrinos, exper- 
imentally testing the nature of neutrinos and lepton-number violation {LV) in general, is 
of fundamental importance. In accelerator-based experiments, neutrinos in the final state 
are undetectable by the detectors, leading to the so-called "missing energy" and therefore 
missing lepton numbers as well. One is thus forced to look for charged leptons in the final 
state. The basic process with AL = 2 can be generically expressed by 

W-W~ ^ £^q, (1.3) 

where W~ is a virtual SM weak boson and £1^2 = e, /U, r. By coupling fermion currents 
to the W bosons as depicted in Fig. [l|, and arranging the initial and final states properly, 
one finds various physical processes that can be experimentally searched for. The best 
known example is neutrinoless double-beta decay (Oz^/3/3) [|l^, |l^, |l^, which proceeds via 
the parton-level subprocess dd — > uu W~*W~* — > uu e~e~. Other interesting classes 
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Figure 1: A generic diagram for AL = 2 processes via Majorana neutrino exchange. 



of LV processes involve tau decays such as t~ — > £~^M^M2 
mesons Mi,M2 are Tr,K, rare meson decays such as 
hyperon decays such as S~ — > T,~^e~e~, E~ — > pfi~ ^~ etc. [21]. One could also explore 
additional processes like e" 



1 13, where the light 
H, 10, [H and 



[|35| , charged-Higgs production e 




neutrino 
top- 
6|, and in 



, /i — > e"*" ||18|, and fi ^ conversion 
One may also consider searching for signals at accelerator and collider experiments via 
e'e- W~W- [H], e+e- ^ N + X 

nucleon scattering vi(j7i)N PtflfX 
quark decays t 
the decay N ^ i^H^ ||]. 

The dynamics for AL = 2 processes as in Eq. ( |1.3D is dictated by the properties of the 
exchanged neutrinos. For a Majorana neutrino that is light compared to the energy scale 
in the process, the transition rates for LV processes are proportional to the product of two 
flavor mixing matrix elements among the light neutrinos and a LV mass insertion 



I2 



m=l 



(1.4) 



where {m)g_^g^ are the "effective neutrino masses" . If the neutrinos are heavy compared to 
the energy scale involved, then the contribution scales as 



3+n 



(1.5) 



where V is the mixing matrix between the light flavor and heavy neutrinos. Unfortunately, 
both situations encounter a severe suppression either due to the small neutrino mass like 
m1^/M^r, or due to the small mixing | Vgjm'V^jm'P- An important observation is that 
when the heavy neutrino mass is kinematically accessible, a process may undergo resonant 
production of the heavy neutrino. The transition rate can be substantially enhanced and 
goes like 

T{N^, ^ i) T{N^, ^ /) 



m 



N„ 



N 



(1.6) 
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where i, f refer to the initial and final state during the transition. 

The possible existence of sterile neutrinos in the mass range relevant for resonant en- 
hancement of AL = 2 processes studied in this paper is motivated in several scenarios. 
Models which implement the see-saw mechanism at low energies have been recently con- 
sidered the neutrino masses generated are accidentally small and active-sterile 
mixing can be as large as few percent. See-saw models at the electroweak scale can explain 
neutrino masses if appropriate symmetries are imposed and at the same time provide an 
appealing mechanism for baryon asymmetry generation via resonant leptogenesis In 
theories with dynamical electroweak symmetry breaking, sterile neutrinos with masses in 
the 100s of MeV to GeV range are invoked to explain light neutrino masses |41]. Sterile 
neutrinos can also play a role in understanding the flavour problem in the leptonic sector. 
It has been shown that mixing with sterile neutrinos can be at the origin of the large angles 
in the neutrino sector |42]. 

Heavy, mostly-sterile neutrinos have been investigated for their role in cosmology 
and astrophysics, in particular in Big Bang Nucleosynthesis, Large Scale Structure for- 
mation [43 1, cosmic microwave background, diffuse extragalactic background radiation, 
supernovae [44| and as dark matter candidates |4^, (for a review on MeV ster- 
ile neutrinos, see Ref. [48|). A keV sterile neutrino is a viable dark matter candidate 
1 45, which can also explain the origin of pulsar kicks [49|. Decays of heavy, mostly- 
sterile neutrinos have been proposed to explain the early ionization of the Universe |^0|| . 
Due to mixing, dark matter sterile neutrinos would decay radiatively contributing to the 
Diffuse Extragalactic Background Radiation and inducing x-ray emission from galaxy clus- 
ters [51, 52 1 . A large coupling between sterile neutrinos and light dark matter scalars can 
be at the origin of neutrino masses and of the observed dark matter abundance [^]. A 
model with sterile neutrinos in the keV-GeV mass range has been proposed to explain the 
dark matter of the Universe as well as baryogenesis [54, 55 [. Its phenomenological and as- 
trophysical signatures have been considered in detail in Refs. [56, 52[. This model assumes 
the existence of one sterile neutrino with keV mass for dark matter and two heavier neu- 
trinos with quasi-degenerate GeV masses for successful baryogenesis. The required mixing 
\Vim'\'^ of the latter neutrinos with the active ones is mass dependent and lies in the range 
10^^^ — 10~^, for a mass of 1 GeV. Additional constraints on the heavy neutrino mass and 
mixing angles can be derived from astrophysical observations. Sterile neutrinos mixed with 
active ones would be efficiently produced in supernovae cores, escaping from it and deplet- 
ing substantially the supernova core energy, and, therefore, might modify the supernova 
evolution. Recently, it was shown that sterile neutrinos in the mass range ~ 0.2 GeV and 
small mixing angle with /i and r neutrinos could enhance the energy transport from the 
core to the stalled shock and favor the supernova explosion They could also explain 
the high velocity of pulsars if the momentum carried away by heavy sterile neutrinos is 
emitted asymmetrically [^^. Detailed reviews and discussions of heavy neutrinos in the 
Early Universe and their present bounds can be found, e.g., in Refs. 



I Mm- 

Cosmological and astrophysical constraints on sterile neutrinos are typically very strong 
but are not as robust as the ones from laboratory searches as they typically depend on the 
production mechanism of sterile neutrinos in the Early Universe and on the cosmological 
evolution. For example, they can be significantly weakened or evaded if the reheating 
temperature is low [47, |5^, if their density in the Early Universe is diluted by entropy 
injected at late times [^] or if they have non-standard interactions. In these cases, much 
larger mixing angles with active neutrinos are allowed by cosmological observations and 
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can be tested in terrestrial experiments. Therefore, it is important to perform experimental 
searches of heavy sterile neutrinos with increased sensitivity and, specifically for Majorana 
neutrinos, to consider AL = 2 processes. If a positive signal is found and is incompati- 
ble with the cosmological and/or astrophysical observations, one would need to consider 
modifications to the standard cosmological scenario and/or would gain new insight on the 
evolution of astrophysical objects. 

In this paper, we study resonant contributions of heavy Majorana neutrinos to AL = 2 
processes involving two charged leptons in accelerator-based experiments. We establish 
our conventions and discuss the current constraints on the mass and mixing of heavy 
neutrinos in Sec. ^. In Section ^ we lay out the general expressions for the heavy neutrino 
contributions to low energy LV decays and study two classes of AL = 2 processes, 

(a) tau decays r~ ^ M2 , 

(b) rare meson decays K+, D+, Df, B+ ^ M2 ■ 

We calculate the enhanced transition rates and branching fractions and compare them 
to the bounds set by direct experimental searches. A non-observation of such AL = 2 
processes places stringent constraints on the mass and mixing of Majorana neutrinos which 
are also presented in this section. The resonant production of Majorana neutrinos at hadron 
colliders, namely the Tevatron and LHC are studied and updated in Sec. |^. We draw our 
conclusions in Sec. We discuss in detail the formalism, the decay modes and the total 
decay width of heavy Majorana neutrinos and the transition rates of LV processes in the 
Appendices. 

2. Majorana neutrinos in extension of the standard model 

To set up our notation and convention, we first discuss the formalism for the simplest 
extension of the SM which includes right handed singlets. Also in this section, we present 
the current constraints on the mass and mixing of a heavy neutrino from various direct 
detection experiments, accelerator searches and electroweak precision constraints. 

2.1 Formalism for Heavy Neutrino Mixing 

The leptonic content in our simplest extension of the SM includes three generations of 
left-handed SM SU{2)l doublets and n right-handed SM singlets 



LaL=\^{'j , NbR, (2.1) 

where a = 1,2,3 and b = l,2,---,n. The gauge-invariant Yukawa interactions lead to 
Dirac masses for the charged leptons and neutrinos after the Higgs field develops a vacuum 
expectation value v. It is also possible for the singlet neutrinos to have a heavy Majorana 
mass term. The full neutrino mass terms as well as the diagonalized eigenvalues can be 
expressed as 

= ^ I E E <b Nm + ML <: <r) +T. Bw Nb,R 1 + h.c. 

b,b'=l J 
3+n \ 

^ v^i^-^^+Yl Mn^, N„,r + h.c. (2.2) 

\m=l m'=4 / 
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with the mixing relations between the gauge and mass eigenstates 



3 3+n 

^ Vam'N'^,^, (2.3) 

m=l m'=A 

UU^ + VV^ = I. (2.4) 
In terms of the mass eigenstates, the gauge interaction Lagrangian can be written as 

(r 3 r 3+n \ 

J2 E t^^'Pl^ + E E N^rPLi + h.c. 
e=e m=l i=e m'=4 / 

(r 3 T 3+n \ 

E E ^^rn ^I^Pl + E E ^T'^^L ^£ + h.c.(2.5) 

£=e m=l ^=e m'=4 / 



g_ 

2 cos 6*^/ 



Further details about the mixing formalism are given in Appendix 

A few important remarks are in order before the detailed considerations. First of all, 
parameterically, the light neutrino masses rriy^ diag are of the order of mag nitude {m^^f/B, 
while the heavy neutrino masses are Mjv^ diag — B. Secondly, the mixing parameters would 
typically scale as C/Tf/ ~ / and V^V ^ m^/MN. Thirdly, the Majorana mass term for the 
flavor states Uai, absent in Eq. ( ^I^ ) and corresponding to the null entry 03x3 in Eq. ([A.^ ), 
may receive non-zero contributions as Majorana masses for the light active neutrinos, for 
instance from higher dimensional AL = 2 operators or in theories with a triplet Higgs 
field. The general formalism presented here remains the same. In this paper, we will take a 
phenomenological approach toward these parameters. We will simply take the masses and 
mixing elements of the heavy neutrino as free parameters, only subject to some constraints 
from experimental observations. The assumption that the masses and mixing elements 
are not rigorously related by the see-saw relations is feasible from a model-building point 
of view, since some fine-tuning or some ansatz of the neutrino mass matrix can always 
alter the general relations. Several scenarios where it is possible to have rather low mass 
of the heavy neutrino were mentioned in the previous section. Here and henceforth, we 
consider the case when only one heavy Majorana neutrino is kinematically accessible and 
denote it by with the corresponding mass and mixing with charged lepton flavors 
Vi4. If we stick with this simple parameterization, the SM Higgs boson will couple to the 
heavy neutrinos as well. We present the couplings in Appendix |A[ When appropriate, we 
will include this effect. As noted above, some fine-tuning would be needed to avoid 
excessive contributions to the light neutrino mass. 

2.2 Current Constraints on N4 Masses And Mixing 

In laboratory searches, no positive evidence of sterile neutrinos has been found so far,^ in 
the mass range of interest, 100 eV-100 GeV.^ 



^Indications of the existence of a neutrino with 17 keV mass were subsequently shown to be non valid. 
For a review, see Ref. |6l| . Studying interactions of neutrinos from it and fi decays, an anomaly in time 
distribution was found [p2|. It could be interpreted as the existence of a neutrino emitted in pion decays 
with mass of 33.9 MeV. Searches for this neutral fermion have not given any positive signature and have 
allowed to constrain the mixing to be |V^4p < 9.2 x 10"® at 95% C.L. 

^For sterile neutrinos with smaller masses a rather complete analysis of the bounds can be found in 
Ref. . See also the implications of the recent results from the MiniBooNE collaboration ^ . 
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A very powerful probe of the mixing of heavy neutrinos with both and z/^ are peak 



searches in leptonic decays of pions and kaons |67]. If a heavy neutrino is produced in 



such decays, the lepton spectrum would show a monochromatic line at 



2mM 



where and mi are respectively the lepton energy and mass, niM is the meson mass. The 
mixing angle controls the branching ratio of this process as: 



T{M+^l+Ni) |y,4|' 



V^mIV , (2.7) 



where p \s a. kinematical factor 



1 + /if + /i| - 2(/i^ + /U4 + pi^i) [fie + IJ.4- - PaY 
P = ^ 7, 7^ ^' (2.^ 

with Hi = mf/mj^. For large m^, the helicity suppression of the tt,K ^ iug decays weakens 
and there is an enhancement for M"*" — > l'^ N/^ by a relative factor ml/mf, reaching up to 
10^ — 10^ compared to that of vr — > eue and K eu^ in the SM, respectively. These bounds 
are very robust because they rely only on the assumption that a heavy neutrino exists and 
mixes with and/or i/^. 

Another strategy to constrain heavy neutrinos mixed with i/g, and v-r, is via searches 
of the products of their decays. If kinematically allowed, A^4 would be produced in every 
process in which active neutrinos are emitted with a branching fraction proportional to 
the mixing parameter |V£4p. They would subsequently decay via Charged Current (CC) 
and Neutral Current (NC) interactions into neutrinos and other "visible" particles, such 
as electrons, muons and pions. Searches for these "visible" decay-products were performed 
and were used to constrain the mixing parameters. In beam dump experiments, A^4 would 
be produced in meson decays, with the detector located far away from the production 
site. The suppression of the flux of A^4 needs to be taken into account if the decay length 
is very short and, therefore, typically both an upper and a lower bound on the mixing 
angle can be set. Otherwise, the production can happen in the detector itself, as for the 
limits obtained from a reanalysis of LEP data, using the possible decays of the [26| 
into heavy neutrinos. In this case, large values of the mixing angle can be excluded. 
These bounds are less robust than the ones previously discussed. In fact, if the heavy 
neutrinos have other dominant decay channels into invisible particles, these bounds would 
be weakened, if not completely evaded. For example, a coupling of the type gNvcf) (see 



Ref. |68, |5^), with (/> a scalar, can induce very fast decays, which might dominate over the 
ones induced by CC and NC interactions. In this case, if the decay length is very short due 
to these strong interactions, the flux of A'4 might be suppressed at the far detector and the 
bound would not apply. If the production happens in the detector itself, the bounds would 
need to be recomputed considering the branching fraction into "visible" channels. Notice 
that we do not report here the bounds from Ref. as they do not apply to the heavy 
neutrinos under consideration. In these analyses it was assumed that heavy neutrinos were 
produced via — > N4^N4^ with the same strength as an active neutrino. In our scheme, 
this would correspond to having mixing angle equal to 1. Then, the search for A^4 decays 
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in the detector was used to constrain the heavy neutrino parameters. These data should 
be reanalyzed considering that the production of A'4 is suppressed by |V£4p. Comparing 
the expected number of events with the backgrounds we estimate that typically bounds of 
order |V£4p < few x 10~^-10~^ could be deduced. However, a detailed analysis should be 
performed and we do not report these limits in our figures. 

For masses above the production threshold, additional constraints can be obtained 
from lepton universality as the decay rates for muons, pions, taus as well as the invisible 
decay width for the boson are modified with respect to the SM predictions ||7^, |72| . 
Flavour changing neutral current processes such as /i — > 67, ee'^e~ and ^u-e conversion 
in nuclei are affected by the existence of heavy sterile neutrinos and strong limits can be 
obtained on the mixing with active neutrinos [^, 74, 75]. These bounds are reported in 
Section p. 2.4 



Finally, in Section 2.2.5 we discuss the very strong constraints on |T44p which can be 



obtained from the non-observation of neutrinoless double beta decay. It should be noticed 
that in the presence of more than one sterile neutrino, possible cancellations between the 
contributions to the decay rate can be achieved and the bounds would be consequently 
much weaker. 

Next, we review the laboratory constraints on the mixing between heavy and active 
neutrinos, depending on flavour and the mass of sterile neutrinos. 

2.2.1 Mixing with i/g 

The mixing parameter Ve4 can be tested in searches of kinks in the /3-decay spectrum, of 
peaks in the spectrum of electrons in meson decays and, finally, of iV4 decays in reactor 
and accelerator neutrino experiments. 

For masses 30 eV ~ 7714 ~ 1 MeV, the most sensitive probe is the search for kinks in 
the /3-decay spectra In the presence of heavy neutrinos mixed with Ug, the Kurie plot 
would be given by the contributions of the decays into light neutrinos as well as into heavy 
ones. This induces a kink in the Kurie plot at the end point electron energy Eg 

2Mi ' ^ ' 

where Mj j are the mass of the initial and final nuclei, respectively, and mg is the electron 
mass. In Fig. |2| we report the most stringent limits, obtained by using different nuclei 
1 76, 77, 78, |8^. In reactors and in the Sun only low mass, mi < few MeV, heavy 



neutrinos can be produced. The constraints obtained by looking for their decays into 
electron-positron pairs |8^ are reported in Fig. ^ with solid (cyan) contour labeled 
Bugey and short dashed (blue) contour labeled Borexino. The region with long dash 
dotted (grey) contour, labelled vr eu, is excluded by peak searches 



For heavier masses peak searches give the most stringent bounds, shown in Fig. |3|. 
Notice that, due to the weakened helicity suppression of the tt decay, the sensitivity on 
Vei increases with m4 till phase space becomes relevant at 7714 > 80 MeV, for vr efg- 
The excluded region, at 90% C.L., from Ref. [^ ], is indicated with the solid (black) line 
labeled tt ev. For heavier masses, stringent bounds are obtained by looking at the 
electron spectrum in K decays |8^] and are indicated by the double dash dotted (purple) 
line labeled K ^ ev \w Fig. Assuming that only CC and NC interactions are at play, 
stringent bounds have been obtained on |V^4p and are reported in Fig. ^ by the rest of the 
contours (except dotted (maroon) line labeled Oz^ /?/?). In particular, the limits at 90% C.L. 
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Figure 2: Bounds on |Ve4p versus 7714 in the mass range 10 eV-10 MeV. The excluded regions 



with contours labeled ^^^Re |76|, , ^^Ni |7| , Q , ^Op and Fermia |80[1 refer to the 

bounds from kink searches. All the limits are given at 95% C.L. except for the ones from Ref. ||8^ 
which are at 90% C.L.. The areas delimited by short dashed (blue) contour labeled Borexino and 
solid (cyan) contour labeled Bugey are excluded at 90% C.L. by searches of A'4 decays from the 
Borexino Counting Test facility |81 and Ref. |Q respectively. The region with long-dash-dotted 
(grey) contour, labelled tt eu, is excluded by peak searches Q. The dotted (maroon) line labeled 
Oi^PP indicates the bound from searches of neutrinoless double beta-decay |83]. 



from Refs. 87, 88], assume the production of in meson decays and look for visible 
channels in a detector located some distance from the source. The limits at 95% C.L. in 
Refs. []89| , |90t| analyse the data from DELPHI and L3 detectors, looking for from Z^- 
decays. In Fig. we also report the excluded region from neutrinoless double beta-decay 



experiments |91, bounded by dotted (maroon) line, valid if the heavy neutrinos are 



Majorana particles (see further). 
2.2.2 Mixing with 

The bounds on |V^4p come from searches of peaks in the spectrum of muons in pion and 
kaon decays and of the decays of A'4 produced in neutrino beams and e^e~ collisions. 

As already discussed in the case of mixing with Vg, peak searches provide very robust 
and stringent bounds, by looking at pion decays for masses up to 34 MeV, and at kaon 



decays for higher masses. A detailed review is given in Figs. 1 and 2 in Ref. |92] and 



for masses larger than 100 MeV the limits are reported in Fig. The other limits on 
|V^4p are found in decay searches and are also shown in Fig. ^. They come from beam 
dump experiments |87, |8^, 93, 94, 95] and from direct A'4 production in the detectors 
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Figure 3: Bounds on |Ve4p versus 7714 in the mass range 10 MeV-100 GeV. The areas with solid 
(black) contour labeled tt eiy and double dash dotted (purple) contour labeled K ^ eiy are 
excluded by peak searches jS^L Limits at 90% C.L. from beam-dump experiments are taken 
from Ref. |86| (PS191), Ref.^ (NA3) and Ref. @ (CHARM). The hmits from contours labeled 
DELPHI and L3 are at 95% C.L. and are taken from Refs. |^ and respectively. The excluded 
region with dotted (maroon) contour is derived from a reanalysis of neutrinoless double beta decay 
experimental data p^ . 



DELPHI II , L3 [|g] and CHARM ||]. 
2.2.3 Mixing with Ur 

Heavy neutrinos mixed with r neutrinos can be produced either via CC interactions if a r 
is produced or in NC interactions. The only limits come from searches of A^4 decays and 
are reported in Fig. |. The bounds at 90% C.L. from CHARM Q and NOMAD 



assume production via D and r decays. The DELPHI bound at 95% C.L. ||89|] assumes 
A^4 production in decays and with respect to the bound on |V'e4p and |VJ^4p there is r- 
production kinematical suppression for low masses which weakens the constraint for masses 
in the range ~ 2-3 GeV. 

2.2.4 Electroweak Precision Tests 

The presence of heavy neutral fermions affects processes below their mass threshold due 
to their mixing with standard neutrinos [^] and significant bounds can be set by precision 
electroweak data. The effective //-decay constant G^, measured in muon decays, is modified 
with respect to the SM value and can be related to the fundamental coupling Gp as: 

= GF^il-\VeA\''){l-\V^A\^) . (2.10) 
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(GeV) 



Figure 4: Limits on |Vp4p versus 1124 in the mass range 100 MeV-100 GeV come from peak 
searches and from A'4 decays. The area with solid (black) contour labeled K ^iu p% is excluded 
by peak searches. The bounds indicated by contours labeled by PS191 lH], NA3 fs^BEBC |93|, 
FMMF 0, NuTeV H and CHARMII |6| are at 90% C.L., while DELPHI g and L3 @ are 
at 95% C.L. and are deduced from searches of visible products in decays. For the beam dump 
experiments, NA3, PS191, BEBC, FMMF and NuTeV we give an estimate of the upper hmit for 
the excluded values of the mixing angle. 



The fjL — e universality test, done by comparing the decay rate of pions into ev and ^v, can 
be used to constrain the ratio 

for m4 > rriTT |7^. The analysis of experimental data leads to = 1.0012 ±0.0016 

|71], which implies |K=4p < 0.004 at 2a for the least conservative case of |VJj4p = 0. For 



m4 > m,-, the ^ — t universality sets limits on: 

l-\VrA? 



1-IK.-|2^ 



(2.12) 



and can be tested by looking at the r leptonic and hadronic decays which give |V^4p — 
\V^A? = 0.0057 ± 0.0065 ]7l| and |K4p - |V"e4p = 0.0054 ± 0.0064 Q. The most con- 
straining bound on |Ft4P is obtained for |Ve4p, |V^4p = and reads |Vr4p < 0.018 at 2a. 
The unitarity constraint on the first row of the CKM matrix [^] reads 

E IK^™I'= 1 IT/ 12 =0-9992 ±0.0011, (2.13) 

i=i,2,3 ^ 
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Figure 5: Bounds on |yr4p versus 7714 from searches of decays of heavy neutrinos, given in Ref. ||97 
(CHARM) and in Ref. [|| (NOMAD) at 90% C.L., and in Ref. ^ (DELPHI) at 95% C.L. 



and translates into a very strong bound on |V^4p, |VJ^4p < 0.0003 (0.0014), at 1 (2)ct, 
which holds for sterile neutrinos heavier than the A baryon. 

In the presence of heavy singlet neutrinos heavier than half the mass, the invisible 
decay rate of would be reduced with respect to the SM one, r|^j^^, as: 



^-(l-^l^e4p-^I^Ml^-^|K4n. (2.14) 



Z— +inv 



By a standard model fit to LEP data, the effective number of neutrinos is now determined 
to be = 2.984 ± 0.008 pi and provid es a bound on |Vg4p similar to but somewhat 



weaker than the ones obtained by lepton-universality. 

A combined analysis of an old set of unitarity bounds which does not include the 
one from the CKM matrix determination, leads to the following limits at 90% C.L. |V^4p < 
0.012, |y^4|2 < 0.0096 and |K4p < 0.016. If the CKM matrix constraint is included and 
partial cancellations between the contributions of different flavors are taken into account, 



a previous combined study [70| then gives the more robust limits at 90% C.L., |VJ=;4| < 



0.0066, |V/,4p < 0.0060 and |y^4p < 0.018. A very recent analysis @ has updated these 
results using the latest electroweak precision data, except for the CKM observables. They 
find at 90% C.L. 

\VeA^ < 0.003, \V^a\^ < 0.003, \Vr4\^ < 0.006 . (2.15) 

If the constraints from CKM observables are included, we expect the bounds to become 
somewhat stronger, given by |Ve4p < 0.002, |y^4|2 < 4 x IO-MK4P < 0.006 |lOO|. In 
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Table 1: Most stringent model- independent constraints on the mixing elements of the heavy 
neutrino from precision electro-weak measurements. The bounds on |V^4p, t = e,^, r and on 
\VeiV^,4:\ at 90% C.L.. See text for details. 



Mixing element 


Range of m4 


EW Measurement 


\yeAV 


m4 > ©(m^) 


< 0.003 [^] 




m4 > 0{m!y) 


< 0.003 g 




m4 > 0{mr) 


< 0.006 [0] 




10 GeV (100 GeV) [1000 GeV] 


< 0.015 (3.5 X 10-"^) [1.2 X 10'^] 



the following, we take the bound |V;,4|^ < 0.0060 as a conservative reference limit on the 
mixing for comparison with the results of our study. 

Indirect limits on the parameters characterizing heavy sterile neutrinos can be obtained 
from searches for flavour changing neutral current processes such as /i — > 67, ee^e~ 
and fi — e conversion in nuclei [75|. The branching fraction for /i — > e7 induced by the 
mixing with heavy singlet neutrinos is given by [^, 74, 75 1: 



Br(/i — > 67) 



3a 
8^ 



m 



T/* 
em' *^/nn' 



m 



W 



where m' indicates the heavy sterile neutrinos with mass w-at^ 
W boson. The function g{x) is given by 



(2.16) 



and mw is the mass of the 



a{x) 



x{l -Qx + 3x2 _^ - Qx^ \n{x)) 
2(1 -x)4 ■ 



(2.17) 



where g{x) goes from to 1 as x varies from to infinity, 
ing fraction is constrained to be Br(^ — > 67) < 1.2 x 10^^^ 



At present 



101] at 



the branch- 
C.L. imply- 
[1.2 X 10"^] for 



ing that, for one extra sterile neutrino, 1 144^^4! < 0.015 (3.5 x 10' 
777-4 = 10 GeV (100 GeV) [1000 GeV]. Similar constraints are imposed by searches for the 
processes fx — e conversion in nuclei and /i ee^e~ [^. The current strongest bound 
comes from the search for fi — e conversion in Ti for which the branching ratio with respect 
to the total nuclear muon capture rate is constrained to be Br(/iTi — > eTi) < 4.3 x 10^^^ 
at 90% C.L. [102]. For one sterile neutrino, this translates into a bound on the follow- 



ing quartic combination of mixing angles |K;4^^4X]£ l^4l | < 1-3 x 10~ (lOOGeV/7774 
which is weaker than the bounds from — > 67 searches but becomes important at very 
high values of the masses, 7*7,4 ^ 10 TeV. In the presence of more than one sterile neutrino, 
partial cancellations between their contributions are possible, with a consequent weakening 
of the bounds. Future more sensitive searches will further improve these limits. 

The most stringent EW precision constraints are compiled in Table || and include the 
bounds reported above from universality tests and lepton flavor changing processes. These 
bounds are obtained barring cancellations between mixing angles and therefore could be 
weakened if some parameters are of the same order. 

2.2.5 Neutrinoless Double Beta Decay (Oi//?/3) 

The most well studied among AL = 2 processes is neutrinoless double beta decay (Oi^/3/3) 
and the constraints from it deserve special attention. The constraints on ]V^4]2 for a wide 
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range of heavy neutrino masses (10 MeV < 1114 < 100 GeV) are shown in Figs. ^ and |3|. 
For heavy neutrinos with mass, it^n^, ^ 1 GeV, the bound is pl| , 34 1 

y < 5 X 10-^ TeV-i. (2.18) 

m' 

The constraint above is very strong and makes it impossible to observe at colhders the 
like-sign dilepton signature with electrons (see Sec. |^. 



3. Lepton-Number Violating Decays 



The key point for the search of lepton-number violating processes in this paper is to consider 
the substantial enhancement via resonant neutrino production. One thus needs to evaluate 
the decay widths of A'4 to various channels. We consider the decay width of the heavy 
Majorana neutrino in two regimes: when the mass is much smaller than that of the W 
boson and when the mass is larger than the mass of the W boson. Based on this, we then 
compute the AL = 2 decay branching fractions for r lepton and K,D,Ds and B mesons. 

3.1 Decay Modes of Heavy Majorana Neutrino 

3.1.1 Decay Modes of Heavy Majorana Neutrino with mass <C mw 

For the LV low energy tau decays and rare meson decays the resonant contribution is from 
a heavy Majorana neutrino with mass of order MeV to GeV. In this section we discuss the 
decay modes of a Majorana neutrino which is lighter than the W boson, so that <C mw- 
The heavy neutrino decays via charged and neutral current interactions to the modes listed 
below. The partial decay widths of the heavy Majorana neutrino with the leading terms 
in mixing and in the massless limit of the final state particles are given below. The full 
detailed expressions for the same are given in Appendix ^ 





= T{N4 - 




(3.1) 


Y^iP 


= r(iV4 - 




(3.2) 


yIV 


= r(iV4 - 




(3.3) 


Y'^eV 


= r(iV4 - 




(3.4) 


Y^i^2'^e2 


= r(iV4 - 




(3.5) 


Y'^llki2 


= r(A^4 - 


i^ij2^2) = l^iiil'^ ml [ai +5^i^2«2], 


(3.6) 






(3.7) 



where and are charged (neutral) pseudoscalar and vector mesons, /m are the 

meson decay constants and Vqq/ are the CKM matrix elements. 
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All the decay modes listed above contribute to the total decay width of the heavy 
Major ana neutrino which is given by: 



e,p ly e,p e,v 

+ ^ 2r^i^2i^^2 _^ ^ Y"'!^^^^ + ^ r'^'^i'''', (3.8) 

ilMil^k) il/2 il 



where ^,£1,^2 = s^fJ-jT- For a Majorana neutrino, the AL = process — > i~P'^ as 
well as its charge conjugate AL = 2 process l^P^ are possible and have the same 

width r^^. Hence the factor of 2 associated with the decay width of this mode in Eq. (p.^). 
Similarly, the AL = and its charge conjugate AL = 2 process are possible for the decay 
modes l^V^ and — > ^^"fg^ffj ^^^d hence have a factor of 2 associated with their 

width in Eq. Q. 

For the low energy LV tau decays and rare meson decays we consider, the mass of 
the heavy neutrino is in the range 140 MeV < 7714 < 5278 MeV. For this mass range we 
list all the possible decay channels for in Table |^ in Appendix |C|. The mass and decay 
constants of pseudoscalar and vector mesons used in the calculation of partial widths given 
in Eqs. (3J) — (3/7) are listed in Table |^ in Appendix^ 



Next, we get a rough numerical estimate of the total width of the heavy neutrino with 
the decay modes discussed in Eqs. (|3.lD — ( ^?7|) . To do this, we consider the massless limit of 
the decay products of the heavy neutrino, include only leading terms in mixing of 0(| V£4p) 
and ignore small factors like vr and |yC'^^^|2 {-^ calculating the partial decay widths. We 
can only get a rough estimate of the width in this approximation, but it is sufficient to 
see that it warrants the use of narrow width approximation. The two body decays of the 
heavy neutrino have a general form 

p26od, ^ ^^/m^V mP ~ ^i^|144p ~ (10-^3 |y,4p) GeV, (3.9) 

where typical values of 771,4 ~ 1 GeV, fu ^ 0.1 GeV and Gp ^ 10^^ GeV~^ have been 
used. The three body decays of the heavy neutrino have a general form 

where typical values of ?n,4 ~ 1 GeV and Gp ~ 10~^ GeV~^ have been used. The total 
width of the heavy neutrino is then given by 

r^, = (number of decay modes) x {V'^^^'^y + T^body^ 

~ 50 X (10-^^ GeV + 10-^3 GeV) |V^4p ~ (10^" GeV. (3.11) 

As shown above, the width of the heavy neutrino ~ 0(10"^^ |V£4p) GeV is much smaller 
than the mass of the heavy neutrino ~ 0(1 GeV) and we can use the narrow width 
approximation to an excellent approximation. 

Now we look at the lifetime of the heavy Majorana neutrino to determine the decay 
length. The lifetime is given by 

1 1 
""^^ " ~ 10-11 |y,4P GeV ' 

~ 10^^ iF^r^ GeV"^ ~ 6.58 x 10"^^ |V^4r^ s, (3.12) 
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which gives a typical decay length ctki^ ~ 1 x 10~^ |V£4|~^ ni. Note that for a very small 
mixing, |V£4p < 0(10"^), the may escape from the detector if it is not much heavier 
than a GeV. We will take this effect into account in the following studies. 

3.1.2 Decay Modes of Heavy Majorana Neutrino with mass 7714 > my/ 

In this section we discuss the decay modes of the Majorana neutrino which is heavier than 
the W gauge boson, so that 7714 > mw. The decay modes of the heavy Majorana neutrino 
are to a or a Z gauge boson plus the corresponding SM lepton. The partial decay widths 
for longitudinal and transverse gauge bosons , in static heavy neutrino frame are 





= r(iV4 - 




Y^Wt 


= r(Ar4 - 


rw+) -- 




= V{Ni - 


viiZl) = 


YUiZt 


= r(iV4 - 





r(iV4 ^ i^Wr 



\Vu\^ mi (l-mf, (3.13) 



= r(iV4 - 



T . 



327r 



\Ve4\ rn4 (1 - nwf 



32-K cos^ 



\Vu\ (1 - ^izf 



(3.14) 
(3.15) 

(3.16) 



where are the masses of the gauge bosons scaled by the mass of the heavy neutrino and 
are given by /ij = m1/m\. To obtain the total decay width for A'4, we sum over the charged 
leptons £ and as discussed earlier include the AL = process A^4 ^~^lt ^^^^ ^ 
charge conjugate AL = 2 process A^4 — > i~^W^j,. Hence the factor of 2 associated with the 
decay width of these modes in the expression for the total width below. 



(3.17) 



In Eqs. ( 3.15| ) — ( 3.17 ), we have used the relation (see Appendix for details) 



El''. 



m4 



m=l 



ll=e 



Since 



UU^ + VV^ = I. 



(3.18) 



Ignoring terms of order |V^4|^ we have 



EK'l\"-^E\^u\'- (3-19) 



In this approximation, the total width of a heavy Majorana neutrino can be written as 

« 2^ \Vi4\ for 7774 > mw, 

riV4 < ^ (3.20) 



oc E I V£4|^ Gpml{f'lj + 7774) for 7774 <C mw, 



where the expression when 7774 <C myy is obtained from Eq. (3.8) and /m are the meson de- 
cay constants. We note that the approximate form of the total width as given in Eq. ( 3.2C| ) 
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Figure 6: (a) Top: decay width and (b) bottom: decay length (normalized by J2e l^4l^) versus 
mass of heavy Majorana neutrino for real and virtual weak bosons with the inclusion of Higgs decay 
channel for uih = 120 GeV . 

is only for intuitive purposes to infer the general behaviour of the total width as a function 
of mass. The precise expressions for the total width of the heavy Majorana neutrino as 
given in Eqs. ( |3.8[ ), ( 3.171 ) and ( p.9| ) have been used in the numerical analysis. 



It should be noted that in the SM, if A^4 is heavier than the Higgs boson, then the 
decay to a Higgs will be present and the partial width is given by 

^ r(7V4 ^ ueH) = |y,4|' ml (1 - iji^f. (3.21) 

In Fig. we plot the decay width of the heavy Majorana neutrino versus its mass normalized 
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Figure 7: (a) Left: branching fractions for decay of heavy Majorana neutrino into W* and Z* 
bosons with varying heavy neutrino mass; (b) right: branching fractions for decay of heavy Majorana 
neutrino into longitudinal and transverse gauge bosons in static heavy neutrino frame with the 
inclusion of Higgs decay channel for mn = 120 GeV . 

by the common factor |V£4|^. We can see in Fig. |6|(a) that for a heavy neutrino with 



mass > m\y, the decay width increases as Gpryi'l as given in Eq. (|330| ). Given the 



rather small mixing parameter, the width remains narrow even for 1114 0(1 TeV). For a 
lighter neutrino with 772,4 n^-W: the width can be very small. The proper decay length is 
presented in Fig. ^(b). We see from this that for 7714 ^ 20 GeV and |F^4p ^ 10~^ from 
Fig. ^, we have cr ~ 1 fim. 

In Fig. 0(a) we plot the branching fractions of the heavy Majorana neutrino decay to 
Wi and Zi' versus varying heavy neutrino mass m^. In Fig. ^(b) we plot the branching 
fractions for the decays into longitudinal and transverse gauge bosons in static heavy 
neutrino frame. When the neutrino mass is large, it mainly decays to longitudinal gauge 
bosons and Br(iV4 ~ Br(7V4 Zv) = Br(iV4 Hv) = 25%. In terms of the 

search at hadron colliders, we prefer to adopt the Wi mode since we wish to reconstruct 
the full event including the lepton number. 

3.2 Lepton-Number Violating Tau Decays 

In this section we examine tau decays into an anti-lepton and two mesons 

r-{pi) ^ ^+(P2) Aff (gi) M^{q2) (3.22) 
which is a process with AL = —2. The decay amplitude for the above process is given by 



iM = 2GlV^^^V^^^fMjM,V:,V,\m4 i ^z-^^. +(91^-92), 

(3.23) 

where V^^^^ and ai'e the quark flavor mixing element and the decay constant for 
the meson Mi respectively. From this decay amplitude, we can calculate the transition 
rate F^y and the branching fraction normalised by the tau decay width. In Appendix 
we give the calculations and the full expressions for the decay branching fraction of the 
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Table 2: Mass and mixing elements of heavy neutrino and the decay mode constraining them with 
the corresponding experimental bounds on branching fractions. Bounds for AL = 2 tau decays are 
from Ref . [|o| 



Mixing element 


Range of (MeV) 


Decay mode 


^exp 




140 - 1637 
140 - 1637 
494 - 1283 


T~ e"*'7r~7r^ 
T- e+R-R- 


2.7 X 10"'' 

1.8 X lO"'^ 
1.5 X lO"'^ 




245 - 1637 
245 - 1637 
599 - 1283 


/i"'"7r~7r~ 
T- fi+R-R- 


0.7 X lO-'' 
2.2 X lO"'^ 
4.8 X 10-"^ 



process ( 3.22 ) in terms of the mass of heavy neutrino, 771,4, and the mixing \VriViii\^. To 
understand the physical picture, we can express the branching fraction in an intuitive form, 
in the massless limit of the final state particles, as 

Br = = r ^ ^ 



3^/-, Ir ^r-a f2 f2 \jrCKMjrCKM\2 \jr t/ 12/1 m| \ / 7714 

-TT[l - -dMiMijf^FJMjphl^Mi ^OV/a I ' '^^l (^1 " ^ j 



~ 10-3 ^yCKMyCKM^2 i^^^^^^i^ (3_24) 

where we have used typical values of rui ~ 1 GeV, /a/. ^ 0.1 GeV, Gp ~ 1 x 10^^ GeV^^ 
and ~ 10^^^ |V£4p GeV. From the simple expression given above one can easily make 
a rough estimate of the required sensitivity and hence the feasibility of observation in terms 
of the mixing parameters for a given model. 

A direct search for LV tau decays has been made at the BaBar detector and the 



limits on the branching fractions were reported in Ref. | 103 |. The experimental limits for 
various decay modes are typically of the order of 10~ , as given in Table |. Fr om the non- 
observation of the LV tau decay modes one can determine bounds on the mixing parameters 
\ViiVri\'^ as a function of the heavy neutrino mass 7714. To do this in a comprehensive 
manner, we carry out a Monte Carlo sampling of the mixing parameters and the mass of 
the heavy neutrino. For simplicity, the mixing elements ^44, V^4 and VV4 are allowed to 
vary in the range from to 1. The ranges of mass sampled for the heavy neutrino are 
listed in Table ^ for the various tau decay modes. We only sample the range of masses that 
lead to a resonant enhancement of the width as the other mass regions have very small 
transition rates as discussed earlier. We then calculate the transition rates and branching 
fractions over the entire range of mixing and mass of the heavy neutrino and the results of 
the Monte Carlo sampling are discussed next. 

The relevant mixing parameters \Ve4Vr4\ and |V^4V^4| are probed as a function of the 
heavy neutrino mass 772-4 and are shown in Fig. ^(a) and Fig. §(b), respectively. Under the 
assumption that the detector was able to reconstruct all the signal events, the region above 
the curves is excluded by the current direct experimental search for LV tau decays. The 
most stringent bound on iVe^Vr^l is of O{10~^) and comes from c'^it^tt^ . The most 

stringent bound on | V^4Vr4| is also of 0(10"^) and comes from r~ fi'^Tr^-K^ . This is three 
orders of magnitude more sensitive than the limits from precision electroweak data which 
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I 1 2 "-i 

constrain the square of the mixing \Vi/^\ to be less than few times 10~ . In the absence 
of detection of LV processes the constraints on mixing from peak searches, accelerator 
experiments, reactor experiments and others (collectively called laboratory constraints here 
and henceforth) described in Fig. |2| — Fig. ^ are also applicable here. In the mass region 
probed by LV tau decays the most stringent current constraints are 

|V^4p < 10~^ — 10~® and |V^4p < 10~^ — 10~^. This would roughly translate into constraints 
on \Ve4Vri\ < 10~^ — 10~^ and |V"^4Ft4| < 10~^ — 10~® which are comparable to the limits 
from LV tau decay modes. We explore more combinations of mixing elements and also 
provide better constraints on mixing in some mass regions. To summarize, the constraints 
on mixing from LV tau decays are always competitive with or better than precision EW 
constraints and laboratory constraints in the corresponding mass region. The experimental 
bounds can improve in future and an order of magnitude improvement in the experimental 
branching fraction will give approximately an order of magnitude improvement in the 
constraints for the mixing parameters |Vf4V^4|. More importantly, a detection in one of 
the laboratory experiments implies the existence of a sterile neutrino while a detection in 
one of the modes studied in our analysis would imply LV and hence the existence of a 
Major ana neutrino. 

It should be noted that the intermediate heavy Majorana neutrino is treated as a real 
particle which propagates before decaying. If it exits the experimental apparatus prior to 
decaying, then the signal corresponding to the AL = 2 process cannot be reconstructed 
and no bound could be deduced from the non-observation of such events. In Figs, ^(a) 
and ^(b) , we provide an estimate of the bound on the mixing parameters which takes into 
account the probability of the heavy Majorana neutrino to decay within the detector of 
size Lexp. This probability is given by 

P = l- exp(-Lexpr^) (3.25) 

and for small masses and/or small mixing parameters and consequently long decay lengths, 
it can be approximated with P ~ LexprAf- We take -Lexp = 10 m, the typical size of the 
detectors used in the experiments under consideration. For simplicity, we take A^4 to be 
relativistic but we keep its gamma factor 7 = 1, as a more precise value requires a full 
understanding of the experimental setup. We assume 1^441 = |V^4| = \Vr4\. An estimate of 
the realistic bound on the mixing parameter {Ve^Vril is then given by 

|K4K4|(= IK4I') = ^\Ve4\lo/{L,^p^N0), (3.26) 

where |K;4|^ is the bound obtained assuming that all the A'^4 decay in the detector and 
discussed above, and r^vo is the decay rate for a fully active heavy Majorana neutrino, i.e. 
when the mixing parameter {Vi^l = 1 for i = e, /x,r. The bounds remain unchanged for 
large values of the mixing angle and / or large values of m^, as the decay length in these cases 
is very short. However, the most sensitive limit on 1 1441^4 1(= |V^4p) coming from r — > evrvr 
searches gets weakened to ~ 4 x 10"'' (4 x 10"^) (1 x lO'^) for 7714 = 0.2 (0.5) (1.0) GeV. 
Similarly, the searches for r — > ^vrvr allows to set a bound on | V^4 1/^41 (= |^4p) which 
weakens to \V^4Vr4\ < 1 x lO""^ (2 x 10"^) (1 x lO^^) for m4 = 300 (600) (900) MeV. A 
detailed analysis taking into account the experimental setup should be performed in order 
to obtain more precise bounds. 
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Figure 8: (a) Left: excluded regions above the curves for {VeiVril versus TO4; (b) right: same as 
(a) but for |V^4T/r4|- The thin black lines correspond to the estimate of the bound (for r — > eirir 
and r fimr) once the probability of N4 decay in the detector is taken into account. 



3.3 Lepton-Number Violating Rare Meson Decays 



We now investigate the LV processes in which a meson decays into two like-sign leptons 
and another meson 

M+{qi) ^ e+ipi) e+{p2) M^{q2). (3.27) 

These decays are similar to the tau decay modes described in the previous section. The 
decay amplitude for the above process is given by 



iM 



V 



CI 



+ ipl 



P2) 



+ {PI 



P2) 



(3.28) 



(3.29) 



where iM^ and iM^ are the decay amplitudes when the meson M2 is a pseudoscalar or 
vector meson respectively and V^^^ and /m; are the quark flavor mixing element and the 
decay constant for the meson Mi respectively. From this decay amplitude, we can calculate 
the transition rate I^^^ and the branching fraction normalised by the decay width of the 
meson Mi. In Appendix^, we give the calculations for the decay branching fraction of the 
process ( f3.27 ) in terms of the mass of heavy neutrino, m^, and the mixing \Vi-^^4Vi^4\. We 
can express the branching fraction in an intuitive form, in the massless limit of the final 
state particles, as 



Br 



r 



Ml 
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Table 3: Same as Table |^ but for AL = 2 rare meson decays. The experimental bounds are from 
Ref. @, the bounds for D+ e+e+TT-{K-) are from Ref. |o|. 



Mixing element 


Range of 7714 (MeV) 


Decay 


mode 






140 - 493 


K+ ^ 




6.4 X 10"^^ 




140 - 1868 


D+ ^ 


e + e^vr^ 


3.6 X 10-6 




494 - 1868 


D+ ^ 




4.5 X 10"*^ 




140 - 1967 


D+ 

s 




6.9 X 10""^ 


\Ve4\^ 

\ ' 6'i 1 


494 - 1967 


Dt 

S 


e+e+K- 


6.3 X 10""^ 




140 - 5278 


B+ ^ 




1.6 X 10-'^ 




494 - 5278 


B+ ^ 




1.0 X 10"*^ 




776 - 5278 


B+ ^ 




2.6 X 10"*^ 




892 - 5278 


B+ ^ 


e+e+K*- 


2.8 X 10"*^ 




245 - 388 


K+ 




3.0 X lO"'-* 




245 - 1763 


D+ 




4.8 X 10"^ 




599 - 1763 


D+ ^ 




1.3 X 10"^ 




881 - 1763 


D+ 


o~ 


5.6 X 10""^ 




997 - 1763 


D+ 




8.5 X 10""^ 


1 1 


245 - 1862 


Dt 

S 




2.9 X 10"^ 




599 - 1862 


Dt 

^ S 




1.3 X 10"^ 




997 - 1862 


Dt 


p+p+K*- 


1.4 X 10"^ 




245 - 5173 


B+ ^ 


P^ p^TT^ 


1.4 X 10"*^ 




599 - 5173 


B+ ^ 


p+p+R- 


1.8 X 10"*^ 




881 - 5173 


B+ ^ 


Li^ 0^ 


5.0 X 10^6 




997 - 5173 


B+ ^ 


p+p+K*- 


8.3 X 10^6 




140 - 493 


K+ ^ 


e'^p'^ir^ 


5.5 X 10-^" 




140 - 1868 


D+ ^ 


e'^p~^'ir~ 


5.0 X 10"^ 




494 - 1868 


D+ ^ 


e+p+R- 


1.3 X lO""^ 




140 - 1862 


Dt^ 




7.3 X 10""^ 




494 - 1967 


Dt^ 


e+p+R- 


6.8 X 10-"^ 


140 - 5278 


B+ ^ 


e'^p'^ir^ 


1.3 X 10-6 




494 - 5278 


B+ ^ 


e+p+R- 


2.0 X 10^6 




776 - 5278 


B+ ^ 


e^p^p" 


3.3 X 10-6 




892 - 5278 


S+ ^ 


e+p+R*- 


4.4 X 10-6 



~ (10-16 GeV) TM, |V™vCKM|2 |v,^4V,,4|, ^ (3.30) 

where we have used typical values of rui ~ 1 GeV, /m, ~ 0.1 GeV, Gj? ~ 1 x 10-^ GeV"^, 
rAr4 ~ IQ-i-i |V£4p GeV and tm-^ is in seconds. Using the values for the lifetimes of the 
mesons in Appendix the branching fractions for the various mesons are given by 

Br(i^) ~ in^f^^y^f^^P \V,,,V,,,l (3.31) 
Br(D, B) ~ 10-^ l^^^f^'^^f^^l' \Vi,^V,,,l (3.32) 
Br(i?,) ~ 10-5 l^^^f^^^f^P |F.,4^..4|. (3.33) 
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As mentioned earlier, with the simple expressions above one can easily make a rough 
estimate of the required sensitivity and hence the feasibility of observation in terms of the 
mixing parameters for a given model. 

Searches for rare meson decay modes have been made in numerous experiments. Table 
m summarizes the current experimental limits on branching fractions given by Refs. |]99| , 
|10^ ]. From the non-observation of these LV rare meson decay modes one can determine 
constraints on mixing parameters |Vfj4V^24l ^ ^ function of the heavy neutrino mass m^. 
To do this in a comprehensive manner, we carry out a Monte Carlo sampling of the mixing 
parameters and the mass of the heavy neutrino similar to tau decay. The mixing elements 
V^4 , V^4 and 1/^4 are allowed to range from to 1 for simplicity. Only the range of mass 
that leads to a resonant enhancement of the width is sampled for the heavy neutrino and 
listed in Table |3| for the various meson decay modes. The transition rates and branching 
fractions are then calculated over the entire range of mixing and mass of the heavy neutrino 
and the results of the Monte Carlo sampling are discussed next. 

1 2 I 

For the various decay modes, the mixing parameters probed are \Vei\ , |Ve4V^4| and 
|y^4|^ depending on the final state leptons. Again, we plot the excluded region of the 
mixing parameters as a function of neutrino mass, as shown in Fig. | for \Ve4\ , Fig. 
for \VeAVf,4\ and in Fig. || for 1 1^,4 p. The regions above the curves are excluded by the 
current direct experimental searches for LV meson decays. First we plot the limits which 
can be derived assuming that all A'^4 decay in the detector and give a positive signature. 
The most stringent constraints are from the lfl2'^~ mode with mixings of O(10~^) 

I 1 2 I 1 2 

excluded for 1^441 , \ye4V^i4\ and {V^^l . This is six orders of magnitude more sensitive 
than the limits from precision electroweak data which constrains the square of the mixing 

I 1 2 '\ 

\Ve4\ to be less than few times 10" . Next in sensitivity are the D and Ds decay modes 
with constraints of order few times 10~^ which are similar to the constraints from precision 
electroweak data. The bounds for the same mixing elements are much weaker in the mass 
range above 2 GeV. Even though the limits are weak in this region, it is important not 
to neglect the experimental study of these processes. It only implies that there is a large 
parameter space available for the mass and mixing of heavy neutrinos. 

As discussed for the AL = 2 tau decays, in the absence of detection of LV processes 
the laboratory constraints on mixing described in Fig. ^ — Fig. |5| are also applicable here. 
In the mass region probed by LV meson decays the most stringent laboratory bounds 
are 1^64 P < 10"^ - 10-^ |y^4p < 10"^ - IQ-^ for ^4 < 2 GeV and |T/^4p < 10"^ for 
m4 > 2 GeV. This would roughly translate into constraints on |T44V^4| < 10~^ — 10~^ 
for m4 < 2 GeV and \VeiV^i\ < 10"^ - lO"*" for 7714 > 2 GeV. It should be noted that, 
if these experiments were able to fully reconstruct the signal, the limits from K meson 
decays would be better than the laboratory constraints by at least an order of magnitude 
in the corresponding mass region. In fact, the constraints on \Ve4,\ from the kaon decay 
mode — > e'^e^7T~ would be more stringent than even the constraints from {)v(3j3 shown 
in Fig. ^. Usually OuPP experiments have the best sensitivity as they have an advantage 
of a large "effective luminosity" resulting from the large number of nuclei available for 
decay. The meson (and tau) experiments on the other hand have a small luminosity 
coming from a limited number of mesons (taus) produced in accelerators compared to 
the number of nuclei in Oi^PP experiments. It is interesting to note that the resonant 
enhancement in the case of the K meson decay is able to match or improve over the large 
"effective luminosity" of Oi^/3/3 experiments. In conclusion, the constraints on mixing from 
LV meson decays are competitive with the precision EW constraints and all the laboratory 
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Figure 9: Excluded regions above the curves for \Ve4\ versus 777.4 from e+e+Mj" searches. 

The thin black line corresponds to an estimate of the bound from K'^ — > e^e^Tr^ once the proba- 
bility of decay of N4 in the detector is taken into account. 



constraints, potentially even Oi^/3/3, in some mass regions. But again, we emphasize that 
the aim of our analysis is to study LV processes and hence Majorana neutrinos. 

We have also taken into account the fact that, for small mixing, only part of the heavy 
sterile neutrinos produced will decay in the detector. We have considered L^xp = 10 m, 
IK4I = |^/i4| = \Vt4\ and the gamma factor of A^4, 7 = 1, for simplicity. In this case, 
as discussed for the AL = 2 tau decays, the bounds get sensibly weakened. An estimate 



of these bounds is reported in Figs, y, |10| and |11| by thin black lines. We see that the 
bounds get significantly weakened by few orders of magnitude for K eevr, K e/xvr and 
K — > /i^uTT and a careful analysis of these searches should be performed to find the detailed 
bounds on the mixing angles. 

The sensitivity of current direct experimental searches are not adequate to constrain 
mixings for some decay modes. The theoretically allowed branching fraction versus mass 
7714 for such modes is given in Fig. 12. As we can deduce from Table |3| and Fig. 12 all the 



modes are very close to start being probed by direct experimental searches. The experimen- 
tal bounds on branching fractions can improve in future and similar to tau decay modes, an 
order of magnitude improvement in the experimental branching fraction will give approx- 
imately an order of magnitude improvement in the constraints for the mixing parameters 

I 1 2 

IViiiVi^il. Currently we do not have any constraints on the mixing parameter |V^4| from 
LV rare meson decay modes. Only very weak constraints for BR(i? Xt^t~) < 0(5%) 
exist in a theoretical analysis |105| ]. The similar signature t^t^M^ is a possible 

decay mode that would bound {Vril'^ and should be pursued. 
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Figure 10: Same as Fig. |9| but for \VeiV^4\ from 
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Figure 11: Same as Fig. |9|but for |V;i4|^ from Mi /i^/i+Afj searches. 



4. Collider Signatures 

In this section we study heavy Majorana neutrinos at hadron colHders. The most distinctive 
channels of the signal involve like-sign di-leptons. It was first proposed in Ref. in the 
context of the left-right symmetric model, and subsequently studied in Ref. 



3C, M, 32, 33| 




m4 (GeV) 

Figure 12: Branching fraction versus heavy neutrino mass 1714 for decay modes tlt^^2 
not yet constrained by direct experimental searches. The regions below the curve are theoretically 
allowed. 




Figure 13: (a) Left: Feynman diagram for like-sign dilepton signature via WW fusion in hadronic 
collisions; (b) right: the exchanged coherent diagram which is same as heavy neutrino production 
and decay. 

We discuss the signatures for a heavy Majorana neutrino and the sensitivity to probe the 
parameters 772.4 and V^/^ at the Tevatron and the LHC. 

As for the production of a heavy Majorana neutrino at hadron colliders, the represen- 
tative diagrams at the parton level are depicted in Fig. |l^, with the exchange of final state 
leptons implied. The first diagram is via WW fusion with a f-channel heavy neutrino A'4 
exchange, directly analogous to the process of Qvf3l3. The second diagram is via s-channel 
A^4 production and subsequent decay. Although in our full calculations, we have coher- 
ently counted for all the contributing diagrams of like-sign dilepton production including 
possible identical particle crossing, it is informative to separately discuss these two classes 
of diagrams due to their characteristically different kinematics. 

The scattering amplitude for the process in Fig. |l^(a) is proportional to Vi-^^V^^^ and 
the cross section can be expressed as 

a {pp ^ W^W^ ^ eflfx) = (2 - 6e,e,) \V,,iV,,^\^ ao{WW), (4.1) 

where <to(1^W^) is the "bare cross section", independent of the mixing parameters. We 
show the bare cross section at the LHC energy of 14 TeV versus the heavy neutrino mass 
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Figure 14: The bare cross section cro(W^T^) versus mass of the heavy neutrino 7714. 



in Fig. 14. This cross section can be at the order of tens of femtobarns. However, due 
to the large suppression of the small flavor mixing to the fourth power, the cross section 
is rather small. This process was calculated in Ref. |20| ] under the effective vector boson 
approximation. The authors of Ref. [^0| obtained significantly more optimistic results than 
ours. Further scrutiny indicated that they missed a factor of and their result 

should be scaled down by this factor. The corresponding curve for the Tevatron is not 



shown in Fig. 14 as the bare cross section is smaller by nearly two orders of magnitude. 
Including the small mixing element (to the fourth power) further reduces the cross section 
drastically with no hope of detection at the Tevatron via this mode. 

By far, the dominant production process of heavy Majorana neutrino in hadronic 
collisions is the diagram shown in Fig. |l3|(b). We calculate the exact process, but it turns 
out to be an excellent approximation to parameterize the cross section as 



aipp ^ if If W^) ^{2- 5,,,,) aipp ^ ifN,)BriN, ^ ifw^) oc i^^^'^y'' . (4.2) 



This observation allows us to study the process in a model-independent way. We can 
rewrite the cross section in a factorized form 

(^{PP ^ if if ^ If £f j f) = (2 - %,,) S,,e, (4.3) 

where cro(-^4)) called the "bare cross section", is only dependent on the mass of heavy 
neutrino and is independent of all the mixing parameters when the heavy neutrino decay 
width is narrow. As seen in Fig. ^, this is indeed the case for 771,4 ^ 1 TeV once we fold in 
the constraints | Vg4p < 0(10"^) from precision EW measurements. We calculate the exact 
cross section for the dilepton production and use the definition Eq. (|4.3|) to find the bare 



cross sections aQ{N4), which are shown in Fig. 15 at the Tevatron and the LHC energies 
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Figure 15: The bare cross section (To(iV4) versus mass of heavy Majorana neutrino TO4 for the 
Tevatron {pp at 1.96 TeV, sohd curve) and the LHC {pp at 10 and 14 TeV, dotted and dashed 
curves, respectively). 



versus the mass of the heavy Majorana neutrino. Due to the fact that the LHC will start 
its operation at 10 TeV, we have calculated the cross sections at both 10 and 14 TeV cm. 
energy. The production rate is increased at the higher energy by a factor of 1.5, 2.0, 2.5 
for 1714 = 100, 550 and 1000 GeV, respectively. We will mainly present our results at 14 
TeV for the rest of the paper. An obvious feature of the cross sections is the transition 
near the W mass. For < Mw the cross section is nearly a constant due to an on-shell 
W production via the Drell-Yan mechanism with its subsequent leptonic decay to £^N4. 
For m4 > Mw the cross section falls off sharply versus and the on-shell decay goes like 
N4 -^£^W^^ £± jii2. 

The flavor information of the final state leptons is parameterized by 

In general the two final state charged leptons can be of any flavor combination, namely, 

e^e^, e^n^, e^T^, n^n^, ^i^t^ and t^t^. (4.5) 

The constraint from Oi^PP as given in Eq. ( |2.18| ) is very strong and makes it difficult to 
observe like-sign di-electrons e^e^. The events with r leptons will be challenging to recon- 
struct experimentally. We will thus concentrate on clean dilepton channels of fJ-^fJ-^ and 
H^e^, although we will comment on our proposal to include the r modes. The correspond- 
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ing mixing parameters in our notation will be 



|4 

'-^MM ~ v^r iTA i2' '-^e/i — 2' 



l^4| _ |K4V'^4| 



respectively. Given the smallness of |V^4p, we can further simplify our study by exploring 
only two cases: an optimistic case |V^4p ^ |^T4p, \Vei\'^ and a generic case |V^4p ~ 
IK4P > IK4P) which lead to 

y^ii\'^ (optimistic) _ J |Ve4p (optimistic) 



- [ ^\V^,\^ (generic) ' '^'^^ " | i|Fe4p (generic) ' ^^'^^ 
4.1 Search for Like-sign Dilepton Signals at the Tevatron 

We now consider the search for A'4 at the Fermilab Tevatron, which is currently running 
at a cm. energy of 1.96 TeV in pp collisions. We concentrate on the clean like-sign fJ-^^J,^ 
mode 

pp n^fi^ jij2 X, (4.8) 

where X is some inclusive hadronic activities common in hadronic collisions. To quantify 
the signal observability, we first impose the basic acceptance cuts on leptons and jets to 
simulate the CDF/DO detector coverage 

p^>5GeV, |7?''|<2.0, p^>10GeV, |rr''| < 3.0. (4.9) 

We also smear the lepton momentum by a tracking resolution and the jet energy by hadronic 
calorimeter resolution as 



1.5xl0->^, ^ = ^=0 3%, (4.10) 

where pj, and Ej are in units of GeV. 

The signal events we are searching for have very unique kinematical features. For the 
purpose of illustration, we choose m4 = 20, 50 GeV (below mw threshold) and 100 GeV 
(above mw). First of all, there are two well-isolated like-sign charged leptons. This is shown 
in Fig. |^(a) by a normalized distribution of the minimal isolation Ai?£j = At]^ + A^^. 
Second, there is essentially no missing transverse energy. However, realistically, the de- 
tectors have finite resolutions as simulated by the Gaussian smearing given in Eq. ( 4.1C| ). 



Consequently, there is always some misbalance in the energy-momentum measurements, 
which is attributed to the missing transverse energies and is plotted in Fig. ^(b). Thirdly, 
due to the existence of an on-shell in the signal process, one would expect to recon- 
struct it by an invariant mass either from the 2£2j system m{Ujj) (in the case of DY 
production) or from the di-jets mpj) (in the case of decay). This is demonstrated in 
Figs. |^(c) and (d), respectively. The above kinematical features motivate us to impose 
the following event selection cuts 

Ai?^™ > 0.5, (4.11) 
60 GeV < either m{iijj) or m{jj) < 100 GeV, (4.12) 
< 20 GeV. (4.13) 
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Figure 16: Normalized distributions a~^da/dX for TO4 = 20, 50 and 100 GeV at the Tevatron for 
(a) upper left: the minimal isolation Ai?^*"; (b) upper right: the missing transverse momentum 
^r^] (c) bottom left: the 2i2j system invariant mass m{i£jj); (d) bottom right: the di-jet invariant 
mass m{jj). 



These cuts are highly efficient in selecting the signal events. We illustrate this in Table 
^, in which we calculate the signal rates with the consecutive cuts for = 60 GeV and 
|V^4|^ = = 5 X 10"3 > 1 144 p. Note that the choice of mixing elements used in the 

illustration is motivated by constraints from precision EW measurements. However this is 
for illustration purposes only and in our full analysis we have kept S^^ as a free parameter. 

At the Tevatron energies, the SM contribution to the like-sign dilepton events is rather 
small. The leading background of this type comes from the top-quark production and its 
cascade decay via the chain 

t W+b £+ui b, (4.14) 
t W~b -^W- cue i+. (4.15) 

The background rates and survival probabilities with the consecutive cuts are also given 
in Table ^ We see that the tt background is essentially eliminated by the selective 
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Table 4: The representative signal and background cross sections at the Tcvatron, for ^J, ^J. jj and 
the efficiencies with the consecutive cuts. For illustration, we have used TO4 — 60 GeV, |VJj4|^ = 

|K4|2 = 5 X 10-3 > IK4I'. 





No cut 


Basic cut (4.9) 


+AR (4.11) 


+m{jj),m{Ujj) (H.12D 


(4.13) 


Signal 
a{Po) 
eff. 


319 


108 
33% 


99 
92% 


96 
97% 


96 
100% 


tt Bkg 
a(fb) 
eff. 


78.4 


58.2 

74% 


1.85 
3.2% 


0.04 
2.2% 


0.005 
12.5% 




Figure 17: (To(A^4) with varying heavy neutrino mass 7714 after all the cuts. The two cases corre- 
spond to muon rapidity acceptance at DO and CDF. 



cuts. We have also considered other SM backgrounds coming from the production of 
W^W^jj, W^Zjj. After the selective cuts, all these backgrounds are negligibly small. 



In Fig. 17, we plot the bare cross section aQ{N4) with the basic cuts of Eq. (|4.9|) as well 



as the selection cuts Eqs. ( 4.11| ) — ( 4.13| ). The reduction in rate is mainly due to the basic 



acceptance cuts. For comparison, we have also included two choices of pseudo-rapidity 
cut |r?(/i)| < 2 and \r]{n)\ < 1.5. We now consider the statistical significance of the signal 
observation. In the absence of background events, we use Poisson statistics to determine 
the search sensitivity. We take a signal with 95% Confidence Level (as this is very close 
to 2a we call it a 2a effect henceforth) to be 3 events. We can thus translate this to the 
sensitivity to the mixing parameter 

(2 - 6i,i,)aoiN4) Se.e^ L>3, (4.16) 

where L is the integrated luminosity. 

The CDF collaboration at the Tevatron has successfully studied the events with like- 



sign dileptons in a different context [lOt]. Given our event selection, in Fig. |l^ we estimate 
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the sensitivity reach for the mixing parameters versus at the 2a (solid curves) and 5a 



(dashed curves) level at the Tevatron. In Figs. 18(a— b) (upper-left and upper-right), the 
sensitivity is shown for S^^ with 2 and 8 fb^^ integrated luminosity. The horizontal dotted 
lines are the constraint on S^^ ~ l^4p < 6 x 10"'^ from an analysis of precision EW 
measurements @. The DELPHI [||] and L3 @ bounds are also given for comparison. 
We find that the Tevatron has the potential to reach the following sensitivity for the mass 
of the heavy neutrino 

J 40 - 130 GeV for 2cr with 2 fb"^; 

\ 10 - 180 GeV (50 - 120 GeV) for 2a {5a) with 8 fb^^ ^ ^ 

Alternatively, the sensitivity for the mixing parameter can be 

f 2 X 10-5 for 2a with 2 fb"^; . . 

\ 5 X 10-6 (2 X 10-5) for 2a {5a) with 8 fb'^ ^ ^ 



Similar to Figs. 18(a— b), Figs. 18(c— d) (lower-left and lower-right) show the results 
for Se^ instead. The lower dotted curve in Fig. |l8|(d) is the bound on ~ from 
Qvf3[3. We have assumed the same detection efficiencies for n and e. With this assumption, 
the slightly better reach for S'g^ compared to S^^ is due the factor of two difference in 
total rate with identical and nonidentical particles as evident from Eq. ( [4.31 ). With 2 fb~^ 
luminosity, the sensitivity to \Vei\'^ is not close to the stringent bound from the Oi^PP 
decay as seen in Fig. ^. We see from Fig. |l8|(d) that with 8 fb-^ luminosity, the Tevatron 
sensitivity for Sefi may reach the level of the current bound from Oz^/3/3. From Eq. ( [4.16| ), it 
is straightforward to obtain future sensitivity to mixing parameters (S^^^, Sefj.) by a simple 
scaling of the luminosity. 

4.2 Search for Like-sign Dilepton Signals at the LHC 

At the LHC with a cm. energy of 14 TeV in pp collisions, we adopt the basic acceptance 
cuts on leptons and jets as 

p^>10GeV, |t/^| < 2.5, |4 > 15GeV, \r]^<2.5. (4.19) 

The efficiency of these cuts increases with heavy neutrino mass and is 50% for 7714 = 
200 GeV and 80% for 7714 = 800 GeV. The smearing parameters to simulate the AT- 
LAS/CMS detectors are flO^ 

^ = 36 X 10-5 p^, ^ = _^ © 5%, (4.20) 



Pt Ej 



where pj, and Ej are in units of GeV. 

We again present the characteristic kinematical distributions for the signal. Fig. p^(a) 
shows the normalized distribution of the minimal isolation ARij. The simulated missing 
transverse momentum after the energy- momentum smearing is plotted in Fig. ^(b). The 
invariant masses of the 2i2j system m{iijj) and the di-jets m(^jj) are demonstrated in 



Figs. 19(c) and (d), respectively. We thus design the selection cuts at the LHC as 

Ai2^™ > 0.5, (4.21) 

60 GeV < either 77i(^^ii) or m{jj) < 100 GeV, (4.22) 

< 25 GeV. (4.23) 
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Figure 18: The Tevatron sensitivity to the mixing parameters versus 7714 (a) upper-left: 2(t and Scr 
sensitivity of S^^ with 2 fb^^ integrated luminosity; (b) upper-right: same as (a) but with 8 fb^"'^ 
integrated luminosity; (c) lower-left: 2cr and Scr sensitivity of 5*6;^ with 2 fb^^ integrated luminosity; 
(d) lower-right: same as (c) but with 8 fb^^ integrated luminosity. The horizontal dotted lines in 
(a) and (b) are the constraint on ~ l^4p < 6 x 10^'^ from an analysis of precision EW 
measurements [Q. The DELPHI and L3 bounds are also given here for comparison. The 
lower dotted curve in (d) is the bound on Se^ — 1^4^ from 



These cuts are highly efficient in selecting the signal events. We illustrate this in Table |5|, 
in which we calculate the signal rates with the consecutive cuts for m4 = 200 GeV and 

I 1 2 I 1 2 I 1 2 

\Vf,i\ = \yri\ = 5 X lO"'^ > \Ve4\ . A gain the choice of mixing elements is motivated 
by constraints from precision EW measurements. However as discussed earlier this is 
for illustration purposes only and in our full analysis we have kept 5^^ and as free 
parameters. 



In Fig. 20, we plot the bare cross section aQ{N4) with the basic cuts of Eq. ( 4.19 ) as well 
as the selection cuts Eqs. ( [4.21 ) — ( [4. 23| ) at 14 TeV. The solid (dotted) curves correspond 
to the bare cross section without (with) the Higgs decay channel for mn = 120 GeV. 
The reduction in rate is mainly due to the basic acceptance cuts. We note that the cross 
section with the cuts at 14 TeV is higher than that at 10 TeV by a factor of 1.4—1.6 for 
m4 = 100 — 500 GeV. The sensitivity reach for the mixing parameters to be presented later 
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Figure 19: Normalized distributions a-'^da/dX for m4 = 60, 100, 200 and 500 GeV at the LHC 
for (a) upper left: the minimal isolation AR^"^; (b) upper right: the missing transverse momentum 
^rp] (c) bottom left: the 2£2j system invariant mass m{££jj); (d) bottom right: the di-jet invariant 
mass m{jj). 



will be scaled down roughly according to this factor for LHC with cm. energy of 10 TeV. 

As discussed in the previous section, a large SM background comes mainly from top 
quark production and decay via the chain decay t ^ 6 — > c f^. Fortunately, after 
all the selective cuts in Eqs. ( 4.19| ) — ( [4.23| ), the top-quark decay background is essentially 



eliminated and has no remaining events for the expected luminosity of 100 fb~^ at LHC. 

There are several other SM backgrounds coming from like-sign W boson production 
at the LHC energies. First of all, the triple gauge-boson production process 



PP 



W^W^W^ l^l^vv jj, (4.24) 



leads to the irreducible background with two like-sign leptons plus jets. Next, the same 
final state can be produced via the process 



PP W^W^ jj l^e^uu jj, (4.25) 
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Table 5: The representative signal and background cross sections at the LHC, for ^Ji jj and 
the efficiencies with the consecutive cuts. For illustration, we have used 7714 = 200 GeV, |V;^4|^ — 
\Vri\^ = 5 X 10-3 > |Ve4|^ and = 120, 300 GeV 





No cut 


Basic cut 


cut 


+Ai? cut 


+m{jj),m{££jj) cut 








(4.23) 


dill) 


(I4.22D 


Signal 












a(fb) 


0.86 


0.42 


0.37 


0.35 


0.33 


eff. 


- 


48% 


88% 


96% 


94% 


ti Bkg 












a(fb) 


29.6 


16.9 


2.7 


0.075 


0.002 


eff. 




57% 


16% 


2.8% 


2.7% 




ruH =120 GeV 










a 


1.01 


0.42 


0.057 


0.052 


0.050 


eff. 


ruH =300 GeV 


42% 


14% 


91% 


96% 


a(fb) 


1.28 


0.58 


0.066 


0.061 


0.058 


eff. 




45% 


11% 


92% 


95% 




ruH =120 GeV 










a(fb) 


4.2 


1.3 


0.29 


0.17 


0.019 


eff. 


ruH =300 GeV 


31% 


22% 


59% 


11% 


a(fb) 


4.4 


1.4 


0.34 


0.19 


0.025 


eff. 




32% 


24% 


56% 


13% 



where the two jets may come from either QCD scattering or from the gauge-boson fusion 
process. However these backgrounds have two missing neutrinos and can be suppressed 
by a combination of cuts on the missing transverse energy and invariant mass. We also 
analysed the backgrounds coming from Z boson production 

PP jjZZ, pp jjZW. (4.26) 

in which some charged leptons are missing in the detection so that they lead to like-sign 
dilepton events. The backgrounds are very small after the cuts. 

We list the number of background events and efficiency of cuts in Table ^ for a lumi- 
nosity of 100 fb"^ at LHC. The total background is about 7 — 8 events for 100 fb~^ at the 
LHC. The main background is from the W^W^W^ channel and can be further suppressed 
if a tighter missing energy cut could be exploited. For instance, the background events may 
be reduced by half, leaving about 3 — 4 events with j^j. < 15 GeV. 

The last but not least important feature of the signal is the direct reconstruction of 
the resonant mass of ^4 in the final state i^jj- This is shown in Fig. 21 for the SM 
background and the signal with 7714 = 200, 400 GeV. We see the effective reconstruction 
of the resonant mass. For a given mass 7714 in the search, one can further make the event 
selection on m{ijj) 

0.8 m4 < m{ijj) < 1.2 ^4, (4.27) 

to estimate the significance of the signal observation. This loose cut has little effect on 
the signal, but reduces the total background to — 4 events for 100 fb"^ in the range of 
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10^ 




m4 (GeV) 



Figure 20: The bare cross section (Jo(-^4) versus heavy neutrino mass 7714 after all the cuts at the 
LHC (14 TeV). The solid (dotted) line correspond to the exclusion (inclusion) of the Higgs decay 
channel for mn = 120 GeV. 



m4 as shown in Fig. We once again adopt Poisson statistics to determine the search 
sensitivity. The number of signal events needed for 2a significance would be 3 — 11; and 
15 — 44 for 5a significance. In Fig. |2^(a) and Fig. |2^(b), we summarize the sensitivity 
for S^^ and Se^ versus m^, respectively. The solid (dashed) curves correspond to 2a {5a) 
limits on Su/ with the exclusion of the Higgs decay channel. The dotted (dash dotted) 
curves are similar but with the inclusion of the Higgs decay channel for niH = 120 GeV. 
The horizontal dotted line corresponds to constraints on |V^4p < 6 x 10~^ from precision 
EW measurements [^0|. In Fig. p3|(b) the dashed line at the bottom corresponds to the 
limit from Qv(3(5. 

I 1 2 I 1 2 I 1 2 '"i 

In the optimistic case, we assume that 11^41 <C |V^4| and S^^ ~ |V^4| < 6 x 10" . 
The detection sensitivity on heavy neutrino mass can be 

f 375 GeV for 2a; 

1 250 GeV for5<.. ^^'^^^ 
Or alternatively, the mixing parameter can be probed to 

^^^^{llZ-^ for 5a.' (^-^^^ 

In particular, even with the very stringent bound on |V^4p from Oz^/3/3 as indicated by 
the dashed curve in Fig. ^(b), one may still have 2a sensitivity if 7714 ~ mw- 

Our calculations for hadron colliders have been based on parton-level simulations. 



A recent study |33| pointed out that there may be other backgrounds to be considered 
when detector effects are included. One of them is the faked like-sign dileptons from 
the bh cascade decay. The other is due to the QCD multi-jet radiation to degrade the 
reconstruction oiW ^ jj. Those backgrounds can not be easily simulated in particular at 
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Figure 21: Invariant mass distributions of m(£jj) for the signal with 7714 = 200, 400 GeV and 
background processes. 



T 1 r 




(GeV) 

Figure 22: Number of background events vs mass of the heavy neutrino, m4. 

the parton-level. A preliminary analysis including full CMS detector simulations cannot 
support their claim |^4|. Nevertheless, we may consider to design more stringent acceptance 
cuts to further discriminate against the backgrounds. First, common wisdom suggests to 
tighten up the charged lepton isolation requirement 

AR'P/" > 0.8, (4.30) 

which would remove the backgrounds from b, c decays substantially, but a full assessment 
can be made only when real data become available and after the detectors are fully un- 
derstood. Next, we may increase the jet threshold to suppress the initial state QCD jet 
radiation to purify the W jj sample. Our estimate based on a PYTHIA simulation 
shows that the kinematics of a DY-type electroweak process can be largely preserved with 
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Figure 23: (a) Left: 2a and ha sensitivity for S*^^ versus m4 at the LHC with 100 fb^^ integrated 
luminosity; (b) right: same as (a) but for Se^^ ■ The solid and dashed (dotted and dash dotted) 
curves correspond to limits with the exclusion (inclusion) of the Higgs decay channel for niH = 120 
GeV. The horizontal dotted line corresponds to the constraint on S^^ ~ l^4p < 6 x 10^"^ from 
precision EW measurements pOl. 
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Figure 24: (a) Left: same as Fig. |g but with the ti ghter cuts of Eq s. (Ug) and ( |4.3l| ); (b) right: 
same as Fig. E3(a) but with the tighter cuts of Eqs. (4.30) and (4.31). 



the appropriate jet threshold. We thus examine the cut 

p>j, > 25 GeV, 



(4.31) 



which results in only about 17% of the events with potential jet contamination. The results 



with the tightened cuts are given in Fig. 24. We see that the stringent cuts severely hurt 



the low mass region, but the effect on the high mass region is modest. 
4.3 Like-sign Dilepton Signals with r in the Final States 

So far we have only presented the results with electron and muon final states and ignored the 
taus. This is due to the experimental challenge of r reconstruction. Given the importance 



-38- 



to cover all the lepton flavors, one must strive to include taus on the search list. Besides the 



known experimental practice for r identification at the Tevatron [108|, there are proposals 



to identify r events in connection with the neutrino sector [109]. The central issue is to 
reconstruct the missing momenta from r decays. We can generalize our requirement for 
the charged leptons to the isolated charged tracks presumably from the r decays (e, /i, or 
one-prong and three-prong charged hadrons) 

PT(track) > lOGeV, |7?(track)| < 2.5. (4.32) 

This assures that the parent taus are very energetic. For events with one r and no other 
sources of missing particles, the missing momentum will be along the direction of the 
charged track. We thus have 

p (invisible) = (track), (4.33) 

where the proportionality constant k is determined from the I^rp measurement by assigning 
= /tpT (track). For events with two taus, we generalize it to 

p (invisible) = Kip (tracki) + K,2P (track2). (4.34) 

As long as the two r tracks are not linearly dependent, ki and K2 can be determined again 
from the I^j, measurement. The missing momenta, as well as the r kinematics, are thus 
fully reconstructed. 

Although we believe that the A'4 signals in the modes of e^r^, /x^r^ and would 
be very promising for observation, the background analyses will be considerably more 
involved due to the complication of r reconstruction. Since our simulations are performed 
at the parton level, we are unable to adequately address the background suppression and 
to quantify the signal observability. We thus leave this for future studies. 



5. Summary and Conclusions 

The observation of a LV process would show that neutrino is a Majorana particle unam- 
biguously. Apart from light neutrinos, LV processes involving SM particles can receive 
a contribution from heavy Majorana neutrinos due to mixing. In fact, this contribution 
can be resonantly enhanced for appropriate masses of the heavy neutrino. In the absence 
of observation of LV interactions, the rates for these processes can constrain the mixing 
elements \Vi-^/iVi^ii\ as a function of the mass 7714 of the heavy Majorana neutrino. We 
considered two classes of LV violating processes: (a) low energy AL = 2 tau decays and 
rare meson decays and (b) collider signals for like-sign dilepton production with no missing 
energy implying the existence of Majorana neutrinos. We emphasize the necessity of in- 
volving two charged leptons and no neutrinos in the initial and final states, to be conclusive 
about lepton-number violation. 

For the low energy interactions we evaluated the transition rates and branching frac- 
tions as a function of the mass and mixing of the heavy neutrinos. We then translated the 
current experimental bounds from direct searches into limits on {Vi^^^Vi^^l as a function of 
the mass 7714 of the heavy neutrino. Amongst the rare meson decays, the ^^^^vr" 
decay mode currently gives the most sensitive experimental limits on \Ve-^4Ve^4\. Poten- 
tially, these constraints are six orders of magnitude more stringent than the constraints 
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from precision electroweak data which hmit |V£4p to few times 10~^. As the intermediate 
heavy sterile neutrino is a real particle which might exit the detector if the decay length is 
longer than the detector size, for very small mixing angles the bounds get weakened but 
arc still much more stringent than the electroweak precision constraints. This effect should 
be taken into account and a detailed analysis of past experimental data is required in order 
to find the precise limits on the mixing angles. Next in sensitivity are the D and Ds meson 
decay modes with constraints of the order of 10~^. Again, these are competitive with if 
not better than constraints from EW precision data. The other processes (in other mass 
ranges) have very weak experimental limits, weaker than EW precision data and essen- 
tially do not impose any meaningful bounds on {Vi-^^Vi^^l. This implies that more accurate 
experimental studies on those rare decays should be strongly encouraged. In particular, 
many interesting processes of D, B decays have not even been experimentally probed as 
well as those with a r lepton in the final state. Among the r-decay modes the best limits 
come from t~ £^-k~-k~ . The other r-decay modes have sensitivity of order 10~^ to 10~^. 
Again, the constraints from r decay modes are competitive with or better than constraints 
from precision EW data by 2 to 3 orders of magnitude. The experimental bound on LV 
processes is expected to improve in the future. The future sensitivity of the square of the 
mixing parameter will increase approximately by an order of magnitude for every order of 
magnitude improvement in experimental bounds on branching fractions. We have shown 
that the low energy AL = 2 r decays and rare meson decays can be very strong probes to 
discover or constrain the mass and mixing of heavy Majorana neutrinos. Even those decay 
modes which do not impose strong constraints should not be neglected. It only implies 
that a large range of the parameter space is available for exploration. 

In addition to analyzing the LV tau and meson decay modes and precision EW mea- 
surements we also compiled the constraints on the mixing elements (| 144^7 l^4p l^4p) 
from peak searches, accelerator experiments, reactor experiments and others - collectively 
called laboratory constraints. In the absence of detection of LV, the laboratory constraints 
and the ones from AL = 2 processes can be compared. The constraints on mixing from 
LV tau decays are always competitive with or better than laboratory constraints in the 
corresponding mass region while the constraints from LV meson decays are competitive 
with laboratory constraints only in some mass regions. We note that we explore more 
combinations of mixing elements and also provide better constraints on mixing in some 
mass regions. More importantly, a detection in one of the experiments analyzed to obtain 
laboratory constraints implies the existence of a sterile neutrino while a detection in one 
of the AL = 2 tau or meson decay modes studied in our analysis would imply LV and 
hence the existence of a Majorana neutrino. We pointed out the fact, often ignored in the 
literature when analyzing low-energy processes, that a heavy neutrino might decay outside 
the detector if it becomes long-lived for low mass (less than a GeV) and/or very small 
mixing. 

For the collider signals of heavy Majorana neutrinos we looked for the definitive lepton- 

number violating like-sign dilepton production and no missing energy. Such signals have 
low backgrounds and have the potential for discovery of heavy Majorana neutrinos. At the 
Tevatron, with the current and future integrated luminosities, we find that the mass of the 
heavy Majorana neutrinos can be probed up to 



m4 ~ 



{ 



10 - 130 GeV (10 - 75 GeV) for 2a {5a) with 2 fb"^ 
10 - 180 GeV (10 - 120 GeV) for 2a (5a) with 8 fb'^ 



(5.1) 
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Alternatively, the sensitivity for the mixing parameter can be 



2 / 2 X 10-5 (10-^) for la (5ct) with 2 fb'^ 
^M/^ ~ I l^Ml ~ I 5 X 10-6 (2 X 10-5) for 2a (5a) with 8 fb-\ ^^^^ 

This will surpass the DELPHI [|9| and L3 [|^, |ll^] 95% C.L. bounds. 

The sensitivity for heavy Majorana neutrinos can be extended significantly at the LHC. 
With 100 fb-^ of integrated luminosity, 

m4 ~ 375 (250) GeV for 2a (5ct); (5.3) 

or alternatively, the mixing parameter can be probed to 

V ~ |V^4p ~ 7 X 10-^ (3 X 10-6) for 2o (5cj). (5.4) 

The sensitivity at LHC will go well beyond the DELPHI and L3 95% C.L. bounds in both 
mass reach and mixing, and beyond the current bound on |V^4p from Ovf]/]. In summary, 
there is a rich avenue of possibilities for discovering or constraining the elusive Majorana 
neutrinos. 
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A. Lepton mixing formalism 

In this appendix, we illustrate the parameterization for the lepton sector, although we 
have not followed the relations literally, assuming that new physics beyond this minimal 
formalism exists. The leptonic content in the theory includes three generations of left- 
handed SM SU{2)i doublets and n right-handed SM singlets: 

i^aL=(^^) , N,R, (A.l) 

where a = 1, 2, 3 and 6 = l,2,3,---,n(n>2 for at least two massive neutrinos). 
The leptonically universal gauge interactions involving neutrinos are of the form 

-C= i^W+Y^v;;^ riaL + \i.c\ + -^Z^Y.^ ^""^-L- (A.2) 
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The gauge-invariant Yukawa interactions are 

(3 3 n 

^ fi, HkR + Y,Y1 fab L7l HN,r I + h.c. (A.3) 
a, 6=1 a=l b=l 

where H is the SM Higgs doublet and H = iT2H*. After the Higgs field develops a vev 
{H) — > \/2 , the Yukawa interactions lead to Dirac masses for the leptons 

(3 3 n \ 

Y,^^abhR + Yjl2^ <b NbR + h.c. (A.4) 
a,b=l a=l b=l J 

where the mass matrices are given by the vev times the corresponding Yukawa couplings 

= fabV/V2. 

The 3x3 mass matrix m' can be diagonalized by two unitary rotations among the 
gauge interaction eigenstates II, Ir 

Ol m} Or = diag(me, m^,m^), la = Oat ^, (A.5) 

where I = e, /i, r are the mass eigenstates, which define the charged lepton flavors. The 
Dirac masses as well as interactions with the Higgs boson for the charged leptons now have 
the standard form 

^ ¥f - 

-Cir = V mWl + — ) £ ^. (A.6) 
^-^ V 

If the Yukawa interactions of Eq. ( [A.3| ) are the whole source for neutrino mass, then we 
would have min(n, 3) massive Dirac neutrinos. 

To complete the neutrino mass sector, there is also a possible heavy Majorana mass 
term 

= o E Bbb' N,,R + h.c. (A.7) 



2 

b,b'=l 



where a charge conjugate state is defined as tp'^ = Cip'^ {if)^ = ijFC)^ and a chiral state 
satisfies {iP^)t = (^-r)^) with r = — , + for The full neutrino mass terms thus read 

^ I 'A n n \ 

-C=2 Y.T.('^<bNbR + N^mlu'aR)+Yl N^B,,, N,,R\+h.c. 

y a=l 6=1 b,b'=l J 



where we have used the identity UaL mab Njjr = N^^ rriba v^r-. 

The mass matrix can be diagonalized by one unitary transformation 



where the mass eigenvalues are of the order 



^^a,^"^. ML,^B. (A.IO) 
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L is a (3 + n) x (3 + n) unitary matrix and can be parameterized as 

\ ^nx3 J^nxn J 

The relation between the gauge interaction eigenstates and the mass eigenstates are given 

by 

'""^-ILf Z) , (A.12) 



with the mass eigenstates Um {m = 1,2,3), N^' {m' = 4, •••,3 + n). The diagonahzed 
(Majorana) mass terms of Eq. ( |A.8 ) thus read 



3+n 



-C = M E ^ + E ^^'L ^rn'R + h.C. , (A.13) 

\r?i=l m'=4 / 

with the mixing relations between the gauge and mass eigenstates 

3 3+ra 3 3+n 

VaL = Uaml^mL + Yl ^am'N^'L^ N^L = Y ^bmV^L + ^ ^bm'N^'L^ (A.14) 

m=l m'=4 m=l m'=4 

3 3+n 3 3+n 



Note that the unitarity condition for L leads to the relations 

UU^ + VV^ = U^U + X^X = h^^, (A.16) 
XX^ + YY^ = V^V + Y^Y = Inxn- (A.17) 

Parametrically, UU^ and Y^Y ~ VV'^ and X'^ X ~ ©(m^/MAr). 

In terms of the mass eigenstates, the gauge interaction lagrangian Eq. ( |A.2| ) can be 
written as 

/ r 3 r 3+n \ 

-C = -j=W; Yl E (U^OlU ^I^PlI + E E i^^OL)m'i N^J^Pl^ + h.C. 

\^=em=l ^=eni'=4 / 

/ 3 3+n 
+ IT^Z^ E iU^U)m,m, lA^l^PLiym, + Y (^^^)m[n.',N^r PlK 



^mi,m2=l m[,m'2=4: 
3 3+n 



2cosvy 



E E {U^VU^^m'^T^l^PLN^'^+i^-c. . (A.18) 



^mi=l m2=4 



To make the couplings more intuitive, we define the combination matrices by 

U^'^ = Olu, V^^ = Olv, U''^ = UW, U'"' = U^U, V^^ = VW. (A.19) 
We thus rewrite the gauge interaction lagrangian by one mixing matrix for each term 

(T 3 T 3+n \ 

E E t^i'Pl^ + E E ^^"PLi + h.C. 
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Figure 25: Feynman rules for the charged current vertices in terms of the neutrino mass eigenstates, 
as given in Eq. (A. 22). 



3+n 



2 cosw 



^mi,m2=l mj,m2=4 
3 3+n 



2 COSVF 



These coupHngs along with the mixing matrices Eq. ( A.19| ) give the most general leptonic 
interactions of the charged and neutral currents in terms of the mass eigenstates. Alter- 
natively, the neutral current interactions can be aligned along with that of the charged 
currents when rotating left-handed neutrinos in the same way as the charged leptons, 



3+n 



l^aL = iOL)ai l^iL, OT V^l = ^iOiU)e^UmL + E (^L^)^™' A^^'L • (^-21) 



m=l m'=4 



It may be convenient in certain practical calculations to rewrite the neutral current inter- 
actions in terms of their flavor eigenstates 

(T 3 r 3+n \ 

E E T^I^PlI + E E K^>1^PlI + h.C. 

l=e m=l l=e m'=A J 



+ V7£^^>^ ( E E ^^rn ^^Pl + E E y^rn' N^'l^ Pl l^i ] + h.C. + .(.A.22) 



r 3 T 3+n 



2 COS 9w 



-e m=l £=e m'=4 



where we have dropped the superscripts for U, V defined in Eq. ( |A.1S| ), for simplicity as 
adopted throughout the text. 

For the reader's convenience, we give most of the corresponding Feynman rules for the 
interaction vertices, listed in Fig. 25 for the charged currents, and in Fig. ^ for the neutral 
currents. The Feynman rules for the other diagrams can be easily deduced from the ones 
that are explicitly written down in Fig. |2^ and Fig. 

Finally, the heavy neutrino interactions with the Higgs boson read 

T 3+n 



= - E E ^fn'PL + h.C. + ... (A.23) 

l=e m'=4 

The corresponding Feynman rule for the interaction vertex is given in Fig. p^. 
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Figure 26: Feynman rules for the neutral current vertices in terms of the neutrino mass eigenstates, 
as given in Eqs. (A. 20). 




Figure 27: Feynman rule for the Higgs vertex in terms of the heavy neutrino mass eigenstates, as 
given in Eq. (A. 23). 

B. General amplitude of AL = 2 processes 

The charged current interaction lagrangian in terms of neutrino mass eigenstates is 

T 3 T 3+n 

-Ccc = -4^;(E E f^t* ^l^'PLi^Y. E ^l^' ^^'PLi) +11.C. (B.l) 

e=em=l e=em'=4 

where Pl = ^(1 — 75). The leptonic AL = 2 subprocess W~W~ — > £^£2 is induced by the 
product of two charged currents 

3 3+n 

Mf:, ^ E ^hmUt^m {hl^PLym){T2YPLiym)+ E Km'VlZ^' {hl>'PLNm'){T2Y PLNm'), 
m=l m'=4 

(B.2) 
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which can be rewritten using charge conjugation as 

3 3+n 
m=l m'=4 

(B.3) 

The Majorana neutrino fields can be contracted to form a neutrino propagator, and the 
transition matrix element is thus given by 



+ y E (^7^^l)^ ^^^^^^^^ in^Pn^W (B.4) 

m'=4 ^^m' ™ 

where (7 is the momentum exchange carried by the neutrino. The g' term vanishes due to 
the chirality flip. Including the crossed diagram {£1 <-> £2) the leptonic amplitude then 
becomes 

^Up = T E ^llmKm ^u^Wl-, + ) PrV2 

2 3+n 

+ 17 2^ l/^im' l^^am' "^A^w X 



2 ^ 

?Tt'=4 



^{-2 2 ' T + ^2 2 ' T (B.5) 

For the light Majorana neutrinos, namely, m = 1,2,3 the masses ~ 0(eV) 

and for the heavy Majorana neutrinos , the masses ^n^, ~ ©(MeV — GeV) for the low 
energy processes we consider. The heavy Majorana neutrino contribution has a resonant 
enhancement when q'^,q''^ ~ and is the dominant one. The light Majorana neutrino 

contribution however encounters a severe suppression due to the small neutrino mass like 
JM^. Hence we can neglect the contributions of the light Majorana neutrinos and 

the Y^^m=i P^^^ amplitude drops out. 

In principle all the heavy Majorana neutrinos will contribute to the amplitude but in 
our analysis we only consider the contribution of one of the heavy neutrinos, in particular 
the lightest one for simplicity. So the amplitude can now be written as 



H^p = T^^i4^^.4 m,-A- + '^^^ |P^^2. (B.6) 

' 2 \ - m| + ^^Ar4 777-4 q'^ - ml + iVN^mi 



We can rewrite the above amplitude as 

■^lep = irVe.AVe^A nii— + —Vi^iVt^i rm— 

2 q'^ - ml + 11^^1714, 2 q''^ - ml + iTj^^mA 

= Mi+M2. (B.7) 

When ~ 7774, Aii has a resonant contribution and when q'"^ ~ m\, M2 has a resonant 
contribution. In general, q ^ q' , and it is convenient to split up the individual resonant 
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contributions by the Single-Diagram-Enhanced multi-channel integration method | 112 |. To 
do this, define the functions 

IMP 



Then the amplitude squared is given by 

2 



(Bi 



(B.9) 



The amplitude squared splits up into the functions fi defined above and the phase space 
integration can be done for each fi separately. This helps to make convenient simplifications 
for the phase space integration and the computation can be carried out in parallel. The 
contributions from each fi can be added up after phase space integration. 



C. Decay modes of heavy Majorana neutrino 

In this section we will discuss in detail the decay modes of the heavy Majorana neutrino 
A^4, with mass much smaller than the mass of the W boson, mw. From EW precision 
measurements the mixing elements |V£4p < O(10~^) and the higher order terms in mixing 
would be very small and can be ignored. Hence the widths are presented only up to leading 
terms in mixing. The charged current and neutral current vertices of N4 with the mixing 



elements are given in Fig. ^ and Fig. ^ With increasing mass of the heavy neutrino new 
decay channels open up and can be classified into two body and three body decays. The 
decay width scales as the third and the fifth power of the mass(m4) for two and three body 
decays respectively. 

1) A'4 — > i~P~^ where £ = e, fi,T and is a charged pseudoscalar meson. This decay 



mode has charged current interactions only as shown in Fig. 25 and the decay width is 
given by 



r^p = r(iV4 ^ rp+) = -^fj> \Vgg,\^ \Vm\^ ml ii{fie,i^p) 



= [{1 + X - y){l + x) - ix]X^l,x,y), 
A(a, b, c) = a"^ + b'^ + - 2ab - 2bc - 2ca, (C.l) 



where fp is the meson decay constant and Vqq' are the CKM matrix elements, fii and fip 
are the masses scaled by the mass of the heavy neutrino and are given by = m1/m\. 

2) A'4 — > v^P^ where = Ve-,T^fi:T^T and P^ is a neutral pseudoscalar meson. This 
decay mode has neutral current interactions only as shown in Fig. ^ and the decay width 
is given by 

T-^p ^ r(iV4 ^ = ^fl |y,4p ml (1 - MP)^ (C.2) 

b47r 

where fp is the meson decay constant, fip \s the mass of the neutral meson scaled by the 
mass of the heavy neutrino and is given by /ip = m?p/m\. The mass of the light neutrino 
~ ©(eV) [^] is much smaller than the mass of A'4 ~ ©(MeV — GeV) and can be neglected 
to a very good approximation. We have set the mass of the light neutrino to zero in the 
expression for the width above and henceforth. 
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3) iV4 



'V~^ where i = e,iJ,,T and is a charged vector meson. This decay mode 



has charged current interactions only as shown in Fig. 25 and the decay width is given by 



r 



[(l + x-y)(l + x + 2y) - 4x]A^ (1, x, y), 



|V£4p m\ l2{fJ.e,iJ,v), 



(C.3) 



where /y is the vector meson decay constant and Vqqi are the CKM matrix elements, ni 
and nv are the masses of the lepton and the vector meson scaled by the mass of the heavy 
neutrino and are given by fii = mf/m\. 

4) ViV^ where = Ve^v^, and is a neutral vector meson. This decay mode 

has neutral current interactions only as shown in Fig. 26 and the decay width is given by 



h{x,y) = (l + 2y)(l-y)At(l,x,y), 



q2 



(C.4) 



where fv is the meson decay constant, /xy is the mass of the neutral meson scaled by 
the mass of the heavy neutrino and is given by /iy = rny/rn^. Ky is the vector coupling 
associated with the meson and is expressed in terms of = sin^ 9^, where 0^ is the 
Weinberg angle. The values of k for the various vector mesons are: k = hx^, for and lv; 



(4 + 

5) iV4 



2 ^1 



for K*^,K"^ and 6: and k 



fi - "^x 



for 0*^,0^" and J/V'. 



-=:*0 



^2 ^^^2 where ^1,^2 = e, ;U,t with li 7^ ^2- This decay mode has charged 



current interactions only as shown in Fig. 25 and the decay width is given by 



(1-^)2 



h{x,y,z) 



12 



ds , 



y^){l + z^-s)\-2{s,x\y^)\-2{l,s,z^), (C.5) 



(a;+j/)2 



where /i (0,0,0) = 1, Xj are the masses scaled by the mass of the heavy neutrino and are 
given by Xi = rrii/m/i. The mass of the light neutrino ~ 0{qY) is much smaller than the 
mass of A'^4 ~ ©(MeV — GeV) and hence can be neglected compared to the mass of A'4. 
We have set the mass of the light neutrino to zero with very good approximation in the 
expression for the width above and henceforth. 

6) A^4 — > vi^l'^'^t where ^1,^2 = e,//, r. Both charged current and neutral current 
interactions as shown in Fig. 25 and Fig. 26 are relevant for this mode and the decay width 
is given by 



96^ 



9l9r + Seie29R]h{xue^,xe2,xe. 



+ 



iaif + {dnf + ^iii2 (1 + 25!)) h (xu,^ , xe^ , xg^ 



/ds 1 1 

^(1 - s)X^{s,y\z^)X^{l,s,x^), 



(C.6) 
(C.7) 
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where /2(0, 0,0) = 1, Ii{x,y,z) has been defined in Eq. ( C.5 ), Xi are the masses scaled by 
the mass of the heavy neutrino and are given by Xj = mi/m^, = — \ + Xw, g\i = x^, and 
Xw = sin^ Oyj = 0.231, where Oyj is the Weinberg angle. 

7) A^4 Vi^vv where V(^^ = Ve^v^^^r- This decay mode has neutral current interactions 



only as shown in Fig. 26. Using the massless approximation for the neutrinos as described 



above the decay width has a simple form given by 

r-.,- = ^ TiN, ^ v,,v,,vi;) = ^m,,? ml (C.8) 

All the decay modes listed above contribute to the total decay width of the heavy 
Major ana neutrino which is given by: 



e,p ly e,p e,v 



where i,£i,i2 = e,ii,T. For a Majorana neutrino, the AL = process A^4 — > £~P^ as 
well as its charge conjugate |AL| = 2 process N4 i~^P~ are possible and have the same 
width, r^^. Hence the factor of 2 associated with the decay width of this mode in Eq. (C.9). 



Similarly, the AL = and its charge conjugate |AL| = 2 process are possible for the decay 
modes iV4 —>■ £~V'^ and — > £^£2^1^ and hence have a factor of 2 associated with their 
width in Eq. (|a9|). 

As mentioned earlier, new channels open with increasing mass of the heavy neutrino. 
For the low energy LV tau decays and rare meson decays we consider, the mass of the 
heavy neutrino is in the range 140 MeV ^ m4 5278 MeV. For this mass range we 
list all the possible decay channels for A'^4 in Table The mass and decay constants of 



pseudoscalar and vector mesons used in the calculation of partial widths given in Eqs. ( |C.l 
are listed in Table |^ in Appendix ^ 



D. Lepton- number violating tau decay 

The decay amplitude for lepton number violating tau decays can be separated into leptonic 
and hadronic parts, 

iM = {Mlep)^,u{Mhadr. (D.l) 

For the tree level amplitude, the hadronic part can be expressed in terms of the decay 
constants of the mesons in a model independent way. The box diagram includes hadronic 
matrix elements which cannot be simplified in terms of decay constants and needs to be 
evaluated in a model dependent way. We expect the tree level amplitude to dominate and 
do not include the box diagram. It has been argued that in certain cases for rare meson 



decays sub- leading contributions may be appreciable |113| , 20 1. Even in such a scenario 
the difference will not be important at the current level of sensitivities and we include the 
more conservative limit from tree level diagrams only. The tau decays and the rare meson 
decays are crossed versions of each other and the above arguments are true for both. 
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Table 6: Decay modes of heavy Majorana neutrino based on its mass 1714. 



Mass of heavy 
neutrino (MeV) 


Decay mode of 
heavy neutrino 


Mass of heavy 
neutrino (MeV) 


Decay mode of 
heavy neutrino 






> nip + nir = 1880 


A^4 fl^T'^Ur + C.C 
A''4 T^H^l/p + C.C 


> 2me = 1.02 


TV 7- -±- 

iV4 V(^e. 


> m^ + m^ = 1920 


A'4 T~-K~^ + C.C 


> nie + = 106 


A'4 fi~e~^Ue + c.c 


> me + niDs = 1970 


A'4 ^ 6^1)+ + C.C 


> m^o = 135 




> mp+ mo = 1980 


A'4 ^ + C.C 


> rUe + m-jr = 140 


N4 — > e~7r"'" + c.c 


> m£)*o = 2010 


A'4 ^ ^'£-0*'^ 


> 2m^ = 211 




> 771-^*0 = 2010 


— *u 
Af4 ^ zv^Ll 


> + m,r = 245 


A''4 ^ fl^TT^ + C.C 


> me + niD* = 2010 


N4 e-D*'^ + c.c 


> rrie + rriK = 494 


A^4 ^ e~K+ + c.c 


> mp +mDs = 2070 


A''4 fi'Df + c.c 


> mrj = 548 


A'4 ^ t'f?? 


> me + m£)* = 2110 


N4 e-D*+ + c.c 


> + rriK = 599 


A^4 fi-K+ + c.c 


> mp+ mo* = 2120 


A^4 ^ ^i-D*+ + c.c 


> nipO = 776 


N4 lygp" 


> mp+moi = 2220 


Af4 ^ fi'D*+ + c.c 


> nie + mp = 776 


JV4 — > e /? + c.c 


> m^ + m/^ = 2270 


A''4 ^ r"K+ + c.c 


> m^; = 783 


A'4 — > z/^u; 


> m-r- + mp = 2550 


A'4 — > T~ p'^ + C.C 


> nip + nip = 882 


A'4 ^ + c.c 


>mr + m*^ = 2670 


A^4 ^ r-K*+ + c.c 


> me + rriK* = 892 


iV4 — > e A ^ + c.c 


> m^^ = 298U 




> 171^*0 = 896 


A^4 ^ i^iK*^ 


> mj/^ = 3100 


N4 viJli^ 


> m^,o = 896 


N4 ^ i/^Z*" 


> 2m^ = 3550 


N4 V^^T^T^ 


> m^' = 958 


A'4 i^ev' 


> m-r- + mo = 3650 


N4 T-D+ + C.C 


> nip + mx* = 997 


N4 ^-K*+ + c.c 


> m-r- + m£)^ = 3750 


A^4 T-D+ + c.c 


> = 1019 




> mr + mo* = 3790 


Af4 ^ r-Z)*+ + c.c 


> nie + "Ir = 1780 


A'4 e^r+z^T- + c.c 
A^4 T~e+fe + c.c 


> mr + mn* = 3890 


Af4 ^ r-D*+ + c.c 


> me + mo = 1870 


A^4 e-D+ + c.c 







The leptonic part of the subprocess r ^ *VF * is obtained by crossing the 



ampUtude in (B.6) 



m4 



-ml + irNim4 

Combining the hadronic and leptonic parts, the decay amphtude for 
is given by 



iM = {Miep)^,uMl,^M'it^ + (Ml 



Mo 



(D.2) 

(D.3) 
(D.4) 



(Pi - gi)^ - m| + irAr4m4 
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M1+M2, 



(D.5) 
(D.6) 



where V^^^ are the quark flavor-mixing matrix elements for the mesons and /m^ are 
meson decay constants. Then the functions, fi and /2 defined in Eq. ( [B.8| ) are given by 



/i 

/2 



FrA 

4 + 1^ 



ial + b^)A + {aia2 + b^)C 



+ (91 ^ 92) 



(4 + + (af + b^)B 

fi{qi ^ 92), 

MPi^Qj) = ^{Pi ■ Qi){P2 ■ q2){qi ■ 92) - 4m|,j^(pi • q2){p2 ■ q2 

- ^'rnlhiPi ■qi){P2-qi) + 2m|f^m|^2(pi -pa), 
B{pi,qj) = A{qi ^ ^2), 
C{pi,qj) = 4(pi •P2)(gi -92)^ -A{pi,qj), 
C{qi ^ 92), 



D{pi,qj) 

Fr 



MiJM2\ 



CKM^rCKM\2\ 
V I 



M2 



VriVul'^ml 



(Pi - qi,2f - ml; b = Tn^tji^. 
The decay width for the LV tau decay is then given by 



fidPSsi + / f2dPS: 



IT 



dPS32 = dPSsiiqi ^ q2), 



'32 



''^^''^-dyidy2dyzdyi 



m 



cl 



(D.7) 
(D.8) 

(D.9) 
(D.IO) 
(D.ll) 
(D.12) 
(D.13) 
(D.14) 



(D.15) 

(D.16) 
(D.17) 



where dPS^i and dPS'32 are the phase space factors obtained by conveniently clustering 
two different sets of particles to enable applying the narrow-width approximation easily. 
yi to 7/4 are rescaled angular variables with integration limits < < 1. As seen in 
Sec. p.l.l| , the width of the heavy neutrino is very small compared to the mass and hence 
we can apply the narrow-width approximation. 



/ 



(m; 



m 



riV4^ 



S{m1_ — m\)dw?^. 



(D.18) 



Applying the narrow-width approximation as described above and integrating over the 
(5— function we get 



fidPSsi 



FrA-K^ 

,2 2 2 



{4 + b'^)A + b'^{B + C + D) 
{4 + b'^)A + b'^B 



I 



X2 {m^, mij^ , mi)\'2{ml, mj,mlj^)j dyidy2dyzdyA, 
f2dPS32 = J fidPSsiiqi ^ q2) 



(D.19) 
(D.20) 
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Now, we can find the decay rate from Eq. ( D.15 ). Normalized to the r decay width 
r,- = G'jpm^/192iT'^, the corresponding branching fraction is Br = TJ^/Tt-. The masses and 
decay constants of mesons are hsted in Table 0. The CKM matrix elements and r mass 
are taken from the Particle Data Group (PDG) [p9[: 



1777 MeV, \Vud\ = 0.9738, |Ks| = 0.2200. 



E. Rare meson decay 

The rare meson decays 

have the same Feynman diagrams as tau decay. The meson M2 can be a pseudoscalar or 
vector meson. The decay amplitude when M2 is a pseudoscalar meson is given by 



(E.l) 



Next we consider the case where M2 is a vector meson. The decay amplitude is given by 



= MX +mX. 



{qi -piY - 77142 + iTAT^ 777-4 
{qi-V2Y - niA? + 7rAr4m4 



(E.2) 



Similar to tau decay, we define the functions /f , f2 and fY , fY for the pseudoscalar 
and vector mesons respectively as given in Eq. ( |B.8| ). ff and fY turn out to have the same 
form and are given below. 



fPiV) 

Jl 




[ a? + 62 


fP{V) 


= fi [Pi ^ 




= 8{pi ■ qi){p2 



{4 + b^)AP(^) + iaia2 + b^)C 



P{V) 



{4 + 62)A^(V) + (a2 + 62)5P(y) 
{qi ■ q2) - 4m|^^ {pi ■ 92) (P2 • ^2) 

- "^l^lhiPl ■ (ll){P2 ■ qi) + '^irLMr^lhiPl ■P2), 

A^(.Pi,qj) = 8{pi ■ qi){p2 ■ q2){qi ■ 92) - 47nif^(pi • q2){p2 ■ ^2) 
+ ^iT^lhiPi ■qi){P2-qi)- '^iriMrrnlf^ipi ■ P2), 

C^(.Pi,qj) = 4(pi ■P2){qi ■q2f -A^{pi,qj), 

C^(.Pi,qj) = ■P2){qi ■ q2f - '^T^lh'mlj^ipi ■P2) - A'^{pi,qj), 



+ {pl 



P2) 



,(E.3) 
(E.4) 

(E.5) 

(E.6) 
(E.7) 
(E.8) 
(E.9) 
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Fm = 4:Gpflf flf \Vj 



MiJM2 \ ^Mi 



CKM^rCKM\2 
V 1 



\Ve^iVi^i\'^ml, 



M2 



ai,2{Pi,qj) = (qi - Pi,2f - ml; b = TNim4^. 
The decay rate for LV rare meson decay is then given by 



LV 



dPS: 



31 



,2 _2 _2 \\h^^2 _2 _2 \dm^i 



dyidy2dyzdyi 



cl 



(E.IO) 
(E.ll) 
(E.12) 



(E.13) 

(E.14) 
(E.15) 



where dPS^i and dPS^2 are the phase space factors obtained by conveniently clustering 
two different sets of particles to enable applying the narrow-width approximation easily. 
yi to ?/4 are rescaled angular variables with integration limits < < 1. The width 
of the heavy neutrino is very small compared to the mass and hence we can apply the 
narrow-width approximation. 



/ 



dm^. 



m 



m 



iy+r%mi 



Tn^-^ 



(5(m^. — ml)dm'l. 



(E.16) 



Applying the narrow-width approximation as described above and integrating over the 
(5— function we get 



/ f^'dPSs, 



( FmA^v)^^ 



N4 



(a2 + b^)AP(v) + b^BP(y) 



f^^''^dPSs2 



A 2 (niM^ , mj^ , 7714) A 2 (ml, m\,m\j^)j dyidy2dy2,dyi, (E.17) 

j f^^''UPS^i{pi^P2), (E.18) 
Now, we can find the decay rate from Eq. ( [E.13 ). The branching fraction is then given by 

The CKM matrix elements and the lifetimes of mesons used in our calculations are 
taken from PDG [99| and are listed below. 



iKfel = 0.00367, l^edl = 0.224, \V,, 



0.996; 



TK = 1.2384 X 10^*^ s, TD = 1.040 x 10"^^ ^ 4 9 ^ iQ-is s, tb = 1.671 x 10"^^ g 

The mass and decay constants of pseudoscalar and vector mesons used in our calculations 
are listed in Table |^. 
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Table 7 : Mass and decay constants of pseudoscalar and vector mesons used. 



Pseudoscalar 
Meson 


Mass 
(MeV)[|9l 


Decay Constant 
(MeV) @ 


Vector 
Meson 


Mass 
(MeV)@ 


Decay Constant 
(MeV)[m 


-r 


139.6 


130.7 


-r 


775.8 


220 




493.7 


159.8 




891.66 


217 




1869.4 


222.6 


m 


D* 


2010 


310 




1968.3 


266 




2112.1 


315 


B± 


5279 


190 [] 




UJ 


782.59 


195 




135 


130 




896.10 


217 


V 


547.8 


164.7 


117] 


<P 


1019.456 


229 




957.8 


152.9 


117] 




2006.7 


310 


Vc 


2979.6 


335.0 






3096.916 


459 [US] 
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Abstract: The Majorana nature of neutrinos can be experimentally verified only via 
lepton-number violating processes involving charged leptons. We study 36 lepton-number 
violating (LV) processes from the decays of tau leptons and pseudoscalar mesons. These 
decays are absent in the Standard Model but, in presence of Majorana neutrinos in the 
mass range ~ 100 MeV to 5 GeV, the rates for these processes would be enhanced due to 
their resonant contribution. We calculate the transition rates and branching fractions and 
compare them to the current bounds from direct experimental searches for AL = 2 tau and 
rare meson decays. The experimental non-observation of such LV processes places stringent 
bounds on the Majorana neutrino mass and mixing and we summarize the existing limits. 
We also extend the search to hadron collider experiments. We find that, at the Tevatron 
with 8 fb~^ integrated luminosity, there could be 2a {5a) sensitivity for resonant production 
of a Majorana neutrino in the jJi^jJir^ modes in the mass range of ~ 10 — 180 GeV (10 — 
120 GeV). This reach can be extended to ~ 10-375 GeV (10-250 GeV) at the LHC of 14 
TeV with 100 fb"^ The production cross section at the LHC of 10 TeV is also presented 
for comparison. We study the //^e^ modes as well and find that the signal could be large 
enough even taking into account the current bound from neutrinoless double-beta decay. 
The signal from the gauge boson fusion channel W^W'^ ^i^t LHC is found to be 

very weak given the rather small mixing parameters. We comment on the search strategy 
when a r lepton is involved in the final state. 
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1. Introduction 

In the Standard Model (SM) of strong and electroweak interactions, neutrinos are strictly 
massless due to the absence of right-handed chiral states {Nn) and the requirement of 
SU{2)l gauge invariance and renormalizability. Recent neutrino oscillation experiments 
have conclusively shown that neutrinos are massive j^] . This discovery presents a pressing 
need to consider physics beyond the SM. It is straightforward to obtain a Dirac mass 
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term mDipiNji + h.c.) for a neutrino by including the right-handed state, just hke the 
treatment for ah other fermions via Yukawa coupUngs to the Higgs doublet in the SM. 
However, a profound question arises: Since Nr is a SM gauge singlet, why should a gauge- 
invariant Majorana mass term ^MNrNr not exist in the theory? In fact, there is strong 
theoretical motivation for the Majorana mass term to exist since it could naturally explain 
the smallness of the observed neutrino masses via the so-called "see-saw" mechanism [01 



(1.1) 

From a model-building point of view, there are many scenarios that could incorporate the 
Majorana mass. Examples include Left-Right symmetric gauge theories |^]; S'O(IO) Su- 
persymmetric (SUSY) grand unification Q and other grand unified theories models 
with exotic Higgs representations [^, ^; R-parity violating interactions (AL = 1) in Super- 
symmetry (SUSY) § and theories with extra dimensions |^]. There are other proposals 
to generate Majorana masses for neutrinos at a higher scale M without relying on the 



right-handed state [10, 11 1. According to the scheme in generating the mass scale M 



in Eq. (|1.1| ), it has been customary to call them Type I Q, Type H [|T^ or Type HI [pT]]. 



Within the context of the SM, there is only one gauge-invariant operator [12| that is 
relevant to the neutrino mass, 

IlH, (1.2) 

where and H are the SM lepton and Higgs doublets, respectively. The constant k is a 
model-dependent effective coupling and A is the new physics cut-off scale. It is a dimension- 
5 non-renormalizable operator, and leads to Majorana neutrino masses of the order kv'^ / A, 
after the Higgs field acquires a vacuum expectation value v, in accordance with the see-saw 
scheme. Higher dimensional operators that give rise to Majorana neutrino masses have also 
been constructed in a model-independent manner [^]. The challenging task is to look for 
experimental evidence to probe the new physics scale A and to distinguish the underlying 
theoretical models mentioned above. 

In the neutrino sector, besides the rich phenomena of neutrino flavor oscillations and 
the possible existence of new sources of CP- violation, lepton number violation by two units 
(AL = 2), as implied by a Majorana mass term, plays a crucial role. Not only may it 
result in important consequences in particle physics, nuclear physics and cosmology but 
it would also guide us in understanding the fundamental symmetries of physics beyond 
the SM. Although the prevailing theoretical prejudice prefers Majorana neutrinos, exper- 
imentally testing the nature of neutrinos and lepton-number violation {LV) in general, is 
of fundamental importance. In accelerator-based experiments, neutrinos in the final state 
are undetectable by the detectors, leading to the so-called "missing energy" and therefore 
missing lepton numbers as well. One is thus forced to look for charged leptons in the final 
state. The basic process with AL = 2 can be generically expressed by 

W-W~ ^ £^q, (1.3) 

where W~ is a virtual SM weak boson and £1^2 = e, /U, r. By coupling fermion currents 
to the W bosons as depicted in Fig. [l|, and arranging the initial and final states properly, 
one finds various physical processes that can be experimentally searched for. The best 
known example is neutrinoless double-beta decay (Oz^/3/3) [|l^, |l^, |l^, which proceeds via 
the parton-level subprocess dd — > uu W~*W~* — > uu e~e~. Other interesting classes 
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Figure 1: A generic diagram for AL = 2 processes via Majorana neutrino exchange. 



of LV processes involve tau decays such as t~ — > £~^M^M2 
mesons Mi,M2 are Tr,K, rare meson decays such as 
hyperon decays such as S~ — > T,~^e~e~, E~ — > pfi~ ^~ etc. [21]. One could also explore 
additional processes like e" 



1 13, where the light 
H, 10, [H and 



[|35| , charged-Higgs production e 




neutrino 
top- 
6|, and in 



, /i — > e"*" ||18|, and fi ^ conversion 
One may also consider searching for signals at accelerator and collider experiments via 
e'e- W~W- [H], e+e- ^ N + X 

nucleon scattering vi(j7i)N PtflfX 
quark decays t 
the decay N ^ i^H^ ||]. 

The dynamics for AL = 2 processes as in Eq. ( |1.3D is dictated by the properties of the 
exchanged neutrinos. For a Majorana neutrino that is light compared to the energy scale 
in the process, the transition rates for LV processes are proportional to the product of two 
flavor mixing matrix elements among the light neutrinos and a LV mass insertion 



I2 



m=l 



(1.4) 



where {m)g_^g^ are the "effective neutrino masses" . If the neutrinos are heavy compared to 
the energy scale involved, then the contribution scales as 



3+n 



(1.5) 



where V is the mixing matrix between the light flavor and heavy neutrinos. Unfortunately, 
both situations encounter a severe suppression either due to the small neutrino mass like 
m1^/M^r, or due to the small mixing | Vgjm'V^jm'P- An important observation is that 
when the heavy neutrino mass is kinematically accessible, a process may undergo resonant 
production of the heavy neutrino. The transition rate can be substantially enhanced and 
goes like 

T{N^, ^ i) T{N^, ^ /) 



m 



N„ 



N 



(1.6) 
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where i, f refer to the initial and final state during the transition. 

The possible existence of sterile neutrinos in the mass range relevant for resonant en- 
hancement of AL = 2 processes studied in this paper is motivated in several scenarios. 
Models which implement the see-saw mechanism at low energies have been recently con- 
sidered the neutrino masses generated are accidentally small and active-sterile 
mixing can be as large as few percent. See-saw models at the electroweak scale can explain 
neutrino masses if appropriate symmetries are imposed and at the same time provide an 
appealing mechanism for baryon asymmetry generation via resonant leptogenesis In 
theories with dynamical electroweak symmetry breaking, sterile neutrinos with masses in 
the 100s of MeV to GeV range are invoked to explain light neutrino masses |41]. Sterile 
neutrinos can also play a role in understanding the flavour problem in the leptonic sector. 
It has been shown that mixing with sterile neutrinos can be at the origin of the large angles 
in the neutrino sector |42]. 

Heavy, mostly-sterile neutrinos have been investigated for their role in cosmology 
and astrophysics, in particular in Big Bang Nucleosynthesis, Large Scale Structure for- 
mation [43 1, cosmic microwave background, diffuse extragalactic background radiation, 
supernovae [44| and as dark matter candidates |4^, (for a review on MeV ster- 
ile neutrinos, see Ref. [48|). A keV sterile neutrino is a viable dark matter candidate 
1 45, which can also explain the origin of pulsar kicks [49|. Decays of heavy, mostly- 
sterile neutrinos have been proposed to explain the early ionization of the Universe |^0|| . 
Due to mixing, dark matter sterile neutrinos would decay radiatively contributing to the 
Diffuse Extragalactic Background Radiation and inducing x-ray emission from galaxy clus- 
ters [51, 52 1 . A large coupling between sterile neutrinos and light dark matter scalars can 
be at the origin of neutrino masses and of the observed dark matter abundance [^]. A 
model with sterile neutrinos in the keV-GeV mass range has been proposed to explain the 
dark matter of the Universe as well as baryogenesis [54, 55 [. Its phenomenological and as- 
trophysical signatures have been considered in detail in Refs. [56, 52[. This model assumes 
the existence of one sterile neutrino with keV mass for dark matter and two heavier neu- 
trinos with quasi-degenerate GeV masses for successful baryogenesis. The required mixing 
\Vim'\'^ of the latter neutrinos with the active ones is mass dependent and lies in the range 
10^^^ — 10~^, for a mass of 1 GeV. Additional constraints on the heavy neutrino mass and 
mixing angles can be derived from astrophysical observations. Sterile neutrinos mixed with 
active ones would be efficiently produced in supernovae cores, escaping from it and deplet- 
ing substantially the supernova core energy, and, therefore, might modify the supernova 
evolution. Recently, it was shown that sterile neutrinos in the mass range ~ 0.2 GeV and 
small mixing angle with /i and r neutrinos could enhance the energy transport from the 
core to the stalled shock and favor the supernova explosion They could also explain 
the high velocity of pulsars if the momentum carried away by heavy sterile neutrinos is 
emitted asymmetrically [^^. Detailed reviews and discussions of heavy neutrinos in the 
Early Universe and their present bounds can be found, e.g., in Refs. 



I Mm- 

Cosmological and astrophysical constraints on sterile neutrinos are typically very strong 
but are not as robust as the ones from laboratory searches as they typically depend on the 
production mechanism of sterile neutrinos in the Early Universe and on the cosmological 
evolution. For example, they can be significantly weakened or evaded if the reheating 
temperature is low [47, |5^, if their density in the Early Universe is diluted by entropy 
injected at late times [^] or if they have non-standard interactions. In these cases, much 
larger mixing angles with active neutrinos are allowed by cosmological observations and 
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can be tested in terrestrial experiments. Therefore, it is important to perform experimental 
searches of heavy sterile neutrinos with increased sensitivity and, specifically for Majorana 
neutrinos, to consider AL = 2 processes. If a positive signal is found and is incompati- 
ble with the cosmological and/or astrophysical observations, one would need to consider 
modifications to the standard cosmological scenario and/or would gain new insight on the 
evolution of astrophysical objects. 

In this paper, we study resonant contributions of heavy Majorana neutrinos to AL = 2 
processes involving two charged leptons in accelerator-based experiments. We establish 
our conventions and discuss the current constraints on the mass and mixing of heavy 
neutrinos in Sec. ^. In Section ^ we lay out the general expressions for the heavy neutrino 
contributions to low energy LV decays and study two classes of AL = 2 processes, 

(a) tau decays r~ ^ M2 , 

(b) rare meson decays K+, D+, Df, B+ ^ M2 ■ 

We calculate the enhanced transition rates and branching fractions and compare them 
to the bounds set by direct experimental searches. A non-observation of such AL = 2 
processes places stringent constraints on the mass and mixing of Majorana neutrinos which 
are also presented in this section. The resonant production of Majorana neutrinos at hadron 
colliders, namely the Tevatron and LHC are studied and updated in Sec. |^. We draw our 
conclusions in Sec. We discuss in detail the formalism, the decay modes and the total 
decay width of heavy Majorana neutrinos and the transition rates of LV processes in the 
Appendices. 

2. Majorana neutrinos in extension of the standard model 

To set up our notation and convention, we first discuss the formalism for the simplest 
extension of the SM which includes right handed singlets. Also in this section, we present 
the current constraints on the mass and mixing of a heavy neutrino from various direct 
detection experiments, accelerator searches and electroweak precision constraints. 

2.1 Formalism for Heavy Neutrino Mixing 

The leptonic content in our simplest extension of the SM includes three generations of 
left-handed SM SU{2)l doublets and n right-handed SM singlets 



LaL=\^{'j , NbR, (2.1) 

where a = 1,2,3 and b = l,2,---,n. The gauge-invariant Yukawa interactions lead to 
Dirac masses for the charged leptons and neutrinos after the Higgs field develops a vacuum 
expectation value v. It is also possible for the singlet neutrinos to have a heavy Majorana 
mass term. The full neutrino mass terms as well as the diagonalized eigenvalues can be 
expressed as 

= ^ I E E <b Nm + ML <: <r) +T. Bw Nb,R 1 + h.c. 

b,b'=l J 
3+n \ 

^ v^i^-^^+Yl Mn^, N„,r + h.c. (2.2) 

\m=l m'=4 / 
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with the mixing relations between the gauge and mass eigenstates 



3 3+n 

^ Vam'N'^,^, (2.3) 

m=l m'=A 

UU^ + VV^ = I. (2.4) 
In terms of the mass eigenstates, the gauge interaction Lagrangian can be written as 

(r 3 r 3+n \ 

J2 E t^^'Pl^ + E E N^rPLi + h.c. 
e=e m=l i=e m'=4 / 

(r 3 T 3+n \ 

E E ^^rn ^I^Pl + E E ^T'^^L ^£ + h.c.(2.5) 

£=e m=l ^=e m'=4 / 



g_ 

2 cos 6*^/ 



Further details about the mixing formalism are given in Appendix 

A few important remarks are in order before the detailed considerations. First of all, 
parameterically, the light neutrino masses rriy^ diag are of the order of mag nitude {m^^f/B, 
while the heavy neutrino masses are Mjv^ diag — B. Secondly, the mixing parameters would 
typically scale as C/Tf/ ~ / and V^V ^ m^/MN. Thirdly, the Majorana mass term for the 
flavor states Uai, absent in Eq. ( ^I^ ) and corresponding to the null entry 03x3 in Eq. ([A.^ ), 
may receive non-zero contributions as Majorana masses for the light active neutrinos, for 
instance from higher dimensional AL = 2 operators or in theories with a triplet Higgs 
field. The general formalism presented here remains the same. In this paper, we will take a 
phenomenological approach toward these parameters. We will simply take the masses and 
mixing elements of the heavy neutrino as free parameters, only subject to some constraints 
from experimental observations. The assumption that the masses and mixing elements 
are not rigorously related by the see-saw relations is feasible from a model-building point 
of view, since some fine-tuning or some ansatz of the neutrino mass matrix can always 
alter the general relations. Several scenarios where it is possible to have rather low mass 
of the heavy neutrino were mentioned in the previous section. Here and henceforth, we 
consider the case when only one heavy Majorana neutrino is kinematically accessible and 
denote it by with the corresponding mass and mixing with charged lepton flavors 
Vi4. If we stick with this simple parameterization, the SM Higgs boson will couple to the 
heavy neutrinos as well. We present the couplings in Appendix |A[ When appropriate, we 
will include this effect. As noted above, some fine-tuning would be needed to avoid 
excessive contributions to the light neutrino mass. 

2.2 Current Constraints on N4 Masses And Mixing 

In laboratory searches, no positive evidence of sterile neutrinos has been found so far,^ in 
the mass range of interest, 100 eV-100 GeV.^ 



^Indications of the existence of a neutrino with 17 keV mass were subsequently shown to be non valid. 
For a review, see Ref. |6l| . Studying interactions of neutrinos from it and fi decays, an anomaly in time 
distribution was found [p2|. It could be interpreted as the existence of a neutrino emitted in pion decays 
with mass of 33.9 MeV. Searches for this neutral fermion have not given any positive signature and have 
allowed to constrain the mixing to be |V^4p < 9.2 x 10"® at 95% C.L. 

^For sterile neutrinos with smaller masses a rather complete analysis of the bounds can be found in 
Ref. . See also the implications of the recent results from the MiniBooNE collaboration ^ . 
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A very powerful probe of the mixing of heavy neutrinos with both and z/^ are peak 



searches in leptonic decays of pions and kaons |67]. If a heavy neutrino is produced in 



such decays, the lepton spectrum would show a monochromatic line at 



2mM 



where and mi are respectively the lepton energy and mass, niM is the meson mass. The 
mixing angle controls the branching ratio of this process as: 



T{M+^l+Ni) |y,4|' 



V^mIV , (2.7) 



where p \s a. kinematical factor 



1 + /if + /i| - 2(/i^ + /U4 + pi^i) [fie + IJ.4- - PaY 
P = ^ 7, 7^ ^' (2.^ 

with Hi = mf/mj^. For large m^, the helicity suppression of the tt,K ^ iug decays weakens 
and there is an enhancement for M"*" — > l'^ N/^ by a relative factor ml/mf, reaching up to 
10^ — 10^ compared to that of vr — > eue and K eu^ in the SM, respectively. These bounds 
are very robust because they rely only on the assumption that a heavy neutrino exists and 
mixes with and/or i/^. 

Another strategy to constrain heavy neutrinos mixed with i/g, and v-r, is via searches 
of the products of their decays. If kinematically allowed, A^4 would be produced in every 
process in which active neutrinos are emitted with a branching fraction proportional to 
the mixing parameter |V£4p. They would subsequently decay via Charged Current (CC) 
and Neutral Current (NC) interactions into neutrinos and other "visible" particles, such 
as electrons, muons and pions. Searches for these "visible" decay-products were performed 
and were used to constrain the mixing parameters. In beam dump experiments, A^4 would 
be produced in meson decays, with the detector located far away from the production 
site. The suppression of the flux of A^4 needs to be taken into account if the decay length 
is very short and, therefore, typically both an upper and a lower bound on the mixing 
angle can be set. Otherwise, the production can happen in the detector itself, as for the 
limits obtained from a reanalysis of LEP data, using the possible decays of the [26| 
into heavy neutrinos. In this case, large values of the mixing angle can be excluded. 
These bounds are less robust than the ones previously discussed. In fact, if the heavy 
neutrinos have other dominant decay channels into invisible particles, these bounds would 
be weakened, if not completely evaded. For example, a coupling of the type gNvcf) (see 



Ref. |68, |5^), with (/> a scalar, can induce very fast decays, which might dominate over the 
ones induced by CC and NC interactions. In this case, if the decay length is very short due 
to these strong interactions, the flux of A'4 might be suppressed at the far detector and the 
bound would not apply. If the production happens in the detector itself, the bounds would 
need to be recomputed considering the branching fraction into "visible" channels. Notice 
that we do not report here the bounds from Ref. as they do not apply to the heavy 
neutrinos under consideration. In these analyses it was assumed that heavy neutrinos were 
produced via — > N4^N4^ with the same strength as an active neutrino. In our scheme, 
this would correspond to having mixing angle equal to 1. Then, the search for A^4 decays 
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in the detector was used to constrain the heavy neutrino parameters. These data should 
be reanalyzed considering that the production of A'4 is suppressed by |V£4p. Comparing 
the expected number of events with the backgrounds we estimate that typically bounds of 
order |V£4p < few x 10~^-10~^ could be deduced. However, a detailed analysis should be 
performed and we do not report these limits in our figures. 

For masses above the production threshold, additional constraints can be obtained 
from lepton universality as the decay rates for muons, pions, taus as well as the invisible 
decay width for the boson are modified with respect to the SM predictions ||7^, |72| . 
Flavour changing neutral current processes such as /i — > 67, ee'^e~ and ^u-e conversion 
in nuclei are affected by the existence of heavy sterile neutrinos and strong limits can be 
obtained on the mixing with active neutrinos [^, 74, 75]. These bounds are reported in 
Section p. 2.4 



Finally, in Section 2.2.5 we discuss the very strong constraints on |T44p which can be 



obtained from the non-observation of neutrinoless double beta decay. It should be noticed 
that in the presence of more than one sterile neutrino, possible cancellations between the 
contributions to the decay rate can be achieved and the bounds would be consequently 
much weaker. 

Next, we review the laboratory constraints on the mixing between heavy and active 
neutrinos, depending on flavour and the mass of sterile neutrinos. 

2.2.1 Mixing with i/g 

The mixing parameter Ve4 can be tested in searches of kinks in the /3-decay spectrum, of 
peaks in the spectrum of electrons in meson decays and, finally, of iV4 decays in reactor 
and accelerator neutrino experiments. 

For masses 30 eV ~ 7714 ~ 1 MeV, the most sensitive probe is the search for kinks in 
the /3-decay spectra In the presence of heavy neutrinos mixed with Ug, the Kurie plot 
would be given by the contributions of the decays into light neutrinos as well as into heavy 
ones. This induces a kink in the Kurie plot at the end point electron energy Eg 

2Mi ' ^ ' 

where Mj j are the mass of the initial and final nuclei, respectively, and mg is the electron 
mass. In Fig. |2| we report the most stringent limits, obtained by using different nuclei 
1 76, 77, 78, |8^. In reactors and in the Sun only low mass, mi < few MeV, heavy 



neutrinos can be produced. The constraints obtained by looking for their decays into 
electron-positron pairs |8^ are reported in Fig. ^ with solid (cyan) contour labeled 
Bugey and short dashed (blue) contour labeled Borexino. The region with long dash 
dotted (grey) contour, labelled vr eu, is excluded by peak searches 



For heavier masses peak searches give the most stringent bounds, shown in Fig. |3|. 
Notice that, due to the weakened helicity suppression of the tt decay, the sensitivity on 
Vei increases with m4 till phase space becomes relevant at 7714 > 80 MeV, for vr efg- 
The excluded region, at 90% C.L., from Ref. [^ ], is indicated with the solid (black) line 
labeled tt ev. For heavier masses, stringent bounds are obtained by looking at the 
electron spectrum in K decays |8^] and are indicated by the double dash dotted (purple) 
line labeled K ^ ev \w Fig. Assuming that only CC and NC interactions are at play, 
stringent bounds have been obtained on |V^4p and are reported in Fig. ^ by the rest of the 
contours (except dotted (maroon) line labeled Oz^ /?/?). In particular, the limits at 90% C.L. 
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Figure 2: Bounds on |Ve4p versus 7714 in the mass range 10 eV-10 MeV. The excluded regions 



with contours labeled ^^^Re |76|, , ^^Ni |7| , Q , ^Op and Fermia |80[1 refer to the 

bounds from kink searches. All the limits are given at 95% C.L. except for the ones from Ref. ||8^ 
which are at 90% C.L.. The areas delimited by short dashed (blue) contour labeled Borexino and 
solid (cyan) contour labeled Bugey are excluded at 90% C.L. by searches of A'4 decays from the 
Borexino Counting Test facility |81 and Ref. |Q respectively. The region with long-dash-dotted 
(grey) contour, labelled tt eu, is excluded by peak searches Q. The dotted (maroon) line labeled 
Oi^PP indicates the bound from searches of neutrinoless double beta-decay |83]. 



from Refs. 87, 88], assume the production of in meson decays and look for visible 
channels in a detector located some distance from the source. The limits at 95% C.L. in 
Refs. []89| , |90t| analyse the data from DELPHI and L3 detectors, looking for from Z^- 
decays. In Fig. we also report the excluded region from neutrinoless double beta-decay 



experiments |91, bounded by dotted (maroon) line, valid if the heavy neutrinos are 



Majorana particles (see further). 
2.2.2 Mixing with 

The bounds on |V^4p come from searches of peaks in the spectrum of muons in pion and 
kaon decays and of the decays of A'4 produced in neutrino beams and e^e~ collisions. 

As already discussed in the case of mixing with Vg, peak searches provide very robust 
and stringent bounds, by looking at pion decays for masses up to 34 MeV, and at kaon 



decays for higher masses. A detailed review is given in Figs. 1 and 2 in Ref. |92] and 



for masses larger than 100 MeV the limits are reported in Fig. The other limits on 
|V^4p are found in decay searches and are also shown in Fig. ^. They come from beam 
dump experiments |87, |8^, 93, 94, 95] and from direct A'4 production in the detectors 
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Figure 3: Bounds on |Ve4p versus 7714 in the mass range 10 MeV-100 GeV. The areas with solid 
(black) contour labeled tt eiy and double dash dotted (purple) contour labeled K ^ eiy are 
excluded by peak searches jS^L Limits at 90% C.L. from beam-dump experiments are taken 
from Ref. |86| (PS191), Ref.^ (NA3) and Ref. @ (CHARM). The hmits from contours labeled 
DELPHI and L3 are at 95% C.L. and are taken from Refs. |^ and respectively. The excluded 
region with dotted (maroon) contour is derived from a reanalysis of neutrinoless double beta decay 
experimental data p^ . 



DELPHI II , L3 [|g] and CHARM ||]. 
2.2.3 Mixing with Ur 

Heavy neutrinos mixed with r neutrinos can be produced either via CC interactions if a r 
is produced or in NC interactions. The only limits come from searches of A^4 decays and 
are reported in Fig. |. The bounds at 90% C.L. from CHARM Q and NOMAD 



assume production via D and r decays. The DELPHI bound at 95% C.L. ||89|] assumes 
A^4 production in decays and with respect to the bound on |V'e4p and |VJ^4p there is r- 
production kinematical suppression for low masses which weakens the constraint for masses 
in the range ~ 2-3 GeV. 

2.2.4 Electroweak Precision Tests 

The presence of heavy neutral fermions affects processes below their mass threshold due 
to their mixing with standard neutrinos [^] and significant bounds can be set by precision 
electroweak data. The effective //-decay constant G^, measured in muon decays, is modified 
with respect to the SM value and can be related to the fundamental coupling Gp as: 

= GF^il-\VeA\''){l-\V^A\^) . (2.10) 
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Figure 4: Limits on |Vp4p versus 1124 in the mass range 100 MeV-100 GeV come from peak 
searches and from A'4 decays. The area with solid (black) contour labeled K ^iu p% is excluded 
by peak searches. The bounds indicated by contours labeled by PS191 lH], NA3 fs^BEBC |93|, 
FMMF 0, NuTeV H and CHARMII |6| are at 90% C.L., while DELPHI g and L3 @ are 
at 95% C.L. and are deduced from searches of visible products in decays. For the beam dump 
experiments, NA3, PS191, BEBC, FMMF and NuTeV we give an estimate of the upper hmit for 
the excluded values of the mixing angle. 



The fjL — e universality test, done by comparing the decay rate of pions into ev and ^v, can 
be used to constrain the ratio 

for m4 > rriTT |7^. The analysis of experimental data leads to = 1.0012 ±0.0016 

|71], which implies |K=4p < 0.004 at 2a for the least conservative case of |VJj4p = 0. For 



m4 > m,-, the ^ — t universality sets limits on: 

l-\VrA? 



1-IK.-|2^ 



(2.12) 



and can be tested by looking at the r leptonic and hadronic decays which give |V^4p — 
\V^A? = 0.0057 ± 0.0065 ]7l| and |K4p - |V"e4p = 0.0054 ± 0.0064 Q. The most con- 
straining bound on |Ft4P is obtained for |Ve4p, |V^4p = and reads |Vr4p < 0.018 at 2a. 
The unitarity constraint on the first row of the CKM matrix [^] reads 

E IK^™I'= 1 IT/ 12 =0-9992 ±0.0011, (2.13) 

i=i,2,3 ^ 
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Figure 5: Bounds on |yr4p versus 7714 from searches of decays of heavy neutrinos, given in Ref. ||97 
(CHARM) and in Ref. [|| (NOMAD) at 90% C.L., and in Ref. ^ (DELPHI) at 95% C.L. 



and translates into a very strong bound on |V^4p, |VJ^4p < 0.0003 (0.0014), at 1 (2)ct, 
which holds for sterile neutrinos heavier than the A baryon. 

In the presence of heavy singlet neutrinos heavier than half the mass, the invisible 
decay rate of would be reduced with respect to the SM one, r|^j^^, as: 



^-(l-^l^e4p-^I^Ml^-^|K4n. (2.14) 



Z— +inv 



By a standard model fit to LEP data, the effective number of neutrinos is now determined 
to be = 2.984 ± 0.008 pi and provid es a bound on |Vg4p similar to but somewhat 



weaker than the ones obtained by lepton-universality. 

A combined analysis of an old set of unitarity bounds which does not include the 
one from the CKM matrix determination, leads to the following limits at 90% C.L. |V^4p < 
0.012, |y^4|2 < 0.0096 and |K4p < 0.016. If the CKM matrix constraint is included and 
partial cancellations between the contributions of different flavors are taken into account, 



a previous combined study [70| then gives the more robust limits at 90% C.L., |VJ=;4| < 



0.0066, |V/,4p < 0.0060 and |y^4p < 0.018. A very recent analysis @ has updated these 
results using the latest electroweak precision data, except for the CKM observables. They 
find at 90% C.L. 

\VeA^ < 0.003, \V^a\^ < 0.003, \Vr4\^ < 0.006 . (2.15) 

If the constraints from CKM observables are included, we expect the bounds to become 
somewhat stronger, given by |Ve4p < 0.002, |y^4|2 < 4 x IO-MK4P < 0.006 |lOO|. In 
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Table 1: Most stringent model- independent constraints on the mixing elements of the heavy 
neutrino from precision electro-weak measurements. The bounds on |V^4p, t = e,^, r and on 
\VeiV^,4:\ at 90% C.L.. See text for details. 



Mixing element 


Range of m4 


EW Measurement 


\yeAV 


m4 > ©(m^) 


< 0.003 [^] 




m4 > 0{m!y) 


< 0.003 g 




m4 > 0{mr) 


< 0.006 [0] 




10 GeV (100 GeV) [1000 GeV] 


< 0.015 (3.5 X 10-"^) [1.2 X 10'^] 



the following, we take the bound |V;,4|^ < 0.0060 as a conservative reference limit on the 
mixing for comparison with the results of our study. 

Indirect limits on the parameters characterizing heavy sterile neutrinos can be obtained 
from searches for flavour changing neutral current processes such as /i — > 67, — > ee^e~ 
and ^ 



e conversion in nuclei [75|. The branching fraction for ^ 
mixing with heavy singlet neutrinos is given by [^, 74, 75 1: 



67 induced by the 



Br(/i ^ 67) = ^ 



m 



N„ 



m 



W 



(2.16) 



where m' indicates the heavy sterile neutrinos with mass m , and mw is the mass of the 
W boson. The function g{x) is given by 



x{l — 6x -I- 3x^ -I- 2x^ — 6x^ ln(2;)) 
2(1 -x)4 - 



(2.17) 



where g{x) goes from to 1 as x varies from to infinity. 



At present, the branch- 
ing fraction is constrained to be Bi{fi 67) < 1.2 x 10^^^ [101] 90% C.L. imply- 
ing that, for one extra sterile neutrino, |V;4V'^*4l < 0.015 (3.5 x 10"^) [1.2 x 10"^] for 
771-4 = 10 GeV (100 GeV) [1000 GeV]. Similar constraints are imposed by searches for 
the processes — e conversion in nuclei and fi ee^e~ [75|. The current strongest 
bound comes from the search for fi — e conversion in Ti for which the branching 
ratio with respect to the total nuclear muon capture rate is constrained to be 
Br(AiTi ^ eTi) < 4.3 x 10^^ ^t 90% C.L. ||l02|| . For one sterile neutrino, this 
translates into a bound on the following quartic combination of mixing angles 



|V£4|^| < 1.3 X lO^'^ (100 GeV/m4)^, which is weaker than the bounds 
from jjL ^ e"f searches but becomes important at very high values of the masses, 
m4 > 10 TeV. In the presence of more than one sterile neutrino, partial cancellations be- 
tween their contributions are possible, with a consequent weakening of the bounds. Future 
more sensitive searches will further improve these limits. 

The most stringent EW precision constraints are compiled in Table |^ and include the 
bounds reported above from universality tests and lepton flavor changing processes. These 
bounds are obtained barring cancellations between mixing angles and therefore could be 
weakened if some parameters are of the same order. 
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2.2.5 Neutrinoless Double Beta Decay {Ou(3(3) 

The most well studied among AL = 2 processes is neutrinoless double beta decay (Oz^/3/3) 
and the constraints from it deserve special attention. The constraints on |V'e4p for a wide 
range of heavy neutrino masses (10 MeV < m4 < 100 GeV) are shown in Figs. ^ and |3|. 



For heavy neutrinos with mass, fnN^, ^ 1 GeV, the bound is pl| , 34 

y < 5 X 10-^ TeV-^ (2.18) 

m' "1 

The constraint above is very strong and makes it impossible to observe at colliders the 
like-sign dilepton signature with electrons (see Sec. |^. 



3. Lepton-Number Violating Decays 

The key point for the search of lepton-number violating processes in this paper is to consider 
the substantial enhancement via resonant neutrino production. One thus needs to evaluate 
the decay widths of A^4 to various channels. We consider the decay width of the heavy 
Majorana neutrino in two regimes: when the mass is much smaller than that of the W 
boson and when the mass is larger than the mass of the W boson. Based on this, we then 
compute the AL = 2 decay branching fractions for r lepton and K,D,Ds and B mesons. 

3.1 Decay Modes of Heavy Majorana Neutrino 

3.1.1 Decay Modes of Heavy Majorana Neutrino with mass <^ m\Y 

For the LV low energy tau decays and rare meson decays the resonant contribution is from 
a heavy Majorana neutrino with mass of order MeV to GeV. In this section we discuss the 
decay modes of a Majorana neutrino which is lighter than the W boson, so that m4 <C mw- 
The heavy neutrino decays via charged and neutral current interactions to the modes listed 
below. The partial decay widths of the heavy Majorana neutrino with the leading terms 
in mixing and in the massless limit of the final state particles are given below. The full 
detailed expressions for the same are given in Appendix ^ 
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where and are charged (neutral) pseudoscalar and vector mesons, /m are the 

meson decay constants and Vgg' are the CKM matrix elements. 

All the decay modes listed above contribute to the total decay width of the heavy 
Major ana neutrino which is given by: 



i,p i,v e,p e,v 

+ Y 2r^i^2f£2 + ^ r'^^i^a^a + ^r'''^i'^^, (3.8) 



where i,£i,£2 = e,fi,T. For a Majorana neutrino, the AL = process i^P^ as 

well as its charge conjugate AL = 2 process N4 i^P~ are possible and have the same 
width r^^. Hence the factor of 2 associated with the decay width of this mode in Eq. j). 
Similarly, the AL = and its charge conjugate AL = 2 process are possible for the decay 
modes N4 £~V^ and l'^l2i'i^ and hence have a factor of 2 associated with their 

width in Eq. (pl) . 

For the low energy LV tau decays and rare meson decays we consider, the mass of 
the heavy neutrino is in the range 140 MeV < < 5278 MeV. For this mass range we 
list all the possible decay channels for in Table |^ in Appendix |C|. The mass and decay 
constants of pseudoscalar and vector mesons used in the calculation of partial widths given 
in Eqs. (3T) — (3^) are listed in Table in Appendix^ 



Next, we get a rough numerical estimate of the total width of the heavy neutrino with 
the decay modes discussed in Eqs. (|3.lD — ( ^?7|) . To do this, we consider the massless limit of 
the decay products of the heavy neutrino, include only leading terms in mixing of 0(1 V£4p) 
and ignore small factors like vr and calculating the partial decay widths. We 

can only get a rough estimate of the width in this approximation, but it is sufficient to 
see that it warrants the use of narrow width approximation. The two body decays of the 
heavy neutrino have a general form 

p26od, ^ 9Mk!i\Ve4\' ~ ^m^|F.4p ~ (10-^3 \Ve4\') GeV, (3.9) 

where typical values of 771,4 ~ 1 GeV, /m ~ 0.1 GeV and Gp ^ 10~^ GeV~^ have been 
used. The three body decays of the heavy neutrino have a general form 

^3body ^ 9>!^\V,,\^ ^ ^i^|F,4|2 ~ (10-13 |F,4|2) QeV, (3.10) 

lOOvrS ' 1000 I ^ I ^"^1 ^ ' ^ ^ 

where typical values of ~ 1 GeV and ~ 10"^ GeV^^ have been used. The total 
width of the heavy neutrino is then given by 

r^, = (number of decay modes) x {T^'""^y + T^^^'^^) 

50 X (10-13 GeV + 10-13 GeV) |V^4p ~ (IQ-" jV^P) GeV. (3.11) 



As shown above, the width of the heavy neutrino ~ 0{W~^^ l^4p) GeV is much smaller 
than the mass of the heavy neutrino ~ 0(1 GeV) and we can use the narrow width 
approximation to an excellent approximation. 
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Now we look at the lifetime of the heavy Majorana neutrino to determine the decay 
length. The lifetime is given by 

1 1 

TNi ~ 



10-11 IF^P GeV ' 



GeV- 



6.58 X 10 



-14 



V, 



s, 



(3.12) 



which gives a typical decay length ct]^^ ~ 1 x 10^ |V^4|~ m. Note that for a very small 
mixing, |V£4p < ©(lO^^), the may escape from the detector if it is not much heavier 
than a GeV. We will take this effect into account in the following studies. 

3.1.2 Decay Modes of Heavy Majorana Neutrino with mass m4 > niw 

In this section we discuss the decay modes of the Majorana neutrino which is heavier than 
the W gauge boson, so that m4 > mw- The decay modes of the heavy Majorana neutrino 
are to a W 01c a Z gauge boson plus the corresponding SM lepton. The partial decay widths 
for longitudinal and transverse gauge bosons , in static heavy neutrino frame are 

^2 
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= r(iV4 - 
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= r(iV4 - 


UiZi) = 




= r(iV4 - 


-> vhZt) = 



= r(iV4 - 
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\Vu\ (1 - ^lw) 
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^ \Va\^ nii (1 



327r cos|(/ 



\Vii\^ nii (1 - p,zf 



(3.13) 

(3.14) 
(3.15) 

(3.16) 



where /ij are the masses of the gauge bosons scaled by the mass of the heavy neutrino and 
are given by /ij = mf/m\. To obtain the total decay width for A'4, we sum over the charged 
leptons £ and as discussed earlier include the AL = process A^4 
charge conjugate AL = 2 process A^4 

decay width of these modes in the expression for the total width below. 



W^rp as well as the 
i~^W£rp. Hence the factor of 2 associated with the 



N4, 



In Eqs. ( |3.15| ) — (3.17), we have used the relation (see Appendix for details) 



m=l e=e ii=e 

Ignoring terms of order {Vul'^ we have 



since 



UU^ + VV^ = I. 



(3.17) 



(3.18) 



i:i^^4| (3.19) 

m I 

In this approximation, the total width of a heavy Majorana neutrino can be written as 

for 1714 > fnw, 



8ttV2 



oc ^ iVul"^ Gpml{fl^ + ml) for 1714, <C mw, 



(3.20) 
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Figure 6: (a) Top: decay width and (b) bottom: decay length (normalized by \Vu\ ) versus 
mass of heavy Majorana neutrino for real and virtual weak bosons with the inclusion of Higgs 
decay channel for mn = 120 GeV . 



where the expression when m4 <C rriw is obtained from Eq. (3^) and /m are the meson de- 
cay constants. We note that the approximate form of the total width as given in Eq. ( 3.2C| ) 
is only for intuitive purposes to infer the general behaviour of the total width as a function 
of mass. The precise expressions for the total width of the heavy Majorana neutrino as 
given m Eqs. Q, and (Q have been used in the numerical analysis. 



It should be noted that in the SM, if A'4 is heavier than the Higgs boson, then the 
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Figure 7: (a) Left: branching fractions for decay of heavy Majorana neutrino into W* and Z* 
bosons with varying heavy neutrino mass; (b) right: branching fractions for decay of heavy Majorana 
neutrino into longitudinal and transverse gauge bosons in static heavy neutrino frame with the 
inclusion of Higgs decay channel for mn = 120 GeV . 

decay to a Higgs will be present and the partial width is given by 

r'^H ^ p(^^ ^ ^ _l_ |^^^|2 ^3 (1 _ ^^)2^ (3 21) 

In Fig. ^ we plot the decay width of the heavy Majorana neutrino versus its mass normalized 
by the common factor |V^4|^. We can see in Fig. |6|(a) that for a heavy neutrino with 
mass 7714 > "7^^, the decay width increases as Gpm^ as given in Eq. (|330| ). Given the 



rather small mixing parameter, the width remains narrow even for 7774 ~ 0(1 TeV). For a 
lighter neutrino with 7774 <C mw, the width can be very small. The proper decay length is 
presented in Fig. ^(b). We see from this that for 7774 < 20 GeV and |V^4p ^ 10~^ from 
Fig. ^, we have cr ~ 1 /um. 

In Fig. 0(a) we plot the branching fractions of the heavy Majorana neutrino decay to 
Wi and Zh' versus varying heavy neutrino mass 7774. In Fig. ^(b) we plot the branching 
fractions for the decays into longitudinal and transverse gauge bosons in static heavy 
neutrino frame. When the neutrino mass is large, it mainly decays to longitudinal gauge 
bosons and Br(iV4 W+l') ~ Br(7V4 Zv) = Br(iV4 Hv) = 25%. In terms of the 
search at hadron colliders, we prefer to adopt the Wi mode since we wish to reconstruct 
the full event including the lepton number. 

3.2 Lepton-Number Violating Tau Decays 

In this section we examine tau decays into an anti-lepton and two mesons 

T'(pi) ^ ^+(P2) Mf (gi) Nq{q2) (3.22) 
which is a process with AL = —2. The decay amplitude for the above process is given by 



{Pi - qif -1^1 + iTN^rui 



+ (qi ^ Q2), 

(3.23) 
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Table 2: Mass and mixing elements of heavy neutrino and the decay mode constraining them with 
the corresponding experimental bounds on branching fractions. Bounds for AL = 2 tau decays are 
from Ref . [|o| 



Mixing element 


Range of (MeV) 


Decay mode 


^exp 




140 - 1637 
140 - 1637 
494 - 1283 


T~ e"*'7r~7r^ 
T- e+R-R- 


2.7 X 10"'' 

1.8 X lO"'^ 
1.5 X lO"'^ 




245 - 1637 
245 - 1637 
599 - 1283 


/i"'"7r~7r~ 
T- fi+R-R- 


0.7 X lO-'' 
2.2 X lO"'^ 
4.8 X 10-"^ 



where V^^^ and /a./^ are the quark flavor mixing element and the decay constant for 
the meson Mj respectively. From this decay amplitude, we can calculate the transition 
rate F^y and the branching fraction normalised by the tau decay width. In Appendix 
we give the calculations and the full expressions for the decay branching fraction of the 
process (|3.22| ) in terms of the mass of heavy neutrino, m^, and the mixing \V-riVi^\^ . To 
understand the physical picture, we can express the branching fraction in an intuitive form, 
in the massless limit of the final state particles, as 

Br = 5^=F- 



~ 10-3 \V^KMyCKM^2 i^^^^^^i^ (3_24) 

where we have used typical values of 777,4 ~ 1 GeV, fu^ ~ 0.1 GeV, Gp ~ 1 x 10^^ GeV^^ 
and Fjv^ ~ 10^^^ |V£4p GeV. From the simple expression given above one can easily make 
a rough estimate of the required sensitivity and hence the feasibility of observation in terms 
of the mixing parameters for a given model. 

A direct search for LV tau decays has been made at the BaBar detector and the 



limits on the branching fractions were reported in Ref. | 103 |. The experimental limits for 
various decay modes are typically of the order of 10~ , as given in Table |. Fr om the non- 
observation of the LV tau decay modes one can determine bounds on the mixing parameters 
\Vi4Vri^ as a function of the heavy neutrino mass 7774. To do this in a comprehensive 
manner, we carry out a Monte Carlo sampling of the mixing parameters and the mass of 
the heavy neutrino. For simplicity, the mixing elements Ve4,V^4 and Vr4 are allowed to 
vary in the range from to 1. The ranges of mass sampled for the heavy neutrino are 
listed in Table § for the various tau decay modes. We only sample the range of masses that 
lead to a resonant enhancement of the width as the other mass regions have very small 
transition rates as discussed earlier. We then calculate the transition rates and branching 
fractions over the entire range of mixing and mass of the heavy neutrino and the results of 
the Monte Carlo sampling are discussed next. 

The relevant mixing parameters \Ve4Vr4\ and |V^4V^4| are probed as a function of the 
heavy neutrino mass 7774 and are shown in Fig. ^(a) and Fig. ^(b), respectively. Under the 
assumption that the detector was able to reconstruct all the signal events, the region above 
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the curves is excluded by the current direct experimental search for LV tau decays. The 
most stringent bound on |VJ;4Vr4| is of 0(10"^) and comes from t~ e^7r~7r~. The most 
stringent bound on | V^4Vr4| is also of 0{W~^) and comes from fi'^n^ir^ . This is three 

orders of magnitude more sensitive than the limits from precision electroweak data which 

I 1 2 

constrain the square of the mixing IVi^l to be less than few times 10~ . In the absence 
of detection of LV processes the constraints on mixing from peak searches, accelerator 
experiments, reactor experiments and others (collectively called laboratory constraints here 
and henceforth) described in Fig. |2| — Fig. ^ are also applicable here. In the mass region 
probed by LV tau decays the most stringent current constraints are 

|V^4p < 10~^— 10^^ and |Vr4p < 10~^ — 10~^. This would roughly translate into constraints 
on |V^4V^4| < 10~^ — 10~^ and |V^4Vr4| < 10~^ — 10~^ which are comparable to the limits 
from LV tau decay modes. We explore more combinations of mixing elements and also 
provide better constraints on mixing in some mass regions. To summarize, the constraints 
on mixing from LV tau decays are always competitive with or better than precision EW 
constraints and laboratory constraints in the corresponding mass region. The experimental 
bounds can improve in future and an order of magnitude improvement in the experimental 
branching fraction will give approximately an order of magnitude improvement in the 
constraints for the mixing parameters 1 144^^41. More importantly, a detection in one of 
the laboratory experiments implies the existence of a sterile neutrino while a detection in 
one of the modes studied in our analysis would imply LV and hence the existence of a 
Major ana neutrino. 

It should be noted that the intermediate heavy Majorana neutrino is treated as a real 
particle which propagates before decaying. If it exits the experimental apparatus prior to 
decaying, then the signal corresponding to the AL = 2 process cannot be reconstructed 
and no bound could be deduced from the non-observation of such events. In Figs, ^(a) 
and ^(b) , we provide an estimate of the bound on the mixing parameters which takes into 
account the probability of the heavy Majorana neutrino to decay within the detector of 
size Lexp- This probability is given by 

P = l- exp(-Lexpr^) (3.25) 

and for small masses and/or small mixing parameters and consequently long decay lengths, 
it can be approximated with P ~ LexpFTv. We take Lexp = 10 m, the typical size of the 
detectors used in the experiments under consideration. For simplicity, we take A'4 to be 
relativistic but we keep its gamma factor 7 = 1, as a more precise value requires a full 
understanding of the experimental setup. We assume 11441 = = iVril- An estimate of 
the realistic bound on the mixing parameter \Ve4Vr4,\ is then given by 

\Ve4Vri\{= \VX) = ^J\Vei\l,/{L,,p^m), (3.26) 

where |Ve4|^ is the bound obtained assuming that all the A'4 decay in the detector and 
discussed above, and Tjyo is the decay rate for a fully active heavy Majorana neutrino, i.e. 
when the mixing parameter \V£4\ = 1 for i = e,^,r. The bounds remain unchanged for 
large values of the mixing angle and / or large values of 772.4, as the decay length in these cases 
is very short. However, the most sensitive limit on |144T44|(= |^4p) coming from r em: 
searches gets weakened to ~ 4 x 10"^ (4 x 10^^) (1 x 10"^) for = 0.2 (0.5) (1.0) GeV. 
Similarly, the searches for r ^vrvr allows to set a bound on |V^4V^4|(= |V^4p) which 
weakens to \Vf,4VrA\ < 1 x 10"^ (2 x lO^^) (1 x lO^^) for 7714 = 300 (600) (900) MeV. A 
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Figure 8: (a) Left: excluded regions above the curves for {VeiVril versus TO4; (b) right: same as 
(a) but for |V^4T/r4|- The thin black lines correspond to the estimate of the bound (for r — > eirir 
and r /ztttt) once the probability of N4 decay in the detector is taken into account. 



detailed analysis taking into account the experimental setup should be performed in order 
to obtain more precise bounds. 



3.3 Lepton-Number Violating Rare Meson Decays 

We now investigate the LV processes in which a meson decays into two like-sign leptons 
and another meson 

^ i+ipi) t{p2) M2iQ2)- (3.27) 

These decays are similar to the tau decay modes described in the previous section. The 
decay amplitude for the above process is given by 



9^2 jrCKMjrCKM f f T/ T/ ^ 



{QI - PlY - + «rAr^m4 



+ {Pl 



P2) 



{qi - PlY - + ^Tat^ 777-4 



+ {Pl 



P2) 



(3.28) 



(3.29) 



where iA4^ and iAi^ are the decay amplitudes when the meson M2 is a pseudoscalar or 
vector meson respectively and V^^^ and /a^. are the quark flavor mixing element and the 
decay constant for the meson Mj respectively. From this decay amplitude, we can calculate 
the transition rate T^^^ and the branching fraction normalised by the decay width of the 
meson Mi. In Appendix^, we give the calculations for the decay branching fraction of the 
process ( p.27| ) in terms of the mass of heavy neutrino, 7714, and the mixing \Vl-^^iV£^i\. We 
can express the branching fraction in an intuitive form, in the massless limit of the final 



- 21 - 



Table 3: Same as Table |^ but for AL = 2 rare meson decays. The experimental bounds are from 
Ref. @, the bounds for D+ e+e+TT-{K-) are from Ref. |o|. 



Mixing element 


Range of 7714 (MeV) 


Decay 


mode 






140 - 493 


K+ ^ 




6.4 X 10"^^ 




140 - 1868 


D+ ^ 


e + e^vr^ 


3.6 X 10-6 




494 - 1868 


D+ ^ 




4.5 X 10"*^ 




140 - 1967 


D+ 

s 




6.9 X 10""^ 


\Ve4\^ 

\ ' 6'i 1 


494 - 1967 


Dt 

S 


e+e+K- 


6.3 X 10""^ 




140 - 5278 


B+ ^ 




1.6 X 10-'^ 




494 - 5278 


B+ ^ 




1.0 X 10"*^ 




776 - 5278 


B+ ^ 




2.6 X 10"*^ 




892 - 5278 


B+ ^ 


e+e+K*- 


2.8 X 10"*^ 




245 - 388 


K+ 




3.0 X lO"'-* 




245 - 1763 


D+ 




4.8 X 10"^ 




599 - 1763 


D+ ^ 




1.3 X 10"^ 




881 - 1763 


D+ 


o~ 


5.6 X 10""^ 




997 - 1763 


D+ 




8.5 X 10""^ 


1 1 


245 - 1862 


Dt 

S 




2.9 X 10"^ 




599 - 1862 


Dt 

^ S 




1.3 X 10"^ 




997 - 1862 


Dt 


p+p+K*- 


1.4 X 10"^ 




245 - 5173 


B+ ^ 


P^ p^TT^ 


1.4 X 10"*^ 




599 - 5173 


B+ ^ 


p+p+R- 


1.8 X 10"*^ 




881 - 5173 


B+ ^ 


Li^ 0^ 


5.0 X 10^6 




997 - 5173 


B+ ^ 


p+p+K*- 


8.3 X 10^6 




140 - 493 


K+ ^ 


e'^p'^ir^ 


5.5 X 10-^" 




140 - 1868 


D+ ^ 


e'^p~^'ir~ 


5.0 X 10"^ 




494 - 1868 


D+ ^ 


e+p+R- 


1.3 X lO""^ 




140 - 1862 


Dt^ 




7.3 X 10""^ 




494 - 1967 


Dt^ 


e+p+R- 


6.8 X 10"^ 


140 - 5278 


B+ ^ 


e'^p'^ir^ 


1.3 X 10-6 




494 - 5278 


B+ ^ 


e+p+R- 


2.0 X 10^6 




776 - 5278 


B+ ^ 


e^p^p" 


3.3 X 10-6 




892 - 5278 


S+ ^ 


e+p+R*- 


4.4 X 10-6 



state particles, as 



■L Ml 



pMi 

= I 

'LV 

~ (10-16 GeV) TM, |V™V™P |V,,4V,,4|, ' (3.30) 

where we have used typical values of 7714 ~ 1 GeV, /m, ~ 0.1 GeV, ~ 1 x 10^^ GeV-^, 
rAr4 ~ 10~^^ |V£4p GeV and tm^ is in seconds. Using the values for the lifetimes of the 
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mesons in Appendix the branching fractions for the various mesons are given by 



MK) - iF^f^^i^^f \V^,.V,,.l (3.31) 
Br(D, B) ~ 10-^ l^igf^^igf^l' \V,,,V,,,l (3.32) 
Bt{D,) 10-5 |y^f^y^f^|2 \Ve,,Ve,,\. (3.33) 

As mentioned ear her, with the simple expressions above one can easily make a rough 
estimate of the required sensitivity and hence the feasibility of observation in terms of the 
mixing parameters for a given model. 

Searches for rare meson decay modes have been made in numerous experiments. Table 
m summarizes the current experimental limits on branching fractions given by Refs. |99 



|10J]. From the non-observation of these LV rare meson decay modes one can determine 
constraints on mixing parameters |V^j4V^24| ^ ^ function of the heavy neutrino mass m^. 
To do this in a comprehensive manner, we carry out a Monte Carlo sampling of the mixing 
parameters and the mass of the heavy neutrino similar to tau decay. The mixing elements 
Ve4 , Vfj_4 and Vr4 are allowed to range from to 1 for simplicity. Only the range of mass 
that leads to a resonant enhancement of the width is sampled for the heavy neutrino and 
listed in Table |3| for the various meson decay modes. The transition rates and branching 
fractions are then calculated over the entire range of mixing and mass of the heavy neutrino 
and the results of the Monte Carlo sampling are discussed next. 

1 2 I 

For the various decay modes, the mixing parameters probed are 1^441 , IVeAVn^l and 

I 1 2 

|V^4| depending on the final state leptons. Again, we plot the excluded region of the 
mixing parameters as a function of neutrino mass, as shown in Fig. |^ for |V^4|^, Fig. 10 
for |Ve4V^4| and in Fig. ^ for |V"^4|^. The regions above the curves are excluded by the 
current direct experimental searches for LV meson decays. First we plot the limits which 
can be derived assuming that all decay in the detector and give a positive signature. 
The most stringent constraints are from the — > £f£2'^~ mode with mixings of O(10~^) 

I 1 2 I 1 2 

excluded for \VeM , \Ve4V^4\ and \V^4\ ■ This is six orders of magnitude more sensitive 
than the limits from precision electroweak data which constrains the square of the mixing 
|V£4|^ to be less than few times 10~^. Next in sensitivity are the D and Dg decay modes 
with constraints of order few times 10"^ which are similar to the constraints from precision 
electroweak data. The bounds for the same mixing elements are much weaker in the mass 
range above 2 GeV. Even though the limits are weak in this region, it is important not 
to neglect the experimental study of these processes. It only implies that there is a large 
parameter space available for the mass and mixing of heavy neutrinos. 

As discussed for the AL = 2 tau decays, in the absence of detection of LV processes 
the laboratory constraints on mixing described in Fig. ^ — Fig. |5| are also applicable here. 
In the mass region probed by LV meson decays the most stringent laboratory bounds 
are |Ve4p < 10"^ - 10-^ |V;,4p < 10"^ - 10"^ for m4 < 2 GeV and |y^4p < 10"'' for 
m4 > 2 GeV. This would roughly translate into constraints on |T44V^4| < 10^^ — 10^^ 
for 7714 < 2 GeV and |Ve4V'^4| < 10^^ - 10^^ for 7714 > 2 GeV. It should be noted that, 
if these experiments were able to fully reconstruct the signal, the limits from K meson 
decays would be better than the laboratory constraints by at least an order of magnitude 

I 1 2 

in the corresponding mass region. In fact, the constraints on \Ve4\ from the kaon decay 
mode — > e^e^vr" would be more stringent than even the constraints from Qv(3(3 shown 
in Fig. ^. Usually dvj3() experiments have the best sensitivity as they have an advantage 
of a large "effective luminosity" resulting from the large number of nuclei available for 
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decay. The meson (and tau) experiments on the other hand have a small luminosity 
coming from a limited number of mesons (taus) produced in accelerators compared to 
the number of nuclei in Oi//3/3 experiments. It is interesting to note that the resonant 
enhancement in the case of the K meson decay is able to match or improve over the large 
"effective luminosity" of ^v(3(3 experiments. In conclusion, the constraints on mixing from 
LV meson decays are competitive with the precision EW constraints and all the laboratory 
constraints, potentially even Oiy(3f3, in some mass regions. But again, we emphasize that 
the aim of our analysis is to study LV processes and hence Majorana neutrinos. 

We have also taken into account the fact that, for small mixing, only part of the heavy 
sterile neutrinos produced will decay in the detector. We have considered Lexp = 10 m, 
1 1 = |V^4| = \Vr4\ and the gamma factor of N4, 7 = 1, for simplicity. In this case, 
as discussed for the AL = 2 tau decays, the bounds get sensibly weakened. An estimate 
of these bounds is reported in Figs. ^, IC and 11 by thin black lines. We see that the 



bounds get significantly weakened by few orders of magnitude for K eevr, K —>■ e//7r and 
K — > /i^uTT and a careful analysis of these searches should be performed to find the detailed 
bounds on the mixing angles. 

The sensitivity of current direct experimental searches are not adequate to constrain 
mixings for some decay modes. The theoretically allowed branching fraction versus mass 
7714 for such modes is given in Fig. 12. As we can deduce from Table |3| and Fig. 12 all the 



modes are very close to start being probed by direct experimental searches. The experimen- 
tal bounds on branching fractions can improve in future and similar to tau decay modes, an 
order of magnitude improvement in the experimental branching fraction will give approx- 
imately an order of magnitude improvement in the constraints for the mixing parameters 

I 1 2 

\Vi-^4Vi^i\. Currently we do not have any constraints on the mixing parameter \Vr4\ from 
LV rare meson decay modes. Only very weak constraints for BR(i? — > Xt~^t~) < 0{5%) 
exist in a theoretical analysis |105| ]. The similar signature — > r"'"r"'"M" is a possible 
decay mode that would bound \ Vr4\'^ and should be pursued. 

4. Collider Signatures 

In this section we study heavy Majorana neutrinos at hadron colliders. The most distinctive 
channels of the signal involve like-sign di-leptons. It was first proposed in Ref. in the 
context of the left-right symmetric model, and subsequently studied in Ref. |3C, 32, 33| 



We discuss the signatures for a heavy Majorana neutrino and the sensitivity to probe the 
parameters 7774 and V£4 at the Tevatron and the LHC. 

As for the production of a heavy Majorana neutrino at hadron colliders, the represen- 
tative diagrams at the parton level are depicted in Fig. |l^, with the exchange of final state 
leptons implied. The first diagram is via WW fusion with a t-channel heavy neutrino A'4 
exchange, directly analogous to the process of dvj3(3. The second diagram is via s-channel 
A^4 production and subsequent decay. Although in our full calculations, we have coher- 
ently counted for all the contributing diagrams of like-sign dilepton production including 
possible identical particle crossing, it is informative to separately discuss these two classes 
of diagrams due to their characteristically different kinematics. 

The scattering amplitude for the process in Fig. |l^(a) is proportional to Vi^^V^^a and 
the cross section can be expressed as 

a {pp ^ W^W^ ^ eflfx) = (2 - 6e,e,) \Vi,AVe,A\^ ao{WW), (4.1) 
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Figure 9: Excluded regions above the curves for \Ve4\ versus 777.4 from e+e+Mj" searches. 

The thin black line corresponds to an estimate of the bound from K'^ — > e^e^Tr^ once the proba- 
bility of decay of N4 in the detector is taken into account. 
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Figure 10: Same as Fig. |^ but for \VeiV^4\ from e+/i+M2~ searches. 



where c7o(T^W^) is the "bare cross section", independent of the mixing parameters. We 
show the bare cross section at the LHC energy of 14 TeV versus the heavy neutrino mass 
in Fig. 14. This cross section can be at the order of tens of femtobarns. However, due 
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Figure 11: Same as Fig. H but for |V^4|^ from 



li^li^M2 searches. 
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Figure 12: Branching fr-action versus heavy neutrino mass 7714 for decay modes —f t[t^^^2 
not yet constrained by direct experimental searches. The regions below the curve are theoretically 
allowed. 



to the large suppression of the small flavor mixing to the fourth power, the cross section 
is rather small. This process was calculated in Ref. |20| ] under the effective vector boson 
approximation. The authors of Ref. [^0| obtained significantly more optimistic results than 
ours. Further scrutiny indicated that they missed a factor of G^m^/S and their result 
should be scaled down by this factor. The corresponding curve for the Tevatron is not 



shown in Fig. 14 as the bare cross section is smaller by nearly two orders of magnitude. 
Including the small mixing element (to the fourth power) further reduces the cross section 
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Figure 13: (a) Left: Feynman diagram for like-sign dilepton signature via WW fusion in hadronic 
collisions; (b) right: the exchanged coherent diagram which is same as heavy neutrino production 
and decay. 
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Figure 14: The bare cross section (To(W^T^) versus mass of the heavy neutrino irn. 

drastically with no hope of detection at the Tevatron via this mode. 

By far, the dominant production process of heavy Majorana neutrino in hadronic 
collisions is the diagram shown in Fig. |l3|(b). We calculate the exact process, but it turns 
out to be an excellent approximation to parameterize the cross section as 



aipp ^ if If W^) « (2 - 5,,,,) a{pp ^ £fN,)BriN, ^ ifw^) oc j^^'^y'ia - (4.2) 



This observation allows us to study the process in a model-independent way. We can 
rewrite the cross section in a factorized form 

a{pp ^ £f £f ^ if ef j f) = (2 - 6e,e,) Se,,, MN^), (4.3) 

where (To(A'4), called the "bare cross section", is only dependent on the mass of heavy 
neutrino and is independent of all the mixing parameters when the heavy neutrino decay 



-27- 




Figure 15: The bare cross section (To(iV4) versus mass of heavy Majorana neutrino TO4 for the 
Tevatron {pp at 1.96 TeV, sohd curve) and the LHC {pp at 10 and 14 TeV, dotted and dashed 
curves, respectively). 



width is narrow. As seen in Fig. ^, this is indeed the case for 7714 <^ 1 TeV once we fold in 
the constraints | Vf4p < O(10~^) from precision EW measurements. We calculate the exact 
cross section for the dilepton production and use the definition Eq. (14.31) to find the bare 



cross sections (To(A^4), which are shown in Fig. 15 at the Tevatron and the LHC energies 
versus the mass of the heavy Majorana neutrino. Due to the fact that the LHC will start 
its operation at 10 TeV, we have calculated the cross sections at both 10 and 14 TeV cm. 
energy. The production rate is increased at the higher energy by a factor of 1.5, 2.0, 2.5 
for m4 = 100, 550 and 1000 GeV, respectively. We will mainly present our results at 14 
TeV for the rest of the paper. An obvious feature of the cross sections is the transition 
near the W mass. For 7714 < Mw the cross section is nearly a constant due to an on-shell 
W production via the Drell-Yan mechanism with its subsequent leptonic decay to i^N^. 
For m4 > M]y the cross section falls off sharply versus 7714 and the on-shell decay goes like 
7V4 ^£^W^^ £^ hj2. 

The flavor information of the final state leptons is parameterized by 

^ilh - |2 ' (4.4) 

In general the two final state charged leptons can be of any flavor combination, namely, 

e^e^, e^fi^, e^r^, fi^fi^, h^t^ and t^t^. (4.5) 

The constraint from {)uf3[3 as given in Eq. ( |2.18| ) is very strong and makes it difficult to 
observe like-sign di-electrons e^e^. The events with r leptons will be challenging to recon- 
struct experimentally. We will thus concentrate on clean dilepton channels of /i^^u^ and 
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//^e^, although we will comment on our proposal to include the r modes. The correspond- 
ing mixing parameters in our notation will be 

respectively. Given the smallness of |K;4p, we can further simplify our study by exploring 
only two cases: an optimistic case |V^4p ^ |^T4p, |V"e4p and a generic case |V^4p « 
|Fr4p » 1144^, which lead to 

e _ / I^mP (optimistic) „ _ j \Vei\'^ (optimistic) , . 

^'^'^ - 1 i|y^4p (generic) ' ^ ' \ \\V,^\'^ (generic) ' ^^''^ 

4.1 Search for Like-sign Dilepton Signals at the Tevatron 

We now consider the search for A'4 at the Fermilab Tevatron, which is currently running 
at a cm. energy of 1.96 TeV in pp collisions. We concentrate on the clean like-sign jJL^^x^ 
mode 

pp fi^fj.^ iu2 X, (4.8) 

where X is some inclusive hadronic activities common in hadronic collisions. To quantify 
the signal observability, we first impose the basic acceptance cuts on leptons and jets to 
simulate the CDF/DO detector coverage 

p^>5GeV, K|<2.0, p^>10GeV, \rf\ <3.0. (4.9) 

We also smear the lepton momentum by a tracking resolution and the jet energy by hadronic 
calorimeter resolution as 

M = i.5.l0-p^. ^ = ^e3%, (4.10) 

where and Ej are in units of GeV. 

The signal events we are searching for have very unique kinematical features. For the 
purpose of illustration, we choose nn = 20, 50 GeV (below mw threshold) and 100 GeV 
(above mw). First of all, there are two well-isolated like-sign charged leptons. This is shown 
in Fig. |^(a) by a normalized distribution of the minimal isolation Ai^^j = y^Ar/^ + A^^. 
Second, there is essentially no missing transverse energy. However, realistically, the de- 
tectors have finite resolutions as simulated by the Gaussian smearing given in Eq. ( 4.1C| ). 



Consequently, there is always some misbalance in the energy-momentum measurements, 
which is attributed to the missing transverse energies and is plotted in Fig. ^(b). Thirdly, 
due to the existence of an on-shell in the signal process, one would expect to recon- 
struct it by an invariant mass either from the 2(1.2 j system m{Ujj) (in the case of DY 
production) or from the di-jets m(^jj) (in the case of A^4 decay). This is demonstrated in 
Figs, ^(c) and (d), respectively. The above kinematical features motivate us to impose 
the following event selection cuts 

Ai?^^'^ > 0.5, (4.11) 
60 GeV < either m{iijj) or m{jj) < 100 GeV, (4.12) 
< 20 GeV. (4.13) 
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Figure 16: Normalized distributions a~^da/dX for TO4 = 20, 50 and 100 GeV at the Tevatron for 
(a) upper left: the minimal isolation Ai?^*"; (b) upper right: the missing transverse momentum 
^r^] (c) bottom left: the 2i2j system invariant mass m{i£jj); (d) bottom right: the di-jet invariant 
mass m{jj). 



These cuts are highly efficient in selecting the signal events. We illustrate this in Table 
^, in which we calculate the signal rates with the consecutive cuts for = 60 GeV and 
|V^4|^ = = 5 X 10"3 > 1 144 p. Note that the choice of mixing elements used in the 

illustration is motivated by constraints from precision EW measurements. However this is 
for illustration purposes only and in our full analysis we have kept S^^ as a free parameter. 

At the Tevatron energies, the SM contribution to the like-sign dilepton events is rather 
small. The leading background of this type comes from the top-quark production and its 
cascade decay via the chain 

t W+b £+ui b, (4.14) 
t W~b -^W- cue i+. (4.15) 

The background rates and survival probabilities with the consecutive cuts are also given 
in Table ^ We see that the tt background is essentially eliminated by the selective 
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Table 4: The representative signal and background cross sections at the Tcvatron, for ^J, ^J. jj and 
the efficiencies with the consecutive cuts. For illustration, we have used TO4 — 60 GeV, |VJj4|^ = 

|K4|2 = 5 X 10-3 > IK4I'. 





No cut 


Basic cut (4.9) 


+AR (4.11) 


+m{jj),m{Ujj) (H.12D 


(4.13) 


Signal 
a{Po) 
eff. 


319 


108 
33% 


99 
92% 


96 
97% 


96 
100% 


tt Bkg 
a(fb) 
eff. 


78.4 


58.2 

74% 


1.85 
3.2% 


0.04 
2.2% 


0.005 
12.5% 




Figure 17: (To(A^4) with varying heavy neutrino mass 7714 after all the cuts. The two cases corre- 
spond to muon rapidity acceptance at DO and CDF. 



cuts. We have also considered other SM backgrounds coming from the production of 
W^W^jj, W^Zjj. After the selective cuts, all these backgrounds are negligibly small. 



In Fig. 17, we plot the bare cross section aQ{N4) with the basic cuts of Eq. (|4.9|) as well 



as the selection cuts Eqs. ( 4.11| ) — ( 4.13| ). The reduction in rate is mainly due to the basic 



acceptance cuts. For comparison, we have also included two choices of pseudo-rapidity 
cut |r?(/i)| < 2 and \r]{n)\ < 1.5. We now consider the statistical significance of the signal 
observation. In the absence of background events, we use Poisson statistics to determine 
the search sensitivity. We take a signal with 95% Confidence Level (as this is very close 
to 2(7 we call it a 2a effect henceforth) to be 3 events. We can thus translate this to 
the sensitivity to the mixing parameter 

(2 - 6i,i,)aoiN4) Se.e^ L>3, (4.16) 

where L is the integrated luminosity. 

The CDF collaboration at the Tevatron has successfully studied the events with like- 



sign dileptons in a different context [lOt]. Given our event selection, in Fig. |l^ we estimate 
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the sensitivity reach for the mixing parameters versus at the 2a (sohd curves) and 



5a (dashed curves) level at the Tevatron. In Figs. 18 (a— b) (upper-left and upper-right), 
the sensitivity is shown for 5"^^ with 2 and 8 fb^^ integrated luminosity. The horizontal 
dotted lines are the constraint on S^^ ~ I^m4|^ < 6 x 10^^ from an analysis of precision 



EW measurements [^]. The DELPHI ||89| and L3 [|90f| bounds are also given for 



comparison. We find that the Tevatron has the potential to reach the following sensitivity 
for the mass of the heavy neutrino 

J 40 - 130 GeV for 2cr with 2 fb"^; 

\ 10 - 180 GeV (50 - 120 GeV) for 2a (5f7) with 8 fb^^ ^ ^ 

Alternatively, the sensitivity for the mixing parameter can be 

f 2 X 10^5 for 2a with 2 fb"^; , . 

'^'^ \ 5 X 10-6 (2 X 10-5) for 2cj (5a) with 8 fb-^ ^ ^ 



Similar to Figs. 18(a— b), Figs. |T8|(c— d) (lower-left and lower-right) show the results 
for instead. The lower dotted curve in Fig. |l8|(d) is the bound on ~ I^e4p from 
Qvf3[3. We have assumed the same detection efficiencies for n and e. With this assumption, 
the slightly better reach for S'g^ compared to S^^ is due the factor of two difference in 
total rate with identical and nonidentical particles as evident from Eq. ( [4.31 ). With 2 fb"^ 
luminosity, the sensitivity to \Vei\'^ is not close to the stringent bound from the Oi^PP 
decay as seen in Fig. ^. We see from Fig. |l8|(d) that with 8 fb-^ luminosity, the Tevatron 
sensitivity for S'e^ may reach the level of the current bound from Oz^/3/3. From Eq. ( [4.16| ), it 
is straightforward to obtain future sensitivity to mixing parameters (S^^^, Sefj.) by a simple 
scaling of the luminosity. 

4.2 Search for Like-sign Dilepton Signals at the LHC 

At the LHC with a cm. energy of 14 TeV in pp collisions, we adopt the basic acceptance 
cuts on leptons and jets as 

p^>10GeV, |t/^| < 2.5, > 15GeV, |r/^| < 2.5. (4.19) 

The efficiency of these cuts increases with heavy neutrino mass and is 50% for 7714 = 
200 GeV and 80% for 7714 = 800 GeV. The smearing parameters to simulate the AT- 
LAS/CMS detectors are [lO^ 

^ = 36 X 10-5 p^, ^ = e 5%, (4.20) 



Pt Ej 



where pj, and Ej are in units of GeV. 

We again present the characteristic kinematical distributions for the signal. Fig. p^(a) 
shows the normalized distribution of the minimal isolation ARij. The simulated missing 
transverse momentum after the energy- momentum smearing is plotted in Fig. ^(b). The 
invariant masses of the 2i2j system m{iijj) and the di-jets m(^jj) are demonstrated in 



Figs. 19(c) and (d), respectively. We thus design the selection cuts at the LHC as 

Ai2^™ > 0.5, (4.21) 

60 GeV < either 77i(^^ii) or m{jj) < 100 GeV, (4.22) 

< 25 GeV. (4.23) 
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Figure 18: The Tevatron sensitivity to the mixing parameters versus 7714 (a) upper-left: 2(t and Scr 
sensitivity of S^^ with 2 fb^^ integrated luminosity; (b) upper-right: same as (a) but with 8 fb^"'^ 
integrated luminosity; (c) lower-left: 2cr and Scr sensitivity of 5*6;^ with 2 fb^^ integrated luminosity; 
(d) lower-right: same as (c) but with 8 fb^^ integrated luminosity. The horizontal dotted lines in 
(a) and (b) are the constraint on Sjjn ~ |V^4p < 6 x lO^'^ from an analysis of precision EW 



measurements |70 . The DELPHI [89] and L3 |90| bounds are also given here for comparison 



The lower dotted curve in (d) is the bound on Sei_i — \Ve4,\ from QvP(3. 



These cuts are highly efficient in selecting the signal events. We illustrate this in Table |5|, 
in which we calculate the signal rates with the consecutive cuts for m4 = 200 GeV and 

I 1 2 I 1 2 I 1 2 

\V^'i\ = I^T4r = 5 X lO"'^ > \Ve4\ . A gain the choice of mixing elements is motivated 
by constraints from precision EW measurements. However as discussed earlier this is 
for illustration purposes only and in our full analysis we have kept 5^^ and S^e as free 
parameters. 



In Fig. 20, we plot the bare cross section (TQ{Nii) with the basic cuts of Eq. ( 4.19 ) as well 
as the selection cuts Eqs. ( [4.21 ) — ( [4. 23| ) at 14 TeV. The solid (dotted) curves correspond 
to the bare cross section without (with) the Higgs decay channel for mn = 120 GeV. 
The reduction in rate is mainly due to the basic acceptance cuts. We note that the cross 
section with the cuts at 14 TeV is higher than that at 10 TeV by a factor of 1.4—1.6 for 
m4 = 100 — 500 GeV. The sensitivity reach for the mixing parameters to be presented later 



-33- 




1.0 

0.8 
0.6 
0.4 
0.2 
0.0 













/ \ 
/ \ 






i ^ 


I 




\ ■ .' 
\ .' ■ . 


.1.1.1 



200 400 600 800 1000 
m{i£jj) (GeV) 



1.0 
0.8 
0.6 
0.4 

0.2 
0.0 




-m^= 60 GeV 
- m^= lOOGeV 
m,i= 200GeV 
■m,,= 500 GeV 



50 100 150 
m(Jj) (GeV) 



200 



Figure 19: Normalized distributions a-'^da/dX for m4 = 60, 100, 200 and 500 GeV at the LHC 
for (a) upper left: the minimal isolation AR^"^; (b) upper right: the missing transverse momentum 
^rp] (c) bottom left: the 2£2j system invariant mass m{££jj); (d) bottom right: the di-jet invariant 
mass m{jj). 



will be scaled down roughly according to this factor for LHC with cm. energy of 10 TeV. 

As discussed in the previous section, a large SM background comes mainly from top 
quark production and decay via the chain decay t ^ 6 — > c f^. Fortunately, after 
all the selective cuts in Eqs. ( 4.19| ) — ( [4.23| ), the top-quark decay background is essentially 



eliminated and has no remaining events for the expected luminosity of 100 fb~^ at LHC. 

There are several other SM backgrounds coming from like-sign W boson production 
at the LHC energies. First of all, the triple gauge-boson production process 



PP 



W^W^W^ l^l^vv jj, (4.24) 



leads to the irreducible background with two like-sign leptons plus jets. Next, the same 
final state can be produced via the process 



PP W^W^ jj l^e^uu jj, (4.25) 
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Table 5: The representative signal and background cross sections at the LHC, for ^Ji jj and 
the efficiencies with the consecutive cuts. For illustration, we have used 7714 = 200 GeV, |V;^4|^ — 
\Vri\^ = 5 X 10-3 > |Ve4|^ and = 120, 300 GeV 





No cut 


Basic cut 


cut 


+Ai? cut 


+m{jj),m{££jj) cut 








(4.23) 


dill) 


(I4.22D 


Signal 












a(fb) 


0.86 


0.42 


0.37 


0.35 


0.33 


eff. 


- 


48% 


88% 


96% 


94% 


ti Bkg 












a(fb) 


29.6 


16.9 


2.7 


0.075 


0.002 


eff. 




57% 


16% 


2.8% 


2.7% 




ruH =120 GeV 










a 


1.01 


0.42 


0.057 


0.052 


0.050 


eff. 


ruH =300 GeV 


42% 


14% 


91% 


96% 


a(fb) 


1.28 


0.58 


0.066 


0.061 


0.058 


eff. 




45% 


11% 


92% 


95% 




ruH =120 GeV 










a(fb) 


4.2 


1.3 


0.29 


0.17 


0.019 


eff. 


ruH =300 GeV 


31% 


22% 


59% 


11% 


a(fb) 


4.4 


1.4 


0.34 


0.19 


0.025 


eff. 




32% 


24% 


56% 


13% 



where the two jets may come from either QCD scattering or from the gauge-boson fusion 
process. However these backgrounds have two missing neutrinos and can be suppressed 
by a combination of cuts on the missing transverse energy and invariant mass. We also 
analysed the backgrounds coming from Z boson production 

PP jjZZ, pp jjZW. (4.26) 

in which some charged leptons are missing in the detection so that they lead to like-sign 
dilepton events. The backgrounds are very small after the cuts. 

We list the number of background events and efficiency of cuts in Table ^ for a lumi- 
nosity of 100 fb"^ at LHC. The total background is about 7 — 8 events for 100 fb~^ at the 
LHC. The main background is from the W^W^W^ channel and can be further suppressed 
if a tighter missing energy cut could be exploited. For instance, the background events may 
be reduced by half, leaving about 3 — 4 events with j^j. < 15 GeV. 

The last but not least important feature of the signal is the direct reconstruction of 
the resonant mass of ^4 in the final state i^jj- This is shown in Fig. 21 for the SM 
background and the signal with 7714 = 200, 400 GeV. We see the effective reconstruction 
of the resonant mass. For a given mass 7714 in the search, one can further make the event 
selection on m{ijj) 

0.8 m4 < m{ijj) < 1.2 ^4, (4.27) 

to estimate the significance of the signal observation. This loose cut has little effect on 
the signal, but reduces the total background to — 4 events for 100 fb"^ in the range of 
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Figure 20: The bare cross section (Jo(-^4) versus heavy neutrino mass 7714 after all the cuts at the 
LHC (14 TeV). The solid (dotted) line correspond to the exclusion (inclusion) of the Higgs decay 
channel for mn = 120 GeV. 



m4 as shown in Fig. We once again adopt Poisson statistics to determine the search 
sensitivity. The number of signal events needed for 2a significance would be 3 — 11; and 
15 — 44 for 5a significance. In Fig. |2^(a) and Fig. |2^(b), we summarize the sensitivity 
for S^^ and Se^ versus m^, respectively. The solid (dashed) curves correspond to 2a {5a) 
limits on Su/ with the exclusion of the Higgs decay channel. The dotted (dash dotted) 
curves are similar but with the inclusion of the Higgs decay channel for niH = 120 GeV. 
The horizontal dotted line corresponds to constraints on |V^4p < 6 x 10~^ from precision 
EW measurements [^0|. In Fig. p3|(b) the dashed line at the bottom corresponds to the 
limit from Qv(3(5. 

I 1 2 I 1 2 I 1 2 '"i 

In the optimistic case, we assume that 11^41 <C |V^4| and S^^ ~ |V^4| < 6 x 10" . 
The detection sensitivity on heavy neutrino mass can be 

f 375 GeV for 2a; 

1 250 GeV for5<.. ^^'^^^ 
Or alternatively, the mixing parameter can be probed to 

^^^^{llZ-^ for 5a.' (^-^^^ 

In particular, even with the very stringent bound on |V^4p from Oz^/3/3 as indicated by 
the dashed curve in Fig. ^(b), one may still have 2a sensitivity if 7714 ~ mw- 

Our calculations for hadron colliders have been based on parton-level simulations. 



A recent study |33| pointed out that there may be other backgrounds to be considered 
when detector effects are included. One of them is the faked like-sign dileptons from 
the bh cascade decay. The other is due to the QCD multi-jet radiation to degrade the 
reconstruction oiW ^ jj. Those backgrounds can not be easily simulated in particular at 
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Figure 21: Invariant mass distributions of m(£jj) for the signal with 7714 = 200, 400 GeV and 
background processes. 



T 1 r 




(GeV) 

Figure 22: Number of background events vs mass of the heavy neutrino, m4. 

the parton-level. A preliminary analysis including full CMS detector simulations cannot 
support their claim |^4|. Nevertheless, we may consider to design more stringent acceptance 
cuts to further discriminate against the backgrounds. First, common wisdom suggests to 
tighten up the charged lepton isolation requirement 

AR'P/" > 0.8, (4.30) 

which would remove the backgrounds from b, c decays substantially, but a full assessment 
can be made only when real data become available and after the detectors are fully un- 
derstood. Next, we may increase the jet threshold to suppress the initial state QCD jet 
radiation to purify the W jj sample. Our estimate based on a PYTHIA simulation 
shows that the kinematics of a DY-type electroweak process can be largely preserved with 



-37- 




Figure 23: (a) Left: 2a and ha sensitivity for S*^^ versus m4 at the LHC with 100 fb^^ integrated 
luminosity; (b) right: same as (a) but for Se^^ ■ The solid and dashed (dotted and dash dotted) 
curves correspond to limits with the exclusion (inclusion) of the Higgs decay channel for niH = 120 
GeV. The horizontal dotted line corresponds to the constraint on S^^ ~ l^4p < 6 x 10^"^ from 
precision EW measurements pOl. 




100 200 300 400 500 600 700 

(GeV) 



10" 



10"- 



CO 



10 



10"- 



5a 



2a 



L3 



100 200 300 400 

(GeV) 



Figure 24: (a) Left: same as Fig. |g but with the ti ghter cuts of Eq s. (Ug) and ( |4.3l| ); (b) right: 
same as Fig. E3(a) but with the tighter cuts of Eqs. (4.30) and (4.31). 



the appropriate jet threshold. We thus examine the cut 

p>j, > 25 GeV, 



(4.31) 



which results in only about 17% of the events with potential jet contamination. The results 



with the tightened cuts are given in Fig. 24. We see that the stringent cuts severely hurt 



the low mass region, but the effect on the high mass region is modest. 
4.3 Like-sign Dilepton Signals with r in the Final States 

So far we have only presented the results with electron and muon final states and ignored the 
taus. This is due to the experimental challenge of r reconstruction. Given the importance 
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to cover all the lepton flavors, one must strive to include taus on the search list. Besides the 



known experimental practice for r identification at the Tevatron [108|, there are proposals 



to identify r events in connection with the neutrino sector [109]. The central issue is to 
reconstruct the missing momenta from r decays. We can generalize our requirement for 
the charged leptons to the isolated charged tracks presumably from the r decays (e, /i, or 
one-prong and three-prong charged hadrons) 

PT(track) > lOGeV, |7?(track)| < 2.5. (4.32) 

This assures that the parent taus are very energetic. For events with one r and no other 
sources of missing particles, the missing momentum will be along the direction of the 
charged track. We thus have 

p (invisible) = (track), (4.33) 

where the proportionality constant k is determined from the I^rp measurement by assigning 
= /tpT (track). For events with two taus, we generalize it to 

p (invisible) = Kip (tracki) + K,2P (track2). (4.34) 

As long as the two r tracks are not linearly dependent, ki and K2 can be determined again 
from the I^j, measurement. The missing momenta, as well as the r kinematics, are thus 
fully reconstructed. 

Although we believe that the A'4 signals in the modes of e^r^, /x^r^ and would 
be very promising for observation, the background analyses will be considerably more 
involved due to the complication of r reconstruction. Since our simulations are performed 
at the parton level, we are unable to adequately address the background suppression and 
to quantify the signal observability. We thus leave this for future studies. 



5. Summary and Conclusions 

The observation of a LV process would show that neutrino is a Majorana particle unam- 
biguously. Apart from light neutrinos, LV processes involving SM particles can receive 
a contribution from heavy Majorana neutrinos due to mixing. In fact, this contribution 
can be resonantly enhanced for appropriate masses of the heavy neutrino. In the absence 
of observation of LV interactions, the rates for these processes can constrain the mixing 
elements \Vi-^/iVi^ii\ as a function of the mass 7714 of the heavy Majorana neutrino. We 
considered two classes of LV violating processes: (a) low energy AL = 2 tau decays and 
rare meson decays and (b) collider signals for like-sign dilepton production with no missing 
energy implying the existence of Majorana neutrinos. We emphasize the necessity of in- 
volving two charged leptons and no neutrinos in the initial and final states, to be conclusive 
about lepton-number violation. 

For the low energy interactions we evaluated the transition rates and branching frac- 
tions as a function of the mass and mixing of the heavy neutrinos. We then translated the 
current experimental bounds from direct searches into limits on {Vi^^^Vi^^l as a function of 
the mass 7714 of the heavy neutrino. Amongst the rare meson decays, the ^^^^vr" 
decay mode currently gives the most sensitive experimental limits on \Ve-^4Ve^4\. Poten- 
tially, these constraints are six orders of magnitude more stringent than the constraints 
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from precision electroweak data which hmit |V£4p to few times 10~^. As the intermediate 
heavy sterile neutrino is a real particle which might exit the detector if the decay length is 
longer than the detector size, for very small mixing angles the bounds get weakened but 
arc still much more stringent than the electroweak precision constraints. This effect should 
be taken into account and a detailed analysis of past experimental data is required in order 
to find the precise limits on the mixing angles. Next in sensitivity are the D and Ds meson 
decay modes with constraints of the order of 10~^. Again, these are competitive with if 
not better than constraints from EW precision data. The other processes (in other mass 
ranges) have very weak experimental limits, weaker than EW precision data and essen- 
tially do not impose any meaningful bounds on {Vi-^^Vi^^l. This implies that more accurate 
experimental studies on those rare decays should be strongly encouraged. In particular, 
many interesting processes of D, B decays have not even been experimentally probed as 
well as those with a r lepton in the final state. Among the r-decay modes the best limits 
come from t~ £^-k~-k~ . The other r-decay modes have sensitivity of order 10~^ to 10~^. 
Again, the constraints from r decay modes are competitive with or better than constraints 
from precision EW data by 2 to 3 orders of magnitude. The experimental bound on LV 
processes is expected to improve in the future. The future sensitivity of the square of the 
mixing parameter will increase approximately by an order of magnitude for every order of 
magnitude improvement in experimental bounds on branching fractions. We have shown 
that the low energy AL = 2 r decays and rare meson decays can be very strong probes to 
discover or constrain the mass and mixing of heavy Majorana neutrinos. Even those decay 
modes which do not impose strong constraints should not be neglected. It only implies 
that a large range of the parameter space is available for exploration. 

In addition to analyzing the LV tau and meson decay modes and precision EW mea- 
surements we also compiled the constraints on the mixing elements (| 144^7 l^4p l^4p) 
from peak searches, accelerator experiments, reactor experiments and others - collectively 
called laboratory constraints. In the absence of detection of LV, the laboratory constraints 
and the ones from AL = 2 processes can be compared. The constraints on mixing from 
LV tau decays are always competitive with or better than laboratory constraints in the 
corresponding mass region while the constraints from LV meson decays are competitive 
with laboratory constraints only in some mass regions. We note that we explore more 
combinations of mixing elements and also provide better constraints on mixing in some 
mass regions. More importantly, a detection in one of the experiments analyzed to obtain 
laboratory constraints implies the existence of a sterile neutrino while a detection in one 
of the AL = 2 tau or meson decay modes studied in our analysis would imply LV and 
hence the existence of a Majorana neutrino. We pointed out the fact, often ignored in the 
literature when analyzing low-energy processes, that a heavy neutrino might decay outside 
the detector if it becomes long-lived for low mass (less than a GeV) and/or very small 
mixing. 

For the collider signals of heavy Majorana neutrinos we looked for the definitive lepton- 

number violating like-sign dilepton production and no missing energy. Such signals have 
low backgrounds and have the potential for discovery of heavy Majorana neutrinos. At the 
Tevatron, with the current and future integrated luminosities, we find that the mass of the 
heavy Majorana neutrinos can be probed up to 



m4 ~ 



{ 



10 - 130 GeV (10 - 75 GeV) for 2a {5a) with 2 fb"^ 
10 - 180 GeV (10 - 120 GeV) for 2a (5a) with 8 fb'^ 



(5.1) 
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Alternatively, the sensitivity for the mixing parameter can be 



2 X 10-5 (10-^) 
5 X 10"^ (2 X 10-5] 



for 2a {5a) with 2 fb'^ 
for 2a {5a) with 8 fb-\ 



(5.2) 



This will surpass the DELPHI ||89| and L3 ||90|, [Tig] 95% C.L. bounds. 

The sensitivity for heavy Majorana neutrinos can be extended significantly at the LHC. 
With 100 fb"^ of integrated luminosity, 



m4 ~ 375 (250) GeV for 2a {5a); 
or alternatively, the mixing parameter can be probed to 



5, 



w, 



Ml 



7 X 10"'^ (3 X lO-*") for 2a {5a). 



(5.3) 



(5.4) 



The sensitivity at LHC will go well beyond the DELPHI and L3 95% C.L. bounds 

in both mass reach and mixing, and beyond the current bound on |144p from Oi^/3/3. In 
summary, there is a rich avenue of possibilities for discovering or constraining the elusive 
Majorana neutrinos. 
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A. Lepton mixing formalism 

In this appendix, we illustrate the parameterization for the lepton sector, although we 
have not followed the relations literally, assuming that new physics beyond this minimal 
formalism exists. The leptonic content in the theory includes three generations of left- 
handed SM SU{2)i doublets and n right-handed SM singlets: 

i^aL=(^^) , N,R, (A.l) 

where a = 1, 2, 3 and 6 = l,2,3,---,n(n>2 for at least two massive neutrinos). 
The leptonically universal gauge interactions involving neutrinos are of the form 

-C= i^W+Y^v;;^ riaL + \i.c\ + -^Z^Y.^ ^""^-L- (A.2) 
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The gauge-invariant Yukawa interactions are 

(3 3 n 

^ fi, HkR + Y,Y1 fab L7l HN,r I + h.c. (A.3) 
a, 6=1 a=l b=l 

where H is the SM Higgs doublet and H = iT2H*. After the Higgs field develops a vev 
{H) — > \/2 , the Yukawa interactions lead to Dirac masses for the leptons 

(3 3 n \ 

Y,^^abhR + Yjl2^ <b NbR + h.c. (A.4) 
a,b=l a=l b=l J 

where the mass matrices are given by the vev times the corresponding Yukawa couplings 

= fabV/V2. 

The 3x3 mass matrix m' can be diagonalized by two unitary rotations among the 
gauge interaction eigenstates II, Ir 

Ol m} Or = diag(me, m^,m^), la = Oat ^, (A.5) 

where I = e, /i, r are the mass eigenstates, which define the charged lepton flavors. The 
Dirac masses as well as interactions with the Higgs boson for the charged leptons now have 
the standard form 

^ ¥f - 

-Cir = V mWl + — ) £ ^. (A.6) 
^-^ V 

If the Yukawa interactions of Eq. ( [A.3| ) are the whole source for neutrino mass, then we 
would have min(n, 3) massive Dirac neutrinos. 

To complete the neutrino mass sector, there is also a possible heavy Majorana mass 
term 

= o E Bbb' N,,R + h.c. (A.7) 



2 

b,b'=l 



where a charge conjugate state is defined as tp'^ = Cip'^ {if)^ = ijFC)^ and a chiral state 
satisfies {iP^)t = (^-r)^) with r = — , + for The full neutrino mass terms thus read 

^ I 'A n n \ 

-C=2 Y.T.('^<bNbR + N^mlu'aR)+Yl N^B,,, N,,R\+h.c. 

y a=l 6=1 b,b'=l J 



where we have used the identity UaL mab Njjr = N^^ rriba v^r-. 

The mass matrix can be diagonalized by one unitary transformation 



where the mass eigenvalues are of the order 



^^a,^"^. ML,^B. (A.IO) 
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L is a (3 + n) x (3 + n) unitary matrix and can be parameterized as 

\ ^nx3 J^nxn J 

The relation between the gauge interaction eigenstates and the mass eigenstates are given 

by 

'""^-ILf Z) , (A.12) 



with the mass eigenstates Um {m = 1,2,3), N^' {m' = 4, •••,3 + n). The diagonahzed 
(Majorana) mass terms of Eq. ( |A.8 ) thus read 



3+n 



-C = M E ^ + E ^^'L ^rn'R + h.C. , (A.13) 

\r?i=l m'=4 / 

with the mixing relations between the gauge and mass eigenstates 

3 3+ra 3 3+n 

VaL = Uaml^mL + Yl ^am'N^'L^ N^L = Y ^bmV^L + ^ ^bm'N^'L^ (A.14) 

m=l m'=4 m=l m'=4 

3 3+n 3 3+n 



Note that the unitarity condition for L leads to the relations 

UU^ + VV^ = U^U + X^X = h^^, (A.16) 
XX^ + YY^ = V^V + Y^Y = Inxn- (A.17) 

Parametrically, UU^ and Y^Y ~ VV'^ and X'^ X ~ ©(m^/MAr). 

In terms of the mass eigenstates, the gauge interaction lagrangian Eq. ( |A.2| ) can be 
written as 

/ r 3 r 3+n \ 

-C = -j=W; Yl E (U^OlU ^I^PlI + E E i^^OL)m'i N^J^Pl^ + h.C. 

\^=em=l ^=eni'=4 / 

/ 3 3+n 
+ IT^Z^ E iU^U)m,m, lA^l^PLiym, + Y (^^^)m[n.',N^r PlK 



^mi,m2=l m[,m'2=4: 
3 3+n 



2cosvy 



E E {U^VU^^m'^T^l^PLN^'^+i^-c. . (A.18) 



^mi=l m2=4 



To make the couplings more intuitive, we define the combination matrices by 

U^'^ = Olu, V^^ = Olv, U''^ = UW, U'"' = U^U, V^^ = VW. (A.19) 
We thus rewrite the gauge interaction lagrangian by one mixing matrix for each term 

(T 3 T 3+n \ 

E E t^i'Pl^ + E E ^^"PLi + h.C. 
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Figure 25: Feynman rules for the charged current vertices in terms of the neutrino mass eigenstates, 
as given in Eq. (A. 22). 



3+n 



2 cosw 



^mi,m2=l mj,m2=4 
3 3+n 



2 COSVF 



These coupHngs along with the mixing matrices Eq. ( A.19| ) give the most general leptonic 
interactions of the charged and neutral currents in terms of the mass eigenstates. Alter- 
natively, the neutral current interactions can be aligned along with that of the charged 
currents when rotating left-handed neutrinos in the same way as the charged leptons, 



3+n 



l^aL = iOL)ai l^iL, OT V^l = ^iOiU)e^UmL + E (^L^)^™' A^^'L • (^-21) 



m=l m'=4 



It may be convenient in certain practical calculations to rewrite the neutral current inter- 
actions in terms of their flavor eigenstates 

(T 3 r 3+n \ 

E E T^I^PlI + E E K^>1^PlI + h.C. 

l=e m=l l=e m'=A J 



+ V7£^^>^ ( E E ^^rn ^^Pl + E E y^rn' N^'l^ Pl l^i ] + h.C. + .(.A.22) 



r 3 T 3+n 



2 COS 9w 



-e m=l £=e m'=4 



where we have dropped the superscripts for U, V defined in Eq. ( |A.1S| ), for simplicity as 
adopted throughout the text. 

For the reader's convenience, we give most of the corresponding Feynman rules for the 
interaction vertices, listed in Fig. 25 for the charged currents, and in Fig. ^ for the neutral 
currents. The Feynman rules for the other diagrams can be easily deduced from the ones 
that are explicitly written down in Fig. |2^ and Fig. 

Finally, the heavy neutrino interactions with the Higgs boson read 

T 3+n 



= - E E ^fn'PL + h.C. + ... (A.23) 

l=e m'=4 

The corresponding Feynman rule for the interaction vertex is given in Fig. p^. 
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Figure 26: Feynman rules for the neutral current vertices in terms of the neutrino mass eigenstates, 
as given in Eqs. (A. 20). 




Figure 27: Feynman rule for the Higgs vertex in terms of the heavy neutrino mass eigenstates, as 
given in Eq. (A. 23). 

B. General amplitude of AL = 2 processes 

The charged current interaction lagrangian in terms of neutrino mass eigenstates is 

T 3 T 3+n 

-Ccc = -4^;(E E f^t* ^l^'PLi^Y. E ^l^' ^^'PLi) +11.C. (B.l) 

e=em=l e=em'=4 

where Pl = ^(1 — 75). The leptonic AL = 2 subprocess W~W~ — > £^£2 is induced by the 
product of two charged currents 

3 3+n 

Mf:, ^ E ^hmUt^m {hl^PLym){T2YPLiym)+ E Km'VlZ^' {hl>'PLNm'){T2Y PLNm'), 
m=l m'=4 

(B.2) 
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which can be rewritten using charge conjugation as 

3 3+n 
m=l m'=4 

(B.3) 

The Majorana neutrino fields can be contracted to form a neutrino propagator, and the 
transition matrix element is thus given by 



+ y E (^7^^l)^ ^^^^^^^^ in^Pn^W (B.4) 

m'=4 ^^m' ™ 

where (7 is the momentum exchange carried by the neutrino. The g' term vanishes due to 
the chirality flip. Including the crossed diagram {£1 <-> £2) the leptonic amplitude then 
becomes 

^Up = T E ^llmKm ^u^Wl-, + ) PrV2 

2 3+n 

+ 17 2^ l/^im' l^^am' "^A^w X 



2 ^ 

?Tt'=4 



^{-2 2 ' T + ^2 2 ' T (B.5) 

For the light Majorana neutrinos, namely, m = 1,2,3 the masses ~ 0(eV) 

and for the heavy Majorana neutrinos , the masses ^n^, ~ ©(MeV — GeV) for the low 
energy processes we consider. The heavy Majorana neutrino contribution has a resonant 
enhancement when q'^,q''^ ~ and is the dominant one. The light Majorana neutrino 

contribution however encounters a severe suppression due to the small neutrino mass like 
JM^. Hence we can neglect the contributions of the light Majorana neutrinos and 

the Y^^m=i P^^^ amplitude drops out. 

In principle all the heavy Majorana neutrinos will contribute to the amplitude but in 
our analysis we only consider the contribution of one of the heavy neutrinos, in particular 
the lightest one for simplicity. So the amplitude can now be written as 



H^p = T^^i4^^.4 m,-A- + '^^^ |P^^2. (B.6) 

' 2 \ - m| + ^^Ar4 777-4 q'^ - ml + iVN^mi 



We can rewrite the above amplitude as 

■^lep = irVe.AVe^A nii— + —Vi^iVt^i rm— 

2 q'^ - ml + 11^^1714, 2 q''^ - ml + iTj^^mA 

= Mi+M2. (B.7) 

When ~ 7774, Aii has a resonant contribution and when q'"^ ~ m\, M2 has a resonant 
contribution. In general, q ^ q' , and it is convenient to split up the individual resonant 
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contributions by the Single-Diagram-Enhanced multi-channel integration method | 112 |. To 
do this, define the functions 

IMP 



Then the amplitude squared is given by 

2 



(Bi 



(B.9) 



The amplitude squared splits up into the functions fi defined above and the phase space 
integration can be done for each fi separately. This helps to make convenient simplifications 
for the phase space integration and the computation can be carried out in parallel. The 
contributions from each fi can be added up after phase space integration. 



C. Decay modes of heavy Majorana neutrino 

In this section we will discuss in detail the decay modes of the heavy Majorana neutrino 
A^4, with mass much smaller than the mass of the W boson, mw. From EW precision 
measurements the mixing elements |V£4p < O(10~^) and the higher order terms in mixing 
would be very small and can be ignored. Hence the widths are presented only up to leading 
terms in mixing. The charged current and neutral current vertices of N4 with the mixing 



elements are given in Fig. ^ and Fig. ^ With increasing mass of the heavy neutrino new 
decay channels open up and can be classified into two body and three body decays. The 
decay width scales as the third and the fifth power of the mass(m4) for two and three body 
decays respectively. 

1) A'4 — > i~P~^ where £ = e, fi,T and is a charged pseudoscalar meson. This decay 



mode has charged current interactions only as shown in Fig. 25 and the decay width is 
given by 



r^p = r(iV4 ^ rp+) = -^fj> \Vgg,\^ \Vm\^ ml ii{fie,i^p) 



= [{1 + X - y){l + x) - ix]X^l,x,y), 
A(a, b, c) = a"^ + b'^ + - 2ab - 2bc - 2ca, (C.l) 



where fp is the meson decay constant and Vqq' are the CKM matrix elements, fii and fip 
are the masses scaled by the mass of the heavy neutrino and are given by = m1/m\. 

2) A'4 — > v^P^ where = Ve-,T^fi:T^T and P^ is a neutral pseudoscalar meson. This 
decay mode has neutral current interactions only as shown in Fig. ^ and the decay width 
is given by 

T-^p ^ r(iV4 ^ = ^fl |y,4p ml (1 - MP)^ (C.2) 

b47r 

where fp is the meson decay constant, fip \s the mass of the neutral meson scaled by the 
mass of the heavy neutrino and is given by /ip = m?p/m\. The mass of the light neutrino 
~ ©(eV) [^] is much smaller than the mass of A'4 ~ ©(MeV — GeV) and can be neglected 
to a very good approximation. We have set the mass of the light neutrino to zero in the 
expression for the width above and henceforth. 
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3) iV4 



'V~^ where i = e,iJ,,T and is a charged vector meson. This decay mode 



has charged current interactions only as shown in Fig. 25 and the decay width is given by 



r 



[(l + x-y)(l + x + 2y) - 4x]A^ (1, x, y), 



|V£4p m\ l2{fJ.e,iJ,v), 



(C.3) 



where /y is the vector meson decay constant and Vqqi are the CKM matrix elements, ni 
and nv are the masses of the lepton and the vector meson scaled by the mass of the heavy 
neutrino and are given by fii = mf/m\. 

4) ViV^ where = Ve^v^, and is a neutral vector meson. This decay mode 

has neutral current interactions only as shown in Fig. 26 and the decay width is given by 



h{x,y) = (l + 2y)(l-y)At(l,x,y), 



q2 



(C.4) 



where fv is the meson decay constant, /xy is the mass of the neutral meson scaled by 
the mass of the heavy neutrino and is given by /iy = rny/rn^. Ky is the vector coupling 
associated with the meson and is expressed in terms of = sin^ 9^, where 0^ is the 
Weinberg angle. The values of k for the various vector mesons are: k = hx^, for and lv; 



(4 + 

5) iV4 



2 ^1 



for K*^,K"^ and 6: and k 



fi - "^x 



for 0*^,0^" and J/V'. 



-=:*0 



^2 ^^^2 where ^1,^2 = e, ;U,t with li 7^ ^2- This decay mode has charged 



current interactions only as shown in Fig. 25 and the decay width is given by 



(1-^)2 



h{x,y,z) 



12 



ds , 



y^){l + z^-s)\-2{s,x\y^)\-2{l,s,z^), (C.5) 



(a;+j/)2 



where /i (0,0,0) = 1, Xj are the masses scaled by the mass of the heavy neutrino and are 
given by Xi = rrii/m/i. The mass of the light neutrino ~ 0{qY) is much smaller than the 
mass of A'^4 ~ ©(MeV — GeV) and hence can be neglected compared to the mass of A'4. 
We have set the mass of the light neutrino to zero with very good approximation in the 
expression for the width above and henceforth. 

6) A^4 — > vi^l'^'^t where ^1,^2 = e,//, r. Both charged current and neutral current 
interactions as shown in Fig. 25 and Fig. 26 are relevant for this mode and the decay width 
is given by 



96^ 



9l9r + Seie29R]h{xue^,xe2,xe. 



+ 



iaif + {dnf + ^iii2 (1 + 25!)) h (xu,^ , xe^ , xg^ 



/ds 1 1 

^(1 - s)X^{s,y\z^)X^{l,s,x^), 



(C.6) 
(C.7) 
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where /2(0, 0,0) = 1, Ii{x,y,z) has been defined in Eq. ( C.5 ), Xi are the masses scaled by 
the mass of the heavy neutrino and are given by Xj = mi/m^, = — \ + Xw, g\i = x^, and 
Xw = sin^ Oyj = 0.231, where Oyj is the Weinberg angle. 

7) A^4 Vi^vv where V(^^ = Ve^v^^^r- This decay mode has neutral current interactions 



only as shown in Fig. 26. Using the massless approximation for the neutrinos as described 



above the decay width has a simple form given by 

r-.,- = ^ TiN, ^ v,,v,,vi;) = ^m,,? ml (C.8) 

All the decay modes listed above contribute to the total decay width of the heavy 
Major ana neutrino which is given by: 



e,p ly e,p e,v 



where i,£i,i2 = e,ii,T. For a Majorana neutrino, the AL = process A^4 — > £~P^ as 
well as its charge conjugate |AL| = 2 process N4 i~^P~ are possible and have the same 
width, r^^. Hence the factor of 2 associated with the decay width of this mode in Eq. (C.9). 



Similarly, the AL = and its charge conjugate |AL| = 2 process are possible for the decay 
modes iV4 —>■ £~V'^ and — > £^£2^1^ and hence have a factor of 2 associated with their 
width in Eq. (|a9|). 

As mentioned earlier, new channels open with increasing mass of the heavy neutrino. 
For the low energy LV tau decays and rare meson decays we consider, the mass of the 
heavy neutrino is in the range 140 MeV ^ m4 5278 MeV. For this mass range we 
list all the possible decay channels for A'^4 in Table The mass and decay constants of 



pseudoscalar and vector mesons used in the calculation of partial widths given in Eqs. ( |C.l 
are listed in Table |^ in Appendix ^ 



D. Lepton- number violating tau decay 

The decay amplitude for lepton number violating tau decays can be separated into leptonic 
and hadronic parts, 

iM = {Mlep)^,u{Mhadr. (D.l) 

For the tree level amplitude, the hadronic part can be expressed in terms of the decay 
constants of the mesons in a model independent way. The box diagram includes hadronic 
matrix elements which cannot be simplified in terms of decay constants and needs to be 
evaluated in a model dependent way. We expect the tree level amplitude to dominate and 
do not include the box diagram. It has been argued that in certain cases for rare meson 



decays sub- leading contributions may be appreciable |113| , 20 1. Even in such a scenario 
the difference will not be important at the current level of sensitivities and we include the 
more conservative limit from tree level diagrams only. The tau decays and the rare meson 
decays are crossed versions of each other and the above arguments are true for both. 
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Table 6: Decay modes of heavy Majorana neutrino based on its mass 1714. 



Mass of heavy 
neutrino (MeV) 


Decay mode of 
heavy neutrino 


Mass of heavy 
neutrino (MeV) 


Decay mode of 
heavy neutrino 






> nip + nir = 1880 


A^4 fl^T'^Ur + C.C 
A''4 T^H^l/p + C.C 


> 2me = 1.02 


TV 7- -±- 

iV4 V(^e. 


> m^ + m^ = 1920 


A'4 T~-K~^ + C.C 


> nie + = 106 


A'4 fi~e~^Ue + c.c 


> me + niDs = 1970 


A'4 ^ 6^1)+ + C.C 


> m^o = 135 




> mp+ mo = 1980 


A'4 ^ + C.C 


> rUe + m-jr = 140 


N4 — > e~7r"'" + c.c 


> m£)*o = 2010 


A'4 ^ ^'£-0*'^ 


> 2m^ = 211 




> 771-^*0 = 2010 


— *u 
Af4 ^ zv^Ll 


> + m,r = 245 


A''4 ^ fl^TT^ + C.C 


> me + niD* = 2010 


N4 e-D*'^ + c.c 


> rrie + rriK = 494 


A^4 ^ e~K+ + c.c 


> mp +mDs = 2070 


A''4 fi'Df + c.c 


> mrj = 548 


A'4 ^ t'f?? 


> me + m£)* = 2110 


N4 e-D*+ + c.c 


> + rriK = 599 


A^4 fi-K+ + c.c 


> mp+ mo* = 2120 


A^4 ^ ^i-D*+ + c.c 


> nipO = 776 


N4 lygp" 


> mp+moi = 2220 


Af4 ^ fi'D*+ + c.c 


> nie + mp = 776 


JV4 — > e /? + c.c 


> m^ + m/^ = 2270 


A''4 ^ r"K+ + c.c 


> m^; = 783 


A'4 — > z/^u; 


> m-r- + mp = 2550 


A'4 — > T~ p'^ + C.C 


> nip + nip = 882 


A'4 ^ + c.c 


>mr + m*^ = 2670 


A^4 ^ r-K*+ + c.c 


> me + rriK* = 892 


iV4 — > e A ^ + c.c 


> m^^ = 298U 




> 171^*0 = 896 


A^4 ^ i^iK*^ 


> mj/^ = 3100 


N4 viJli^ 


> m^,o = 896 


N4 ^ i/^Z*" 


> 2m^ = 3550 


N4 V^^T^T^ 


> m^' = 958 


A'4 i^ev' 


> m-r- + mo = 3650 


N4 T-D+ + C.C 


> nip + mx* = 997 


N4 ^-K*+ + c.c 


> m-r- + m£)^ = 3750 


A^4 T-D+ + c.c 


> = 1019 




> mr + mo* = 3790 


Af4 ^ r-Z)*+ + c.c 


> nie + "Ir = 1780 


A'4 e^r+z^T- + c.c 
A^4 T~e+fe + c.c 


> mr + mn* = 3890 


Af4 ^ r-D*+ + c.c 


> me + mo = 1870 


A^4 e-D+ + c.c 







The leptonic part of the subprocess r ^ *VF * is obtained by crossing the 



ampUtude in (B.6) 



m4 



-ml + irNim4 

Combining the hadronic and leptonic parts, the decay amphtude for 
is given by 



iM = {Miep)^,uMl,^M'it^ + (Ml 



Mo 



(D.2) 

(D.3) 
(D.4) 



(Pi - gi)^ - m| + irAr4m4 
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M1+M2, 



(D.5) 
(D.6) 



where V^^^ are the quark flavor-mixing matrix elements for the mesons and /m^ are 
meson decay constants. Then the functions, fi and /2 defined in Eq. ( [B.8| ) are given by 



/i 

/2 



FrA 

4 + 1^ 



ial + b^)A + {aia2 + b^)C 



+ (91 ^ 92) 



(4 + + (af + b^)B 

fi{qi ^ 92), 

MPi^Qj) = ^{Pi ■ Qi){P2 ■ q2){qi ■ 92) - 4m|,j^(pi • q2){p2 ■ q2 

- ^'rnlhiPi ■qi){P2-qi) + 2m|f^m|^2(pi -pa), 
B{pi,qj) = A{qi ^ ^2), 
C{pi,qj) = 4(pi •P2)(gi -92)^ -A{pi,qj), 
C{qi ^ 92), 



D{pi,qj) 

Fr 



MiJM2\ 



CKM^rCKM\2\ 
V I 



M2 



VriVul'^ml 



(Pi - qi,2f - ml; b = Tn^tji^. 
The decay width for the LV tau decay is then given by 



fidPSsi + / f2dPS: 



IT 



dPS32 = dPSsiiqi ^ q2), 



'32 



''^^''^-dyidy2dyzdyi 



m 



cl 



(D.7) 
(D.8) 

(D.9) 
(D.IO) 
(D.ll) 
(D.12) 
(D.13) 
(D.14) 



(D.15) 

(D.16) 
(D.17) 



where dPS^i and dPS'32 are the phase space factors obtained by conveniently clustering 
two different sets of particles to enable applying the narrow-width approximation easily. 
yi to 7/4 are rescaled angular variables with integration limits < < 1. As seen in 
Sec. p.l.l| , the width of the heavy neutrino is very small compared to the mass and hence 
we can apply the narrow-width approximation. 



/ 



(m; 



m 



riV4^ 



S{m1_ — m\)dw?^. 



(D.18) 



Applying the narrow-width approximation as described above and integrating over the 
(5— function we get 



fidPSsi 



FrA-K^ 

,2 2 2 



{4 + b'^)A + b'^{B + C + D) 
{4 + b'^)A + b'^B 



I 



X2 {m^, mij^ , mi)\'2{ml, mj,mlj^)j dyidy2dyzdyA, 
f2dPS32 = J fidPSsiiqi ^ q2) 



(D.19) 
(D.20) 
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Now, we can find the decay rate from Eq. ( D.15 ). Normalized to the r decay width 
r,- = G'jpm^/192iT'^, the corresponding branching fraction is Br = TJ^/Tt-. The masses and 
decay constants of mesons are hsted in Table 0. The CKM matrix elements and r mass 
are taken from the Particle Data Group (PDG) [p9[: 



1777 MeV, \Vud\ = 0.9738, |Ks| = 0.2200. 



E. Rare meson decay 

The rare meson decays 

have the same Feynman diagrams as tau decay. The meson M2 can be a pseudoscalar or 
vector meson. The decay amplitude when M2 is a pseudoscalar meson is given by 



(E.l) 



Next we consider the case where M2 is a vector meson. The decay amplitude is given by 



= MX +mX. 



{qi -piY - 77142 + iTAT^ 777-4 
{qi-V2Y - niA? + 7rAr4m4 



(E.2) 



Similar to tau decay, we define the functions /f , f2 and fY , fY for the pseudoscalar 
and vector mesons respectively as given in Eq. ( |B.8| ). ff and fY turn out to have the same 
form and are given below. 



fPiV) 

Jl 




[ a? + 62 


fP{V) 


= fi [Pi ^ 




= 8{pi ■ qi){p2 



{4 + b^)AP(^) + iaia2 + b^)C 



P{V) 



{4 + 62)A^(V) + (a2 + 62)5P(y) 
{qi ■ q2) - 4m|^^ {pi ■ 92) (P2 • ^2) 

- "^l^lhiPl ■ (ll){P2 ■ qi) + '^irLMr^lhiPl ■P2), 

A^(.Pi,qj) = 8{pi ■ qi){p2 ■ q2){qi ■ 92) - 47nif^(pi • q2){p2 ■ ^2) 
+ ^iT^lhiPi ■qi){P2-qi)- '^iriMrrnlf^ipi ■ P2), 

C^(.Pi,qj) = 4(pi ■P2){qi ■q2f -A^{pi,qj), 

C^(.Pi,qj) = ■P2){qi ■ q2f - '^T^lh'mlj^ipi ■P2) - A'^{pi,qj), 



+ {pl 



P2) 



,(E.3) 
(E.4) 

(E.5) 

(E.6) 
(E.7) 
(E.8) 
(E.9) 
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Fm = 4:Gpflf flf \Vj 



MiJM2 \ ^Mi 



CKM^rCKM\2 
V 1 



\Ve^iVi^i\'^ml, 



M2 



ai,2{Pi,qj) = (qi - Pi,2f - ml; b = TNim4^. 
The decay rate for LV rare meson decay is then given by 



LV 



dPS: 



31 



,2 _2 _2 \\h^^2 _2 _2 \dm^i 



dyidy2dyzdyi 



cl 



(E.IO) 
(E.ll) 
(E.12) 



(E.13) 

(E.14) 
(E.15) 



where dPS^i and dPS^2 are the phase space factors obtained by conveniently clustering 
two different sets of particles to enable applying the narrow-width approximation easily. 
yi to ?/4 are rescaled angular variables with integration limits < < 1. The width 
of the heavy neutrino is very small compared to the mass and hence we can apply the 
narrow-width approximation. 



/ 



dm^. 



m 



m 



iy+r%mi 



Tn^-^ 



(5(m^. — ml)dm'l. 



(E.16) 



Applying the narrow-width approximation as described above and integrating over the 
(5— function we get 



/ f^'dPSs, 



( FmA^v)^^ 



N4 



(a2 + b^)AP(v) + b^BP(y) 



f^^''^dPSs2 



A 2 (niM^ , mj^ , 7714) A 2 (ml, m\,m\j^)j dyidy2dy2,dyi, (E.17) 

j f^^''UPS^i{pi^P2), (E.18) 
Now, we can find the decay rate from Eq. ( [E.13 ). The branching fraction is then given by 

The CKM matrix elements and the lifetimes of mesons used in our calculations are 
taken from PDG [99| and are listed below. 



iKfel = 0.00367, l^edl = 0.224, \V,, 



0.996; 



TK = 1.2384 X 10^*^ s, TD = 1.040 x 10"^^ ^ 4 9 ^ iQ-is s, tb = 1.671 x 10"^^ g 

The mass and decay constants of pseudoscalar and vector mesons used in our calculations 
are listed in Table |^. 
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Table 7 : Mass and decay constants of pseudoscalar and vector mesons used. 



Pseudoscalar 
Meson 


Mass 
(MeV)[|9l 


Decay Constant 
(MeV) @ 


Vector 
Meson 


Mass 
(MeV)@ 


Decay Constant 
(MeV)[m 


-r 


139.6 


130.7 


-r 


775.8 


220 




493.7 


159.8 




891.66 


217 




1869.4 


222.6 


m 


D* 


2010 


310 




1968.3 


266 




2112.1 


315 


B± 


5279 


190 [] 




UJ 


782.59 


195 




135 


130 




896.10 


217 


V 


547.8 


164.7 


117] 


<P 


1019.456 


229 




957.8 


152.9 


117] 




2006.7 


310 


Vc 


2979.6 


335.0 






3096.916 


459 [US] 
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